CATALOGED  '3V  4  03l  46 


403  146 

Vidya  Technical  Report  No.  63-95 
January  9, 1963 

THEORETICAL  INVESTIGATION  OF  DYNAMIC  STABILITY  DERIVATIVES 
OF  DUCTED  PROPaLERS 

by 

A.  R.  Kriebei 
A.  H.  Sack$ 

J.  N.  Nielsen 


Prepored  Under  Navy,  Bureau  of  Novol  Weopons 
Confroef  NOw  62-01 29-c 


Vidya  Project  No.  167/C 


f" 


V  I  D  Y  A 

The  word  Vidya,  taken  from  the  Vedanta  philosophy  of  the  Hindus, 
means  knowledge.  The  symbol  used  to  denote  the  Vidya  organization 
is  the  letter  from  Sanskrit,  the  ancient  language  of  India. 


I 

1 


VlOYA  REPORT  NO.  63-95 
January  9,  1963 


VlOYA  COPY  NO. 


THEORETICAL  INVESTIGATION  OF  DYNAMIC  STABILITY 
DERIVATIVES  OF  DUCTED  PROPELLERS 

by 

A.  R.  Kritb«l 
A.  H.  Socks 
J.  N.  Niolton 


Prepared  Under  Navy  Contract  NOw  62-0129-c 
Administered  Under  Technical  Direction  of  the  Bureau  of  Naval  Weapon 
Fluid  Mechanics  and  Flisht  Dynamics  Branch,  RRRE-4 


Vidya  Project  No.  167/C 


VIDYA 


A  DIVISION  OF 


Itek 


CORPORATION 


1450  PAGI  MILL  ROAD  •  PALO  ALTO,  CALIFORNIA 
TILi  DAvinfort  1-2455  TWXt  415  492-9270 


SUMMARY 


The  static  and  dynamic  stability  derivatives  of  a  ducted  pro¬ 
peller  are  predicted  by  theoretical  analysis.  The  analysis  is  based 
on  the  theory  presenced  in  Reference  1,  and  the  ducted  propeller  is 
initially  represented  by  a  short,  straight,  thin  ring  surrounding  a 
uniformly  loaded  actuator  disk.  Then  the  interference  effect  of  the 
ring  duct  on  an  actual  propeller  with  finite  number  of  blades  is 
calculated  in  order  to  determine  the  contribution  of  the  propeller 
to  the  stabilit/  derivatives  of  the  ducted  propeller. 

The  analysis  is  applied  to  the  case  of  a  torpedo-like  config¬ 
uration  consisting  of  a  body  of  revolution  with  a  rear-mounted  ducted 
propeller  in  axicti  flow.  The  natural  frequency  and  damping  of  the 
configuration  is  calculated,  and  the  individual  contributions  of 
the  duct,  the  propeller,  the  hull,  and  a  fin  (added  for  static  sta¬ 
bility)  are  determined.  It  is  found  that  the  ducted  propeller  is 
a  more  effective  oscillation  damper  than  is  a  fin  of  the  same  pro¬ 
jected  area  and  that  the  effectiveness  Increases  with  propeller  disk 
loading.  The  contribution  of  the  ducted  propeller  to  the  dynamic 
stability  of  the  torpedo  is  found  to  be  dominated  by  the  static  sta¬ 
bility  derivative  of  the  isolated  ducted  propeller. 

The  theory  for  the  ducted  propeller  is  compared  with  available 
experimental  da  a  in  covering,  steady  axial  flow,  and  steady  flow 
at  angle  of  attack.  It  is  found  that  the  important  derivative  Cjj 

CL 

is  predicted  well  for  small  angles  of  attack  by  the  present  theory. 


li 


The  division  of  thrust  between  duct  and  propeller  is  predicted  mod¬ 
erately  well  for  ducts  of  small  chord/diameter  ratios  with  no  flow 
separation.  Experimental  data  suitable  for  comparison  with  experi¬ 
ment  are  extremely  sparse  since  ficw  separation  is  generally  present 
at  the  duct  leading  edge. 

A  method  for  reducing  ducted  propeller  data  is  presented  and 
shown  to  give  better  agreement  with  theory  than  conventional  methods. 
A  qualitative  analysis  of  real  slipstream  effects  is  presented,  and 
an  experimental  program  is  outlined  for  obtaining  some  badly  needed 
data  on  ducted  propellers  without  leading-edge  flow  separation. 
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1 .  INTRODUCTION 


The  benefits  to  be  gained  by  ducting  or  shrouding  a  propeller 
for  static  operation  and  low>speed  flight  were  demonstrated  in  1931 
by  the  experiments  of  Stipa  in  Italy.  Since  that  time,  and  partic¬ 
ularly  in  recent  years  when  hovering  flight  and  vertical  take-off 
became  feasible  for  conventional  aircraft,  ducted  propellers  have 
received  considerable  attention  from  numerous  investigators.  A 
vast  number  of  individual  experimental  programs  have  been  under¬ 
taken,  as  reviewed  in  Reference  2,  and  a  number  of  VTOL  (vertical 
take-off  and  landing)  vehicles  were  designed  and  built  employing 
ducted  propellers. 

Because  of  its  complexity,  however,  the  theory  of  the  ducted 
propeller  has  lagged  f€u:  behind  experiment.  Only  now  is  a  rigorous 
self-consistent  theory  emerging,  even  for  predicting  the  aerodynamic 
forces  and  moments  on  ducted  propellers  in  steady  axial  flow.  This 
theory  (Refs.  3  and  4)  is  based  on  the  classical  method  of  singu¬ 
larities  in  three  dimensions  and  is  thus  far  restricted  to  small 
angles  of  attack  and  light  propeller  loadings  (i.e.,  small  pertur¬ 
bations).  However,  an  approximate  theory  for  all  angles  of  attack 
has  been  developed  by  Burggraf  in  Reference  1  in  which  the  ducted 
propeller  is  represented  by  a  short,  thin,  straight  duct  surrounding 
a  uniformly  loaded  actuator  disk.  This  representation  permits  the 
assumption  of  nearly  two-dimensional  flow  over  each  chordwise  strip 
of  the  duct. 

As  for  dynamic  stability  derivatives,  the  situation  is  far 
worse.  At  this  time,  there  appear  to  be  neither  esqperimental  data 


nor  any  theoretical  analysis  for  the  dynamic  stability  derivatives 
of  an  Isolated  ducted  propeller. 

The  purpose  of  the  present  investigation  is  to  provide  a  the¬ 
oretical  analysis  for  predicting  the  forces,  moments,  and  stability 
derivatives  of  an  isolated  ducted  propeller  for  hovering  and  for 
forward  flight  at  all  angles  of  attack.  In  order  to  accomplish 
this  end,  we  shall  first  introduce  the  assvunptions  made  by  Burggraf 
in  Reference  1.  We  shall  then  extend  his  analysis  to  cover  the 
cases  of  pitching  and  plunging  flight,  and,  finally,  we  shall  gen¬ 
eralize  the  analysis  to  include  the  interference  effect  of  the 
duct  on  a  propeller  with  finite  number  of  blades.  In  the  course 
of  the  analysis,  we  shall  investigate  the  approximations  made  in 
Reference  1,  particularly  the  averaging  processes  used  in  calculat¬ 
ing  certain  induced  velocities,  and  we  shall  investigate  the  effect 
of  altering  the  propeller  location  within  the  duct. 

The  theory  developed  for  the  ducted  propeller  will  be  applied 
herein  to  a  torpedo-like  configuration  in  order  to  investigate  its 
dynamic  stability  and  to  get  some  insight  into  the  effectiveness  of 
the  ducted  propeller  as  a  stabilizer  and  oscillation  deunper.  For 
this  problem,  the  thrust  interference  between  the  hull  and  the 
ducted  propeller  will  also  be  investigated  by  approximate  methods. 

In  addition  to  the  basic  analysis,  a  qualitative  analysis  of 
real  slipstream  effects  will  be  presented,  and  the  latter  will  be 
used  in  developing  a  rational  method  for  comparison  of  theoretical 
and  experimental  results. 
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2.  APPROACH  TO  THE  OUCTED'PROPELLER  ANALYSIS 

The  primary  objective  of  this  study  is  the  prediction  of  the 
stability  derivatives  of  a  ducted  propeller  without  restriction  on 
either  angle  of  attack  or  propeller  disk  loading.  This  is  done  by 
finding  the  duct  derivatives  in  the  presence  of  the  propeller  first 
and  then  the  propeller  derivatives  in  the  presence  of  the  duct. 

2.1  The  Duct  in  the  Presence  of  the  Propeller 

The  aerodynamic  forces  and  moments  on  the  duct  in  the  presence 
of  the  propeller  are  found  by  representing  the  propeller  as  a  uni¬ 
formly  loaded  actuator  disk.  The  net  duct  reaction  is  expressed 
as  a  thrust  force  '^^^(P)  ’  ®  normal  force  Nj^(p)»  ®  pitching 

moment  about  the  midchord  diameter  shown  in  Figure  1. 

For  a  given  set  of  design  variables  c/D  and  y/Vn  , 

^P(D) 

these  three  components  can  be  written  as  dimensionless  coefficients 
in  terms  of  the  flight  conditions  as  follows : 


gplpj  . 

(“o-  ^ 

Da  N 

•  V  ) 

O  ' 

(1) 

'^D(P)  « 

(°’')d(P) 

(a  ^ 

1  “o  *  V 
\  o 

Da  ^ 

’  V  J 
o 

(2) 

DAq 

o 

\  "*)d(p) 

(a  ^ 

\  o’ 

Da 

(3) 

Thus  there  are  nine  stability  derivatives  (with  respect  to  a^, 
q,  and  d)  for  the  duct  in  the  presence  of  the  propeller,  since 
is  assumed  constant.  The  way  in  which  these  are  defined  will  be 
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illustrated  now  for  the  three  moment  derivatives.  But  since  a 

is  not  restricted  to  smctll  values,  an  a  coefficient  will  be  used 

rather  than  an  a  derivative.  Thus,  partial  differentiation  of 

(C  )  r>/B\  with  V  constant  and  a  ■  a  ■  constant^,  and  then 
in  JLi  V  “  /  o  o 

setting  q  ■  d  ■  0,  gives  the  q  and  d  stability  derivatives: 


Hence,  for  small  rates  of  pitching  q  and  plunging  d,  the  in¬ 
stantaneous  duct  moment  is  given  by 


'"b{P)  "b(P) 


(a^,0,0) 


D(P' 


(6) 


The  first  term  is  the  static  moment  coefficient.  The  static  coeffi¬ 
cient  and  the  q  derivative  are  found  in  the  present  analysis  by 
considering  the  ducted  propeller  in  the  pure  pitch  flight  condition, 
and  the  d  derivative  is  estimated  independently  by  use  of  the 
apparent  mass  concept.  For  pure  pitch,  the  ducted  actuator  disk 
moves  as  though  it  were  mounted  rigidly  to  an  arm  with  fixed 

radius  R  as  shown  in  Figure  2.  The  arm  rotates  at  constant 

3  .  .  -  —  —  —  -  -  - 

^Note  that  a  differs  from  a  only  for  plunging  flight 
(See  Fig.  1) 


angular  velocity  q  .  Both  and  are  constant  and 

The  method  of  singularities  is  used  to  find  the  velocity  distribu¬ 
tion  W  over  the  duct  surface  as  induced  by  all  of  the  bound  and 
shed  vorticity  for  the  ducted  actuator  disk  shown  in  Figure  2 .  The 
moment  coefficient  for  this  motion 


is  then  found  by  integration  of  the  static  pressure  ^Pq  "*■  2^ ”  2^ 

over  the  duct  surface.  Then  by  setting  q  ■  0  (but  holding 

constant),  one  obtains  D(P)  0,  0)  .  By  differentiating 

(C„)  _/_v  with  respect  to  q  as  described  above,  one  obtains 
m  D(P) 

The  total  surface  velocity  distribution  W  on  the  duct 
during  pure  pitch  is  considered  to  be  composed  of  six  parts  as 
shown  in  Figure  2.  That  is, 

W“w.  +  w+w^  +  w,  +w  +w  (7) 

h  a  q  h^  q^ 

The  first  component  Wj^  is  induced  by  the  actuator  disk  for  no 
motion  of  the  duct  (i.e.,  for  hovering  flight).  The  second  two 
components  are  due  to  the  two  components  of  the  duct  motion:  a 
translational  motion  at  a  ,  which  produces  w  ,  and  a  rota- 

tional  q  about  the  midchord  duct  diameter  which  produces  w^.  The 
w^  and  Wq  components  include  the  velocity  induced  by  the  duct 
wake  which  is  generated  in  each  case. 
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It  is  assumed  that  the  propeller  thrust  coefficient  C_ 

P  (D) 

is  sufficiently  large  to  produce  a  large  jet  velocity  Vj  com¬ 
pared  with  the  lateral  component  of  flight  speed.  That  is, 

Vj  »  sin  (8) 

In  this  case,  the  jet  extends  downstream  nearly  aligned  with  the 
duct  axis  as  shown  in  Figure  2  for  any  combination  of  the  duct 
motions.  Furthermore,  the  duct  motions,  their  vorticity  distribu¬ 
tions,  and  the  reactions  which  they  produce  on  the  duct  are  all 
additive  (See  Eq.  6) .  The  propeller  vorticity  y  and  the  duct  wake 
vortex  filaments  dP^  both  lie  on  the  cylindrical  surface  r  “  D/2 
and  are  mutually  perpendicular  (Fig.  2) . 

The  duct  wake  vortex  fil2unents  produce  a  counter  rotating 
swirl  distribution  (Fig.  6) .  There  is  no  swirl  generated  directly 
by  the  propeller  in  the  present  analysis  since  it  is  represented 
by  an  ideal  actuator  disk.  However,  the  effect  of  real  propeller 
slipstream  effects,  including  swirl,  is  considered  qualitatively 
in  Section  4.2.2. 

The  bound  vorticity  distributions  on  the  duct  corresponding 
to  hovering , angle  of  attack,  and  pitching,  are  composed  of 
Yjj,  and  y^,  which  are  calculated  with  the  assumption  of  two-dimensional 
flow  over  local  chordwise  strips  of  the  duct  for  c  «  D.  But 
since  the  vorticity  is  actually  distributed  on  a  ring  rather  than 
a  flat  plate,  there  are  three  additional  self-induced  axial  vel¬ 
ocity  components  w.  ,  w  ,  and  w  .  These  are  analogous  to  the 

flj  ot. 


self-induced  axial  velocity  (or  propagation  speed)  of  a  ring  vortex 
Since  the  self-induced  radial  velocity  distribution  is  small  for 
00  >  D  »  c,  it  is  neglected.  In  order  to  make  the  analysis 
tractable,  it  is  found  necessary* to  approximate  the  three  self- 
induced  components  since  it  is  not  feasible  to  express  them  analyti 
cally  (See  Appendix  D)  and  to  incorporate  them  into  the  subsequent 
analysis.  The  approximation  cannot  be  made  by  simply  concentrating 
the  bound  vorticity  on  the  duct  into  a  single  ring  filament  of 
strength  ^  ^  such  a  ring  induces  an  infinitely 

large  axial  velocity  upon  itself.  The  method  employed  here  and  in 
Reference  1  is  to :  (1)  average  the  bound  vorticity  over  the  duct 

chord,  (2)  find  the  self-induced  velocity  distribution  due  to  the 
averaged  vorticity,  and  (3)  then  average  the  self-induced  velocity 
over  the  duct  chord  and  use  these  values,  w,  ,  w„  ,  and  w  ,  in  the 

hi 

subsequent  analysis.  The  accuracy  of  this  averaging  process  is 
assessed  in  Appendix  O. 

2.2  The  Propeller  in  the  Presence  of  the  Duct 

The  stability  derivatives  for  the  isolated  propeller  are  de¬ 
fined  just  as  for  the  duct  for  the  three  components  of  the  net 
aerodynamic  force  Tp,  Np,  Yp,  and  for  the  two  components  of  the 
moment  Mp  and  N*p  which  are  about  propeller  diameters.  The 
interference  of  the  duct  on  the  propeller  due  to  a  and  q  is 
found  to  produce  only  pitching  moments  ^(D-*P)q* 

The  thrust  interference  produced  by  the  duct  in  steady 

axial  flight  (i.e.,  for  ■  0)  is  not  predicted  here  since  it  is 
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conaidered  to  be  part  of  the  initial  design  process.  Thus  in  the 
presence  of  the  duct  each  component  of  the  propeller  reaction 
(  is  the  same  as  for  the  isolated  propeller  (  )p  except 

for  the  pitching  moment  which  is  given  by  Mp(D)  “  Mp  +  Mjj_^p. 

The  a  and  q  derivatives  are  obtained  from  Reference  5  for 
the  isolated  propeller  at  small  a^.  The  d  derivative  and  the 
duct  interference  will  be  estimated  in  the  present  analysis. 
Although  the  propeller  thrust  is  generally  greater  than  the  duct 
thrust,  the  flight- induced  forces  and  moments  on  the  propeller  in 
rho  presence  of  the  duct  are  found  to  be  small  compared  with  those 
for  the  duct  for  typical  conditions  except  for  .  There- 

fore,  a  more  refined  analysis  of  the  propeller  reactions  does  not 
appear  to  be  warranted  except  possibly  for  the  duct  interference 
moment  on  the  propeller  due  to  pitch,  .  This  moment  may 

be  appreciable  for  isolated  ducted  propellers,  such  as  in  VTOL 
applications,  but  even  it  is  negligible  when  the  ducted  propeller 
is  combined  with  a  long  hull,  as  for  a  torpedo. 
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3 .  ANALYSIS 

3.1  Stability  Derivatives  of  a  Duct  in  the  Presence  of  an 

Actuator  Disk 

3.1.1  Discussion 

The  purpose  here  is  to  determine  the  force  and  moment  coef¬ 
ficients  and  the  stability  derivatives  or.  a  th.1n  duct  of  uniform 
diameter  with  an  actuator  disk  located  in  the  exit  plane  (Fig.  1(a)) 
as  in  Reference  1.  (The  results  are  actually  valid  for  any  actuator 
disk  location  within  the  duct  as  shown  in  Appendix  C.)  The  coef¬ 
ficients  and  derivatives  here  pertain  to  the  aerodynamic  force  and 
moment  acting  on  the  duct  in  the  presence  of  the  actuator  disk. 

The  initial  assumptions  are: 

(a)  The  fluid  is  inviscid  and  incompressible,  and  a  velocity 
potential  represents  the  actual  flow  field  in  which  the  boundary 
layers  are  negligibly  thin  and  unseparated. 

(b)  The  chord-to-diameter  ratio  of  the  duct  c/D  is  small 
enough  so  that  the  flow  over  each  chordwise  strip  is  nearly  two 
dimensional. 

(c)  The  lateral  component  of  flight  speed  is  much  smaller 
than  the  slipstream  (or  jet)  velocity;  that  is,  Vj  »  sin 

in  Figure  2.  In  this  case  the  slipstream  extends  nearly  axially 
downstream  as  shown. 

(d)  The  propeller  is  represented  by  an  ideal  actuator  disk 
which  produces  a  uniform  loading  or  pressure  rise  Ap  over  the 
disk  but  no  swirl . 

(e)  The  diameter  of  the  slipstream  is  considered  constant 
and  equal  to  the  duct  diameter  for  the  purpose  of  calculating  the 
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duct  surface  velocity  W.  This  assumption  is  discussed  on  page  130 
of  Reference  6  and  in  Section  4.2.1. 

With  the  above  assumptions ,  the  method  of  singularities  is 
used  in  Appendix  A  to  determine  the  duct  surface  velocity  distribu¬ 
tion  W  for  the  pure-pitch  flight  condition  (Fig.  2) .  The  pres¬ 
sure  distribution  corresponding  to  W  is  integrated  over  the  duct 
surface  to  obtain  the  net  force  and  moment  on  the  duct  in  Appendix  B. 
In  the  present  section,  we  shall  find  the  resulting  force  and  moment 
coefficients  for  the  duct  in  the  presence  of  the  actuator  disk  both 
in  steady  flight  at  angle  of  attack  and  in  hovering.  Then  we  shall 
find  the  pitching  (q)  stability  derivatives  and  finally  the  plunging 
(d)  derivatives. 

3.1.2  Static  (a)  coefficients 

Let  us  find  the  static  (a)  coefficients  of  the  duct  by  the 
procedure  in  Section  2,  from  the  results  in  Appendix  B.  These 
coefficients  give  the  net  aerodyneunic  reaction  on  the  duct  when 
q  "  d  “  0,  that  is,  for  steady  flow  at  angle  of  attack  a^. 

There  is  a  normal  force  ^^(p)  ^  pitching  moment  ^**n)d(P) 

which  are  associated  with  the  nonsingular  portion  of  the  pressure 
distribution,  that  is,  excluding  leading-edge  suction.  The  latter, 
in  turn,  gives  rise  to  a  thrust  force  ^^(p)  ®  pitching  moment 

C*^}d(P) *  shall  now  consider  these  four  components  in  succession. 
The  general  expression  for 


('=h)d(P)  (“o'  ^  •  “) 
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corresponding  to  Equation  (1)  is  given  by  Equation  (B.33)  of 
Appendix  B.  When  q  *  0,  we  obtain  Equation  (B.36) .  By  sub¬ 
stituting  Equations  (A. 53)  and  (A. 77)  for  w.  and  w  and  by 

"i  “i 

dropping  terms  higher  than  second  order  in  c/D,  we  find  the 
static  coefficient 

(<=h)d(p)  r" 

For  the  special  case  of  a  ring  wing  at  small  angle  of  attack  (i.e., 
for  7  ■  0  and  «  1) ,  this  result  reduces  to  that  of  Refer¬ 


ence  7. 


The  general  expression  for  ^**u}d(P)  given  by  Equation  (B.35) 

and  when  q  •  0  by  Equation  (B.37) .  Substituting  Equations  (A. 53) 

and  (A. 77)  for  w.  and  w  and  then  retaining  only  the  lowest 
"i  “i 

order  terms  in  c/D,  we  find 

/c  'N  (a  ,0,0)  « -S—  (  +  2w  cos  a_)  sin  a  (10) 

The  duct  thrust  Td(p)  expressed  by  Equation  (B.41) ,  and 

when  q  ■  0  we  obtain  the  thrust  coefficient 


(‘=t)d(p)  -  ^(™)  (*"  ^  -  “) 


^  2D 
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The  expression  for  (Mt)j3(p)  is  given  in  Equation  (B.42), 
and  when  q  -  0,  we  find 


(a^,0,0) 


4c7/DV^ 

TTC 


1  + 


2D 


(vn  ^  -  2) 


Sin  a 


(12) 


Thus,  except  for  the  case  of  very  low  disk  loading  (7/V^  -♦0),  we 
find  that  C  in  Equation  (10)  is  much  smaller  (higher  order  in 

c/D)  than  Cj^  in  Equation  (12).  Therefore,  we  shall  take  the 
total  pitching  moment  on  the  duct  to  be  equal  to  the  component 
arising  from  leading -edge  suction,  so  that 
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)  «c 
(13) 


3.1.3  Thrust  ratio  in  hovering  flight 

For  hovering  flight,  the  duct  thrust  is  given  by  Equation 
(B.41)  with  ■  0.  That  is, 

O  ' 


''D(P)h 


CD 

2 


£ll 


(in  ip 


(14) 


Also,  for  hovering  flight,  application  of  Bernoulli's  equation 
far  ahead  and  far  behind  the  propeller  disk  gives 

7  -  Vj  -  V  2Ap/p 


(15) 
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Thus,  the  actuator -disk  loading  is  given  by 


A 


Ap 


(16) 


and  we  find  that  the  thrust  ratio  for  hovering  flight  is  therefore 


'^P(D)h 


16D 

c 


(17) 


3.1.4  Pitching  (q)  derivatives 

We  shall  now  find  the  q  derivatives  of  the  duct  by  using 
the  procedure  of  Section  2  with  the  results  of  Appendix  B.  Again, 
we  shall  consider  in  succession  the  q  stability  derivatives  of 
each  of  the  four  duct  reaction  components 

**D(P)’  (“n)d(P)*  '^D(P)’  (*S’)d(P) 

Thus,  by  differentiating  Equation  (B.33)  with  respect  to 

Dq/V^,  holding  and  constant,  we  obtain  Equation  (B.38) 

for  ( C  ]  .  But  if  only  the  two  lowest  order  terms  in  c/D 

\  “q/D(P) 

are  retained.  Equation  (B.38)  reduces  to 


.  DIS- 


16D 

c 


1^  +  cos  a 


(18) 


By  the  same  procedure.  Equation  (B.35)  gives  (B.39)  and 
we  get 


■  'VC 

2VD^ 

o 


(19) 
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Similarly,  if  Equation  (fi.41)  is  differentiated  with  respect 
to  holding  and  constant,  and  if  q  is  then  set 
equal  to  zero,  we  obtain 


(%)i 


2  3 

IT  c  Sin  a 


"■>>  (l  ^  II)' 

Differentiating  Equation  (B.42)  in  the  same  way  gives 

Tryc^  (in  ^  -  2) 


2v_D'’  r  1  + 


0 


TTC 


2d} 


(20) 


(21) 


Now,  since 


is  of  higher  order  in  c/D  than 
/D(P)  \  ■"q/D(P) 

retention  of  terms  only  through  second  order  in  c/D  gives  for 


C'"»q)D( 


the  q  derivative  of  the  pitching  moment 


yc 


(P) 


2V  D 
o 


(22) 


3.1.5  Plunging  (a)  derivatives 

In  order  to  estimate  the  plunging  (d)  derivatives  of  the  duct, 
we  shall  use  an  approximate  analysis  which  is  independent  of  the 
foregoing  results  for  a  and  q  and  which  is  based  on  the  appar¬ 
ent  mass  concept ,  as  mentioned  in  Section  2 .  Again ,  we  shall  make 
use  of  the  initial  assumption  that  c  «  D  and  utilize  the  flat- 
plate  formulas  for  each  chordwise  strip  of  the  duct. 

We  shall  estimate  the  d  derivatives  of  the  duct  normal  force 
and  the  duct  pitching  moment,  but  not  of  the  thrust  force.  The 
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reason  for  this  is  that  only  the  normal  force  and  moment  deriva¬ 
tives  appear  independently  in  the  equations  for  the  natural 
frequency  and  damping  of  a  vehicle  at  small  angle  of  attack. 

The  approximate  method  used  to  find  the  duct  d  derivatives  seems 
justified  by  the  fact  that  they  are  found  to  represent  only  about 
1  percent  of  the  total  damping  of  a  torpedo-like  configuration, 
to  which  the  static  normal-force  derivative  C„  of  the  duct  is 
actually  the  major  contributor. 

The  plunging  stability  derivatives  of  the  isolated  duct  will 
be  calculated  by  considering  the  flight  condition  illustrated  in 
Figure  1(a)  where  the  duct  is  situated  in  a  steady  free  stream  at 
small  angle  of  attack  (a^  «  1)  and  oscillates  with  a  normal 
velocity  conponent 

V  ■  v„  sin  vt  (23) 

m 

and  a  simultaneous  axial  velocity  component 

u  ■  u„  sin  vt  (24) 

m 

The  oscillations  are  assumed  to  be  of  small  amplitude  and  of  low 
frequency.  In  addition,  the  oscillations  are  taken  to  be  normal 
to  with  u  and  v  in  phase  with  one  another.  That  is, 

u  ■  -  V  sin  (25) 

so  that  V  -  0,  and  d  is  given  by 

«  • 

_ y _  a  _  .JL 

cos 

o  o  o 


a  » 


(26) 
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The  instantaneous  flow  condition  corresponds  to  the  free 
stream  V  at  a  as  shown  in  Figure  1(a).  The  distribution  of 
bound  and  trailing  vorticity  from  the  duct  now  consists  of  two 
parts : 

(a)  The  steady-state  distribution  corresponding  to  the 
instantaneous  flow  condition  (V  at  a  in  Section  3.1.2)  which 
is  in  phase  with  velocity  v. 

(b)  A  distribution  which  is  in  phase  with  acceleration  v. 

Within  the  framework  of  this  analysis,  the  force  and  moment 
acting  on  each  chordwise  strip  of  the  duct  correspond  with  those 
on  a  flat  plate  in  the  same  orientation  performing  the  same  motion. 
Now,  the  force  and  moment  on  an  oscillating  flat  plate  in  a  uniform 
stream  at  small  angle  of  attack  can  be  found  in  Reference  8 
(pp.  293-303),  and  it  can  be  shown  that  for  any  combination  of 
u^  and  at  low  amplitude  and  frequency,  there  are:  (1)  a 

lift  force  acting  at  the  quarter  chord  which  corresponds  with  the 
instantaneous  flow  condition  (V  at  a)  and  which  is  in  phase 
with  the  velocity  v^,  and  (2)  a  force  normal  to  the  plate  acting 
at  its  midchord  and  in  phase  with  acceleration  v^.  At  low  fre¬ 
quency  the  latter  force  per  unit  span  is  approximately  given  by 
the  quasi-steady  formula 

Nf  -  -  ^  pVf  (27) 

and  evidently  it  does  not  depend  upon  either  the  free-stream  condi¬ 
tion  or  the  axial  component  of  acceleration  (V,  a,  or  u^) .  The 
factor  (ircV4)p  is  recognized  as  the  apparent  mass  of  a  flat 
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p]ate  of  chord  c  undergoing  normal  oscillations  v^.  For  stream- 
wise  oscillations  u^,  the  apparent  mass  is  clearly  zero. 

Applying  the  above  results  to  each  chordwise  strip  of  the 
isolated  duct  and  integrating  round  the  duct ,  one  finds  a  normal 
force  acting  at  the  midchord  diameter,  in  phase  with  v,  which  is 
given  by 


N, 


8 


pv 


TT^C^DV^dp 

8 


(28) 


Thus,  the  apparent  mass  of  the  duct  (for  ❖  at  small  a^)  is 
that  of  a  toroidal  volume  of  fluid  with  the  dimensions  shown  in 
Sketch  A. 


Sketch  A.-  Apparent  mass  of  a  duct  for 
lateral  oscillation  (v) . 
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Now,  as  described  in  Section  2,  the  plunging  derivative  can  be 
found  by  differentiating  Equation  (28)  with  respect  to  D^/Vq 
holding  constant.  This  gives,  for  «  1, 


(29) 


It  will  be  assumed  further  that  the  propeller  interference  has  no 
effect  on  N_,  so  that  when  a  «  1, 

Urn  O 


(30) 


This  seems  to  be  a  reasonable  approximation  since  the  apparent 
mass  of  the  duct  is  independent  of  the  free  stream  (Eq.  (27)), 
and  since  the  corresponding  apparent  mass  of  the  propeller  is 
much  smaller  than  that  of  the  duct.  Also,  since  acts  at 


the  midchord  diameter,  we  find  that  the  corresponding  moment 
vanishes.  That  is. 


(31) 


3.2  Stability  Derivatives  of  an  Isolated  Propeller 

The  duct  stability  derivatives  have  been  calculated  in 
Section  3.1  by  representing  the  propeller  as  a  uniformly  loaded 
actuator  disk,  for  which  the  stability  derivatives  are  all  zero. 
We  expect  the  actual  propeller  derivatives  to  be  small  compared 
with  those  of  the  duct,  and  this  will  be  verified  here  (with  the 
exception  of  )  for  a  typical  example.  To  calculate  the 
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propeller  stability  derivatives,  we  shall  first  use  the  analysis 
of  Reference  5  to  calculate  the  derivatives  of  a  simple  paddle-wheel 
type  of  isolated  propeller  (having  constant  blade  pitch  and  chord) . 
Then  the  duct  interference  on  the  propeller  will  be  estimated  by 
approximating  the  change  in  flow  due  to  the  duct.  These  effects 
will  be  added,  and  finally  the  resulting  propeller  derivatives 
will  be  compared  with  those  of  the  duct. 

The  analysis  of  Reference  5  is  based  upon  blade  element  theory 
in  which  each  element  operates  in  a  uniform  induced  flow  field 
which  is  found  from  simple  momentum  theory.  Empirical  constants 
are  introduced  such  as  a  "side  wash  factor"  and  a  "spinner  factor," 
and  the  theory  is  limited  to  small  angle  of  attack  a^.  The  theory 
gives  the  a  derivatives  of  either  a  single  or  a  dual  propeller, 
and  the  q  derivatives  of  a  single  propeller. 

3.2.1  Static  (a)  derivatives 

The  static  derivatives  of  a  single  isolated  propeller  at  a 
small  constant  angle  of  attack  a^,  (p.  11,  Ref.  5)  give  rise  to 
a  normal  force  Np  and  a  yawing  moment  Np*.  For  dual  counter- 

CL  Gt 

rotating  propellers,  Np  is  increased  from  18  to  32  percent  for 

a 

typical  aircraft  propellers,  and  Np*  is  zero.  Specific  values 
have  been  calculated  for  the  following  assumed  conditions: 
c,  ■  blade  section  lift  curve  slope  ■  2ir 

3  ■  blade  pitch  angle  ■  const 


b  •  blade  chord  ■  const 
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V 

■  advance  ratio  ■  0.1 

a  -  4^  ■  propeller  solidity  (at  3/4  D)  *  0.1 
^  (where  B  -  number  of  blades) 

C„  ■  propeller  thrust  coefficient  ■  2.17 

The  following  factors,  defined  in  Reference  5,  are  used: 

a  "  function  of  T  ■  0.392  for  T  ■  J  C_  -  0.85 

c  c  o  1  p 

f (a)  “  function  of  a  ■  1.52 
k  ■  side-wash  factor  ■  0.4 

cl 

kg  ■  spinner  factor  ■  1.14 

3 

“  2  V  sin  p 

For  three  values  of  p,  the  static  stability  derivatives  are  given 
in  Table  I  for  either  single  or  dual  propellers  (whichever  gives 
the  larger  value  of  the  derivative  in  question) . 

The  approximate  values  of  ^  in  Table  I  are  from  the 

following  simple  equation  (p.  15,  Ref.  5)  which  neglects  the 
effects  of  induced  velocity 

^  »  f(a)kgal^  -  0.816  sin  p  (32) 

It  can  be  seen  in  Table  I  that  this  equation  gives  approximate 
values  about  three  times  greater  than  the  more  refined  analysis 
for  the  present  case. 

3.2.2  Pitching  (q)  derivatives 

Equations  for  the  q  derivatives  of  a  single  isolated  pro¬ 
peller  are  given  on  page  13  of  Reference  5.  The  pitching  propeller 
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develops  a  side  force  Yp  and  a  pitching  moment  Mp,  and  we  obtain 
the  values  in  Table  II  for  the  same  conditions  assumed  in  Section  3.2.1 
above . 

The  approximate  values  in  Table  II  are  from  Equations  (49) , 
Reference  5,  which  neglect  the  effects  of  induced  velocity.  These 
equations  can  be  written  as 


■  -  ^7r(l  +  a)  a  cos  3  -  -  0.33  cos  3 


(S). 


-  2(1  +  a)a  =  -  0.28 


(33) 

(34) 


where  a  •*  0.85  and  a  “  0.1,  as  before.  It  can  be  seen  in  Table  II 
that  these  equations  give  approximate  values  about  twice  as  large 
as  those  from  the  more  refined  analysis  for  the  present  case. 

For  comparison  with  Equation  (34) ,  the  simple  blade  element 


method  (used  in  Ref.  9)  was  employed  to  estimate 


for  a 


paddle-wheel  type  of  propeller  O  ■  const,  b  »  const,  c 


I, 


2v)  . 


This  method  assumes  a  very  low  advance  ratio  (oiD  »  V^)  ,  neglects 


induced  flow  effects,  and  gives 


32 


TTO 


(?) 


(35) 


For  the  previously  assumed  conditions,  we  have 

1 


(?) 


(36) 


so  that  Equation  (35)  gives 


-  0.29 


(37) 
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This  result  is  in  good  agreement  with  Equation  (34) .  This  latter 
method,  however,  gives  no  side  force. 

3.2.3  Plunging  (d)  derivatives 

The  (d)  derivative  of  a  plunging  isolated  propeller  will  be 
estimated  here  by  use  of  the  apparent  mass  method,  applying  the 
flat-plate  result  given  by  Equation  (27)  to  each  blade  element. 
Consider  one  pair  of  propeller  blades,  shown  in  Sketch  B,  having 
constant  chord  b,  constant  pitch  P,  and  located  at  azimuth 
angle  0.  If  the  propeller  oscillates  laterally  (v  at  low  fre¬ 
quency  and  amplitude)  then  the  component  of  v  normal  to  the 
blade  surface  is  ^  sin  p  cos  0  and  there  is  a  force  in  phase 
with  V  acting  at  the  center  of  each  blade  surface  and  normal  to 
it  (Eq.  (27)).  Thus,  the  total  force  on  each  blade  is 

^  pv  sin  p  sin  0 


Sketch  B.-  Oscillating  propeller  nomenclature. 
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Hence,  the  net  force  on  the  blade  pair  normal  to  the  duct  axis  is 
simply  (Sketch  B) 

N  -  2N^  sin  3  sin  (^  ■  -  -  pv  sin®  3  sin®  0  (39) 

For  B/2  pairs  of  blades,  we  have  the  total  normal  force  on  the 
propeller  (integrated  over  0  to  give  the  time  average  value) 

Np  -  ttN  j  “  -  (sin®  P)pv  (40) 

When  «  1  such  that  d  =  -  v/V^,  differentiating  with  respect 
to  Dd/V^  (holding  constant)  gives 

sin=  e  (41) 

3.3  Interference  of  Duct  on  Propeller 

Here  we  shall  estimate  the  aerodynamic  force  generated  by  the 
duct  on  the  propeller  blades.  We  shall  approximate  this  effect 
by  first  estimating  the  velocity  profiles  induced  by  the  duct  at 
the  propeller  plane  and  then  estimating  the  effect  produced  on 
the  propeller  blade  elements.  For  configurations  which  include  a 
long  hull,  the  total  propeller  reactions  including  interference 
will  be  relatively  unimportant  (see  Section  3.5).  For  the  isolated 
ducted  propeller  (Section  3.4),  the  propeller  pitching  moment  due 
to  pitch  is  the  only  interference  effect  which  is  not 

small  compared  with  the  duct  reaction  itself.  For  configurations 
where  this  is  significant  (e.g.,  flying  platforms  and  VTOL  designs), 
the  present  estimate  of 


should  be  considered  a  first 
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approximation  to  the  first  step  of  an  iterative  process.  That  is, 
duct  interference  causes  the  propeller  disk  loading  to  be  non- 
uniform,  which  changes  the  bound  vorticity  and  net  reaction  on 
the  duct.  This,  in  turn,  changes  the  propeller  disk  loading,  and 
so  on.  This  process  is  beyond  the  scope  of  the  present  study. 

Here  the  duct  interference  is  considered  to  be  due  only  to 
the  axial  velocity  component  induced  by  the  duct  at  the  propeller 
location.  The  effects  of  the  radial  and  tangential  velocity  com¬ 
ponents  induced  by  the  duct  are  probably  small  and  are  not  considered. 
The  duct  interference  on  the  propeller  is  due  to  three  sources  which 
will  be  estimated  in  succession  below;  .  (1)  the  bound  vorticity 
(2)  the  bound  vorticity  due  to  the  free  stream  at  angle 

of  attack  a^,  and  (3)  the  bound  vorticity  due  to  the  pitching 

motion.  It  can  be  seen  in  Figure  2  that  the  trailing  vortex  fila¬ 
ments  shed  from  the  duct  do  not  directly  induce  an  axial  component 
of  velocity  at  the  propeller,  since  they  are  aligned  with  the  duct 
axis.  However,  they  do  induce  an  upwash  on  the  duct  and  thereby 
influence  the  bound  vorticity  and  y^  which  in  turn  causes 

duct  interference  upon  the  propeller. 

Interference  effects  for  the  plunging  (d)  motion  cannot  be 
found  with  the  apparent  mass  method  employed  here.  However,  with 
a  long  hull  (Section  3.5),  the  apparent  mass  of  the  ducted  propeller 
is  found  to  give  a  negligibly  small  force  and  moment  on  the  hull 
(Tables  V  and  VI) .  For  lateral  (v)  oscillations  of  an  isolated 
ducted  propeller,  the  apparent  mass  of  the  propeller  is  much  smaller 
than  that  of  the  duct  (Section  3.4).  Therefore,  interference  is 
probably  unimportant  here  also. 
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3.3.1  Velocity  profiles  induced  by  Yjj  +  Y 

ty 

within  the  duct  by  Yjj  +  Y*  We  have  already  found  the  axial  velocity 


Let  us  consider  the  axial  velocity  profiles  (u  +  u  )  induced 

Yh  Y 


at  the  inner  duct  surface;  it  is  ''h  ^  '^h  ’  where  is  given  by 

Equations  (A. 15)  plus  (A. 47)  and  w,  is  given  by  Equations  (D.5) 

"i 

through  (D.7),  evaluated  on  page  D-4. 

At  points  away  from  the  duct  surface,  the  elementary  velocity 

du  induced  by  the  bound  vorticity  on  each  ring  element  of  the 
'h 

duct  is  given  by  the  results  on  page  305  of  Reference  6.  The  vor¬ 
ticity  on  each  ring  element  can  be  found  from  Equation  (A. 50)  in 
terms  of  y.  Combining  these  results  and  summing  over  the  duct 
surface,  we  obtain  the  following  expression  for  the  total  velocity 
induced  (at  r,  6)  by  the  duct  bound  vorticity  (at  S'). 


u 


'>'h 


_cl_ 

2Tr^D 


K(k)  - 


1  + 


ZXr/K  -  1) 


(t-;j 


E(k)  >  dO' 


(42) 


where 


si  ( 


cos^  6-2  cos  6  cos  6'  +  cos^  6') 


(43) 


k^  - 


4(D^/c^)r/R 


cos®  6-2  cos  6  cos  6*  +  cos®  6'  + 


5(s-') 


(44) 
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For  various  radial  locations  (r  <  0.8R),  Equation  (42)  was  evaluated 
on  the  Vidya  IBM  1620. 

The  ratio  u^/y  for  the  axial  velocity  component  u^  induced 

by  a  vortex  cylinder  of  strength  y  can  be  obtained  directly  from 

Table  15  on  page  328  of  Reference  6.  The  combined  dimensionless 

axial  velocity  profiles  (u  +  u  )/y  induced  by  +  Y  as  found 

Yh  Y  n 

by  the  procedure  described  above  are  plotted  in  Figure  3(a)  for 
three  axial  stations  in  a  duct  with  c  ~  0.2D  and  with  the  pro¬ 
peller  in  the  exit  plane.  These  profiles  do  not  depend  upon  the 
location  of  the  propeller  in  the  duct  if  it  is  represented  as  an 
actuator  disk  (see  Appendix  C) .  Therefore,  the  profiles  show  the 

individual  velocity  components  u  and  u  ■  7/2  at  three  alter- 

Yh  Y 

nate  propeller  locations  within  the  duct. 

The  results  of  this  section  (Fig.  3(a))  show  that  and  y 

induce  a  nearly  constant  axial  velocity  at  the  propeller  which 

reduces  the  effective  angle  of  attack  of  the  free  stream  (a^  as 

seen  by  the  propeller)  and  causes  the  blade  pitch  p  (or  speed  cd) 

to  be  higher  for  a  fixed  thrust  coefficient  C_  .  Only  the  effect 

^P(D) 

of  yj^  is  directly  due  to  the  duct  since  y  is  shed  by  the  propeller. 
However,  the  effect  of  both  yj^  and  y  are  considered  here  to  be  a 
part  of  the  initial  propeller  design  problem  rather  than  thrust 
force  interference. 


3.3.2  Interference  produced  by  y 


Let  us  now  consider  the  axial  velocity  u  induced  by  y 


The  axial  velocity  induced  at  the  inner  duct  surface  by  y 


a 


a 
can 


be  found  from  w  +  w 
a  a- 


Thus,  Equations  (A. 66)  and  (D.12)  give 
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(.Xl 


.  EsL.  f  2S. 


e 


2TrD  V  c  2  *  in  ^-  gJ 


(45) 


It  can  be  seen  (Eq.  (A. 64))  that  r  is  proportional  to  cos  <p. 

Now,  the  radial  variation  of  u..  will  be  estimated  near  the 

T'a 

center  line  of  the  duct  by  considering  to  be  concentrated  at 

a  single  ring  vortex.  For  small  values  of  c/D  and  r/R,  this 

approximation  yields  accurate  values  of  u  .  We  can  therefore 

'a 

write  (using  pp.  305  to  307  of  Ref.  6) 

aw 


u 


2wD 


f  _£e: 

i  Au* 


■  1  -  cos(0  -  0')]  -t  1  -  R'}  cos  d0' 

Va 


+  R'*  -  2R'  +  1  +  2R'  [l  -  cos(0  -  0')j 


(46) 


where  R'  ■  r/R. 

Integration  of  Equation  (46),  by  use  of  Equation  (A. 54)  and 
the  results  on  page  306,  Reference  6,  gives 


VL 


2irDV^ 


a _ 

\fR'  k 


(3k*  -  4)K  +  (4  -  k*)E 


. 


^  Al  -k* 


E  -  2K 


(47) 


In  the  plane  of  the  vortex  ring  (^p  *  0) ,  we  obtain 


ML 


(R*  +  1)‘ 


(48) 
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and 


2XR 


-  4)K  +  (4  -  k®)E  +  2 


R‘^  +  1 
R'®  -  1 


E 


2 


R'  -  1 
R'  +  1 


(49) 


Thus*  evaluation  of  Equation  (4,9)  gives  the  following  variation  of 


(2TrD/r  )u  over  the  propeller  disk. 


a  Y 


a 


These  values  and  those  from  Equation  (45)  for  c  >■  0.2D  are 

plotted  in  Figure  3(b).  These  profiles,  when  added  to  those  in 

Figure  3(a)  give  the  duct-induced  velocity  distribution  over  the 

propeller  disk  at  angle  of  attack. 

We  shall  now  estimate  the  propeller  interference  produced  by 

u  by  making  the  rough  approximation  that  the  velocity  profiles 
'a 

are  linear  with  r/R,  as  shown  by  the  dashed  lines  in  Figure  3  (b)  . 
We  shall  see  later  that  the  final  results  are  insensitive  to  this 
approximation.  Thus,  the  assumed  velocity  profile  is  from  Equation 


£]S_ 


2irD 


COS  0 


(50) 


(45) 
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We  shall  also  assume  that  the  propeller  consists  of  B  blades  of 
constant  chord  b  and  that  its  peripheral  speed  cur  is  much 


larger  than  the  axial  velocity  u 


The  angle  of  attack  on  a 


'a 


blade  element  changes  by  an  amount  (“U  /use) ,  which  causes  a 


'a 


pitching  moment  (but  no  net  force)  on  each  pair  of  blades.  The 
elementary  pitching  moment  on  a  pair  of  blades  is  given  by 


dM 


(D— P) 


-  2b 


a 


^  (cur)^r  cos  dr 


(51) 


We  shall  assume  here  that  c 


2ir  as  before.  Then 


a 


dM 


(D-»P) 


-b  £ 


a 


(cos'  <,)  tan  I  t  in  ife 


(52) 


where  B  corresponds  to  the  axial  location  of  the  propeller  in 
the  duct.  Integration  gives  the  moment  per  pair  of  blades  as 


M', 


(D-P) 


a 


bpeuD^r, 

— 22 -  (cos^  (^) 


0? 


e  ^  ,  16D 

tan  —  +  In  - 

2  c 


(53) 


Thus,  for  B/2  pairs  of  blades,  the  total  pitching  moment  (time 
averaged)  on  the  propeller  due  to  u  is 


'a 


1 

“(D-P)^  “  2  17  j  *^(D-P)^ 


d()> 


128 


(54) 
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Dif ferentiating  with  respect  to  and  using  Equation  (A. 64) 

for  we  find  that  the  interference  of  the  duct  on  the  pro¬ 

peller  produces  the  following  interference  pitching  moment. 


3 


32 


-ff) 


tan  f-  +  In 

2  C 


-  c  ^ 

V„ 


cos  a. 


(55) 


If,  in  lieu  of  Equation  (50),  we  had  assumed  that  u  -  ^u  =  const, 
integration  over  r  would  have  given  a  result  only  4/3  times  larger. 
Thus,  the  interference  effect  is  quite  insensitive  to  the  profile 
shape.  For  small  angles  of  attack  (a^  «  1),  we  find  that  the 
ratio  of  interference  to  duct  moment  (Eq.  (13)  and  (55))  is  given  by 


^a): 


.  32^  Xo  2  g  tan  f  ^  in  ^ 


'd(p) 


128 


cuD\ 

^oJ 


in  ^-2 


(56) 


where  =  2,  3,  or  4  for  the  propeller  located  at  the  duct  lead¬ 
ing  edge,  midchord,  or  trailing  edge,  respectively  (see  Eq.  (D.12)). 

Thus,  (^C  is  probably  generally  small  compared  with  fc 

V  "'a^^P  \  "‘a/D(P) 

except  for  propeller  locations  near  the  duct  leading  edge. 

3.3.3  Interference  produced  by 

The  interference  from  the  duct  bound  vorticity  y^  will  be 
estimated  now  by  the  same  technique  used  above  for  y^.  At  the 
inner  duct  surface,  the  axial  velocity  induced  by  y^  is 
^u  ]  *  w  +  w  .  As  above,  we  shall  assume  that  the  axial 
velocity  profile  is  linear  with  radius,  so  that,  using  Equations 


(A. 88),  (A. 90),  (A. 94),  (D.16),  (D.17),  and  (D.18),  we  find 
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where  is  evaluated  on  page  D-7 . 

By  analogy  with  Equations  (50)  and  (54) ,  it  can  be  seen  that 
the  pitching  moment  on  the  propeller  due  to  vu  is 
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(D— P), 


-  igi)  4 


JL 


1  + 


—  tan  ^ 
TTC  2 

2D 


+  G, 


(58) 


or 


JL 


1  + 


—  tan  T  + 

TTC  2  2 

20 

(59) 


Thus,  using  Equation  (22)  for  the  duct  moment  gives  the  following 
moment  ratio  for  pitching: 


+  4 


JL. 


1  + 


TTC 

2D 


(60) 


In  Section  3.4.2  it  will  be  shown  that  this  ratio  is  large  for  the 
conditions  assumed  there.  The  tan (6/2)  term  in  Equations  (60) 
and  (57)  is  comparatively  small  except  near  the  leading  edge  of  the 
duct  where  it  dominates.  This  causes  the  interference  moment  to 
reverse  sign  with  variation  of  propeller  location  0. 
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Both  7  and  7  and  the  axial  velocity  profiles  induced  by 

Ur  T. 

them  are  all  proportional  to  cos  0.  This  causes  only  pitching 
moments  on  the  propeller  due  to  a  or  q  interference  so  that 
Np(D)  “  Np  and  “  Yp.  The  theoretical  pitching  moments  are 

infinitely  large  when  the  propeller  is  at  the  duct  leading  edge 
because  of  the  mathematical  singularity  there  which  invalidates  the 
assumption  that  the  angle  of  attack  on  the  propeller  blade  element 
changes  by  only  a  small  amount  due  to  duct  interference. 

3.4  Comparison  of  Duct  and  Propeller  Reactions  for  an  Isolated 
Ducted  Propeller 

The  static  coefficients  and  stability  derivatives  of  the  duct 
and  propeller  as  determined  above  will  now  be  compared  for  a  spe¬ 
cific  configuration  in  order  to  illustrate  their  relative  magnitudes 
and  overall  importance  for  an  isolated  ducted  propeller  (and  hence 
for  certain  VTOL  configurations).  It  is  found  that  all  of  the 
propeller  reactions  are  small  compared  with  the  duct  reactions 
except  the  propeller  pitching  moment  due  to  pitching.  In  Section 
3.5.4,  it  is  shown  that  even  this  component  is  negligible  when  the 
ducted  propeller  is  combined  with  a  long  hull.  It  should  be  noted, 
however,  that  two  of  the  reaction  components  for  an  isolated  propel¬ 
ler  are  perpendicular  to  the  motion  (a  yawing  moment  due  to 
and  a  side  force  due  to  q) .  The  results  given  below  are  for  small 
angles  of  attack,  an  aircraft  type  propeller  which  is  moderately 
loaded,  and  a  duct  with  c  ■  0.2D.  Duct  interference  on  the  pro¬ 
peller  causes  only  pitching  moments  with  a  and  0  (Section  3.3). 
The  duct  moments  are  about  its  midchord  diameter  and  the  propeller 
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moments  are  about  its  own  diauneter.^  The  propeller  conditions 

assumed  for  Section  3.2  are  also  applied  here;  for  exaunple, 

V 

0  “  1  and  C_  “2.17.  From  Equation  (170),  we  can  evaluate 

^  Tp (D) 

v/V  in  terms  of  C„  to  find  that  these  conditions  lead  to 

°  ^P(D) 

y/V  “  0.78. 

*  o 


3.4.1  Static  (a)  derivatives 


First,  consider  the  normal  force  derivatives  of  the  duct  and 
propeller  at  small  angle  of  attack  a^.  For  the  duct  in  the 
presence  of  an  actuator  disk  (Eq.  (9)  with  c  “  0.2D) ,  we  have 


/C„  \  -  1.91  +  0.845,'^  “  2.57  (61) 

V  D(P)  '^o 


For  the  propeller  (from  Table  I) 


C„  \  “  0.05,  0.10,  0.16  for  p  -  10°,  20°,  30°  (62) 

V  P(D) 


Thus,  the  normal  force  on  the  duct  is  much  larger  than  on  the 
propeller. 

Now,  let  us  compare  the  pitching  moments  on  the  duct  and 
propeller  and  the  yawing  moment  on  a  single  propeller-  The  duct 
moment  derivative  is  (Eq.  (13)). 

(c  \  -  1.45  ^  “  1.13  (63) 

V  D(P)  '^o 


^For  small  values  of  c/D  we  shall  assume  that  the  axes  of 
pitch  and  moment  for  the  duct  and  propeller  coincide. 


From  Table  I ,  the  propeller  yawing  moment  derivative  is 


*N  "0.4  for  10  <  P  <  30°  for  a  single  propeller 

\^aj  P(D) 

■  0  for  a  dual  propeller 

The  propeller  pitching  moment  due  to  duct  interference  is  given 
by  Equation  ( 56)  as 


(^”*0)  (D— P) 

(^"a)  D(P) 


»,  0.451,  0.015 


(65) 


with  the  propeller  located  at  the  leading  edge,  midchord,  and 

trailing  edge  of  the  duct,  respectively.  Thus,  the  propeller 

interference  moment  is  smaller  than  the  duct  moment  unless  the 

propeller  is  located  near  the  leading  edge  of  the  duct. 

The  thrust  ratio  from  Equation  (11)  and  with  C„  “  2.17  is 

^P(D) 


-  0.20  (66) 

^P(D) 

so  that  the  duct  carries  only  16.7  percent  of  the  total  thrust. 

3.4.2  Pitching  (q)  derivatives 

The  pitching  derivatives  will  now  be  compared  at  •  0. 
The  duct  normal  force  derivative  is  from  Equation  (18) 
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The  side  force  derivative  of  a  single  propeller  is  from  Table  II 


/C„  N  -  -0.15  for  p  «  10° 

V  ^q>'P(D) 

-  -0.12  for  p  -  30° 

(68) 

Thus,  the  duct  normal  force  is  about  five  times  larger  than  the 
side  force  on  the  propeller. 

Now  let  us  compare  the  pitching  moment  on  the  duct  and  single 
propeller.  For  the  duct  (Eq.  (22))  we  have 

/C  \  -  +  0.0376  if-  -  0.029  (69) 

V  "*q;D(P)  '^o 

For  the  propeller  (Table  II) 

/'C™  \  -  -0.17  for  10  <  p  <  30°  (70) 

V  q> 

But  due  to  duct  interference  on  the  propeller  (Bq.  (60)),  we  get 

-12.0,  -10.2  (71) 

(S)d(p) 

for  leading-edge,  raidchord,  and  trailing-edge  propeller  locations, 
respectively.  Thus,  the  net  propeller  pitching  moment  reverses 
sign  with  propeller  position  as  mentioned  below  (Eq.  (60)).  It 
is  generally  large  compared  with  that  of  the  duct  and  it  is  dom¬ 
inated  by  duct-propeller  interference. 

3.4.3  Plunging  (d)  derivatives 

For  the  duct  normal  force  the  plunging  derivative  is  from 


Equation  (30) 
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0.13 


(72) 


This  force  acts  at  the  duct  midchord,  so  that  /c_  N  “0. 

The  plunging  derivative  of  the  propeller  normal  force  is  given 
by  Equation  (41)  .  Previously,  we  have  assumed 


.  4Bb  .  i 
°  37715  10 

Here,  we  shall  also  assume  a  three-bladed  propeller  (B 
that  D  ■  12.7b.  Then  we  have 


(73) 
3)  so 


\  -  0.058  sin®  P,  and  /C  N  -  0  (74) 

(,^d;P(D)  V  "d)  P(D) 


Thus,  the  apparent  mass  of  the  propeller  is  considerably  smaller 
than  that  of  the  duct,  and  the  interference  effect  of  the  duct  on 
the  propeller  is  probably  negligible  as  assumed  in  Section  3.3. 


3.4.4  Summary 

Specific  calculations  have  been  made  for  a  moderately  loaded, 
aircraft  type,  ducted  propeller  of  low  solidity.  The  propeller 
could  be  either  single  or  dual  for  a,  but  only  single  for  q- 

The  results  of  this  section  indicate  that  for  an  isolated 
ducted  propeller  of  this  type  (and  hence  for  certain  VTOL  con¬ 
figurations)  ,  all  of  the  forces  and  moments  on  the  duct  are  larger 
than  on  the  propeller  except  for  (1)  the  pitching  moment  due  to 
pitch,  and  (2)  the  thrust  force. 

The  motion  studied  consisted  of  small  rates  of  pitching  and 
plunging  about  a  straight-line  flight  path  aligned  with  the  duct 
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axis  because  of  the  analytical  restrictions  used  to  obtain  the  duct 
plunging  derivatives  and  all  of  the  propeller  stability  derivatives. 
The  midchord  duct  diameter  was  taken  as  the  axis  of  pitching  and  the 
center  of  duct  moments.  The  propeller  pitching  and  moments  were 
taken  about  its  own  diameter.  The  duct  and  propeller  axes  coincide 
for  c  «  D.  It  should  be  borne  in  mind  that  two  of  the  propeller 
reaction  components  were  perpendicular  to  their  duct  counterparts 
(a  yawing  moment  due  to  a  and  a  side  force  due  to  q) . 

The  duct  interference  on  the  propeller  markedly  increased  its 
stabilizing  moment  for  pitching  if  it  was  not  located  near  the  duct 
leading  edge.  This  can  be  significant  for  VTOL-type  applications. 

3.5  Calculative  Example  of  a  Torpedo-Like  Configuration 

In  this  section  the  foregoing  results  will  be  employed  to 
estimate  the  dynamic  stability  of  a  complete  torpedo-like  configu¬ 
ration.  The  purpose  is  to  determine  the  relative  importance  of  the 
ducted  propeller  derivatives  for  this  type  of  vehicle. 

We  shall  consider  small  rates  of  pitching  and  plunging  about 
an  equilibrium  condition  of  steady  axial  flight  with  the  thrust  of 
the  ducted  propeller  equal  and  opposite  to  the  hull  drag.  It  is 
found  that  this  condition  requires  a  very  lightly  loaded  propeller 
for  the  assumed  configuration  (Fig.  4)  where  the  propeller  diam¬ 
eter  is  taken  equal  to  the  maximum  hull  diameter.  With  this  assump¬ 
tion,  the  reference  areas  (A)  of  the  hull  and  ducted  propeller 
derivatives  are  the  same  so  that  they  are  directly  comparable.  The 
hull  shape  will  be  chosen  for  mathematical  simplicity  rather  than 
from  practical  considerations;  however,  for  the  present  purpose 
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it  is  probably  sufficiently  representative  of  actual  hull  shapes. 

For  purposes  of  the  present  stability  analysis,  the  hull  and  ducted 
propeller  will  be  assumed  to  be  isolated  from  one  another  (i.e., 
no  duct-hull  interference) ,  but  thrust  equilibrium  will  be  imposed. 
This  assumption  will  be  examined  qualitatively  at  the  end  of  the 
section. 

The  method  of  approach  will  be  to  estimate  first  the  isolated 
hull  stability  derivatives  about  its  centroid  by  use  of  slender 
body  theory  and  then  to  calculate  the  Isolated  ducted  propeller 
derivatives  about  the  hull  centroid  by  use  of  the  previous  analysis. 
The  static  stability  of  the  configuration  will  be  examined,  and 
sufficient  fins  will  be  added  to  the  duct  to  provide  a  10  percent 
margin  of  static  stability  (neglecting  aspect  ratio  and  interfer¬ 
ence  effects) •  Finally,  the  dynamic  stability  of  the  entire  con¬ 
figuration  will  be  considered  by  shifting  the  axis  of  pitch  and 
center  of  moments  from  the  centroid  to  the  assumed  location  of  the 
center  of  gravity  of  the  hull(0.4L  in  Fig.  4)  Interference  between 
the  hull  and  ducted  propeller  will  then  be  considered  briefly  in 
order  to  furnish  a  better  understanding  of  the  results  obtained. 

3.5.1  Isolated  hull  derivatives  about  hull  centroid 

The  stability  derivatives  of  an  isolated  slender  hull  (Fig.  5) 
will  be  estimated  by  using  slender  body  theory  and  by  assuming  that 
the  hull  boundary  layer  remains  attached  and  axially  symmetric. 

The  hull  probably  actually  oscillates  within  its  boundary  layer 
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to  some  extent.  Thus  at  small  angle  of  attack  the  latter  assump¬ 
tion  probably  leads  to  overestimation  of  the  effect  of  an  attached 
boundary  layer  and  the  hull  stability  derivatives.  However,  at 
higher  angles  of  attack,  which  are  not  considered  here,  the  boundary 
layer  actually  separates,  and  the  normal  force  on  the  hull  becomes 
much  greater  than  indicated  by  the  above  method.  The  total  drag  on 
the  hull  will  be  taken  equal  to  the  friction  drag,  and  its  form 
drag  will  be  neglected.  This  seems  justified  for  the  present  case 
where  the  boundary  layer  is  found  to  be  thin.  This  condition  is 
believed  to  be  generally  true  for  submerged  torpedoes  which  have 
little  wave  drag. 

Static  (g)  Derivatives 

The  normal  force  on  a  slender  hull  at  small  angle  of  attack  is 
given  by  Reference  10  (page  68) 

(75) 

a 


where  r  is  the  base  radius.  Also,  the  center  of  pressure  is  at 
L 

(Ref.  10) 


a  ■ 


L 


Vol 


(76) 


where  Vol  is  the  volume  of  the  body  and  L  is  its  length  (Fig.  5) . 
For  a  slender  body  with  a  pointed  base  at  small  angle  of  attack,  it 
is  assumed  that  the  effective  base  radius  r^  is  simply  6^^* »  the 
displacement  thickness  of  the  hull  boxindary  layer  at  the  base 
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(See  Ref.  10,  page  68) .  If  the  boundary  layer  is  axially  sym¬ 
metric,  then  the  friction  drag  of  the  body  is  given  by 
where  9.  is  the  momentum  thickness  of  the  boundary  layer  at  the 

ij 

base.  Thus,  if  we  neglect  form  drag  (assuming  the  boundary  layer 
be  thin) ,  then  the  total  drag  is  equal  to  the  friction  drag  and 
we  have 

X  -  2Tre^%^  (77) 


It  is  assumed  here  that  the  hull  is  a  body  of  revolution  (Fig.  5) 
for  which 


(78) 


and  the  hull  volume  is  from  Appendix  E  (Eq.(E.7)  with 

Vol  - 

15  m 


-  0) 
(79) 


Consequently,  from  the  above  equations,  the  axial  and  normal  forces 
of  the  hull  are  given  by 


o  in 


-  2 


m 


N, 


-  -o  ^ 

2  Op 


o  in 


m 


(80) 


(81) 


and,  for  no  boundary- layer  volume,  the  center  of  pressure 

8r  ^ 

3.  •  1  -  in 

^  156  *2 


J.9  a  u 


(82) 
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The  boundary  layer  profile  over  the  axisymmetric  hull  will  be  taken 
equal  to  that  given  by  the  l/7-power  law  (p.  432,  Ref.  11)  for  a 
turbulent  boundary  layer  over  a  flat  plate  of  length  L  in  uniform 
two-dimensional  flow.  Thus,  at  the  hull  base  (Fig.  5),  the 
boundary- layer  profile  is  given  by 

1/7 


t#) 


(83) 


for  which  (Ref.  11) 

6j^*  -  1.29  -  0.37  L  (Re^)  ®  (84) 

where  the  length  Reynolds  niomber  of  the  hull  is  RSj^  “  (pV^L)/pt. 
For  a  hull  moving  through  water  with  kinematic  viscosity  fx/p  ■ 
10~®  (ft^sec)  ,  we  can  write 

R  -  lO^V  L  (85) 

®L  ° 

Now,  if  L  “  10  ft  and  ■  40  (ft/sec)  ,  we  have  Rej^  ■  4x10^. 
Thus,  the  hull  boundary  layer  is  mostly  turbulent  (as  assumed), 
since  transition  to  turbulence  occurs  at  a  length  Reynolds  number 
of  approximately  2x10®  for  a  flat  plate® (or  at  5  percent  of  the 
hull  length) .  Substitution  of  the  assumed  values  (Fig.  4)  into 
Equation  (84)  gives 

6^*  -  0.111  ft  -  1.34  in.  (86) 

-  0.0865  ft  (87) 


^Kuethe,  A.  M.  and  Schetzer,  J.  D.:  Foundations  of  Aerodynamics 
J.  Wiley  and  Sons,  Inc.,  New  York,  p.  283. 
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and  the  drag  coefficient  and  normal  force  derivative  of  the  hull 
are,  from  Equations  (80)  and  (81) , 

-  0.0600  (88) 

f  C„  \  -  0.0994  (89) 

V  H 

Neglecting  the  boundary- layer  volume,  we  have  from  Equation  (82) 

I  -  -  9.78  (90) 


However,  with  the  displacement  thickness  of  the  assumed  hull 
boundary  layer  Included,  the  total  volume  of  hull  and  boundary 
layer  from  Equation  E.9  in  Appendix  E  is  given  by 

Vol*  -  0.709(7rr  ®L)  (91) 

in 


and  we  find  that 


a*/L  -  -13.6 


(92) 


We  shall  use  this  value  to  obtain  the  pitching  moment  derivative 
about  the  hull  centroid.  Thus  (see  Fig.  5)  we  have 

'  AiV(l  -  »*)  5  -  0.0994(10  .  13.6  X  20)  -  28.0 

V  Ot/  n  O  O 


(93) 


Pitching  (q)  derivatives 

The  stability  derivatives  for  the  hull  pitching  about  its 
centroid  are  given  by  the  following  formulas,  taken  from  p.  371  of 


Reference  10: 
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(94) 


+J5 


/ 

-k 


77  iT  T- 


(95) 


where  A  and  C  are  apparent  mass  coefficients  defined  in 

as  aa 

Reference  10,  and  the  different  reference  lengths  appearing  in  the 
above  formulas  require  conversion  to  the  present  system  of  nomen¬ 
clature.  Both  of  the  above  derivatives  are  zero  for  the  assumed 
hull  shape  if  we  neglect  the  boundary- layer  thickness.  With  the 
previous  assumptions  regarding  the  boundary  layer,  the  above  formulas 
become  (with  r,L  “  6*^^  and  “  -0.144  from  Appendix  E,  Eq.  (E.13) 


0.0990 


0.050  +  0.576  -  0.526 


(96) 

(97) 


The  above  q  derivatives  for  the  hull  are  due  entirely  to  the 
boundary  layer.  Now,  converting  to  the  present  system  of  nomen¬ 
clature,  we  obtain 


H 


0.0990  ^ 


0.495 


(98) 


H 


T  ^ 

0.526  — 
D® 


52.6 


(99) 
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Plunalna  (d)  PerlvatlveB 

The  plunging  (d)  stability  derivatives  for  the  normal  force 

Njj  and  pitching  moment  M'^  about  the  centroid  of  the  isolated 

hull  are,  from  pages  369  and  374  of  Reference  10, 

+J5 


where  (Fig.  5)  is  opposite  in  sign  from  Z  of  Reference  10. 
For  the  assumed  hull  shape  (neglecting  the  displacement  thickness 
of  the  boundary  layer) ,  the  above  formulas  become 


4  2lol.. 

TTT 

m 


2.13 


(102) 

(103) 


Now,  approximating  the  effect  of  the  boundary  layer  as  before  the 

effective  value  of  B  is 

es 


_  Vol* 

Trr 

in 


0.709 


(104) 


and  the  effective  value  of  C  from  Equation  E.13  in  Appendlk 

22 


E  is 
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C*  -  -0.144 

(105) 

22 

Therefore,  Equations 

(100)  and  (101)  for  the 

hull  with  boundary 

layer  become 

(106) 

/^C  .  ^  -  -0.576 

(107) 

\  aJn 

Thus,  for  pure  plunging  motion,  all  of  the  moment  and  a  small  part 
of  the  normal  force  on  the  hull  are  due  to  its  boundary  layer.  Now, 
converting  to  the  present  system  of  nomenclature,  we  obtain  for  d 
derivatives  of  the  hull  with  boundary  layer 

The  above  results  for  the  hull  are  summarized  in  Table  V. 

3.5.2  Isolated  ducted-propeller  derivatives  about  the  hull  centroid 
Here  we  shall  determine  the  derivatives  of  the  ducted  propeller 
mounted  at  the  base  of  the  hull  in  Figure  4  by  use  of  the  results 
of  Section  3.  In  order  to  do  this,  we  must  assume  that  the  ducted 
propeller  is  Isolated  from  the  hull  but  is  in  thrust  equilibrium 
with  it.  That  is,  we  shall  neglect  all  interference  between  the 
hull  and  ducted  propeller,  including  the  effect  of  the  hull  boundary 
layer  flowing  into  the  ducted  propeller.  Neglecting  this  latter 
effect  is  probably  justified,  since  the  displacement  thickness  of 
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the  boundary  layer  at  the  hull  base  was  found  to  be  1.34  inches, 
which  Is  reasonably  small  compared  with  the  assumed  duct  radius 
of  6  inches. 

We  shall  assume  that  the  propeller  is  located  toward  the  rear 
of  the  duct,  so  that  the  duct  interference  on  the  propeller  from 
q  is  not  large  and  destabilizing.  (See  Eg.  (71)) 

Thrust  Equilibrium 

The  propeller  disk  loading  is  from  Equation  (170)  with  ”  0 


also,  from  Equation  (11) ,  we  find  that  the  duct  thrust  (for  “ 
0  and  c  -  0.2D)  is 


'^D(P) 


0.724 

1  +  2{V^/y) 


^P(D) 


(111) 


Thus,  the  total  thrust  of  the  ducted  propeller  is  given  by 


(112) 


For  thrust  equilibrium  of  the  entire  vehicle,  we  have  '^D(p)  "*■ 
Tp(D)  “  X.  Using  Equation  (88)  for  the  drag  coefficient  of  the 
hull,  we  obtain 


Cy  -  -  1.724  ^  +  2  ^  -  0.060 

o 


(113) 
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so  that 


^  -  0.050  (114) 

''o 

Thus,  for  the  assumed  configuration,  the  slipstream  velocity  Vj 
is  only  5  percent  greater  than  the  flight  speed  and  the 

propeller  disk  loading  is  low.  We  shall  now  determine  the  stabil¬ 
ity  derivatives  of  the  isolated  ducted  propeller  about  the  hull 
centroid  for  this  case,  namely  ■  0,  c  ■  0.2D,  and  7  ■  0.05  V^. 


Static  (g)  Derivatives 

The  normal  force  and  pitching  moment  derivatives  of  the  duct 
are  found  from  (Eqs.  (9)  and  (13))  to  be 


^C„  -  1.91  +  0.845  ^  -  1.95  (115) 

V*’a^(P)  ^o 

/^c  ^  -  1.45  ^  -  0.0725  (116) 

(^"'a;D(P)  ^o 


The  above  pitching  moment  is  about  the  midchord  duct  diameter,  so 
that  the  moment  about  the  centroid  of  the  hull  due  to  the  duct 
(Fig.  4)  is  given  by 


LN, 


M' 


D(P) 


a 


»t(P) 


b(P) 


Si 


a 


Thus  we  find  that 


-  0.07  -  19.5  -  -19.4 


D(P) 


(118) 
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where  the  first  small  term  is  the  duct  moment  about  its  own  center 
and  the  second  is  the  moment  about  the  hull  centroid  due  to  the 
duct  normal  force. 

The  duct  a  derivatives  are  listed  in  Table  III  together 

with  the  propeller  derivatives  estimated  from  the  results  of  Section 

3.4.1  for  either  a  dual  or  single  aircraft  type  propeller  located 

near  the  rear  of  the  duct  with  v  •  0.05  V  .  It  can  be  seen  that 

^  o 

the  normal  force  on  the  propeller  is  much  smaller  than  on  the  duct 
and  that  the  pitching  moment  about  the  hull  centroid  due  to  the 
propeller  is  much  smaller  than  that  due  to  the  duct.  It  should  be 
noted,  however,  that  there  is  a  smaller  yawing  moment  on  the 
propeller  which  is  perpendicular  to  the  duct  moment. 

Pitching  (q)  Derivatives 

The  pitching  derivatives  of  the  duct  for  rotation  about  the 
hull  centroid  are  (from  Eqs.  (9)  and  (18))  : 

^  “  0.675  ^  +  0.126  +  (^.91  +  0.845  )-  9.91 

(119) 


and  from  Equations  (13)  and  (22) 


0.0376  + 


0.0725 


0.364 


(120) 


The  last  term  in  each  expression  is  due  to  the  induced  angle  of 
attack  which  is  introduced  at  the  ducted  propeller  by  rotation 
about  the  hull  centroid. 
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-  ak_ 

2^0 


(121) 


The  pitching  moment  about  the  centroid  of  the  hull  (denoted 
by  a  prime)  due  to  the  duct  is 


“'D(P)q  “  ”D(P)q  “  2  *^D(P)q 

so  that 

C„.  \  -  0.364  -  99.1  -  -98.7 

q>(P) 

Note  that  the  duct  moment  about  its  own  center  ^(p)  negli- 

gible,  just  as  was  *^(p)  • 

ct 

The  duct  q  derivatives  are  listed  in  Table  IV  for  comparison 
with  the  propeller  derivatives  as  estimated  from  the  results  of 
Section  3.4.2,  for  an  aircraft  type  propeller  located  near  the  rear 
of  the  duct,  with  v  0.05  V  . 

It  can  be  seen  that  the  forces  and  moments  acting  at  the  hull 
centroid  due  to  the  propeller  are  small  compared  with  those  due  to 
the  duct.  It  should  be  noted,  howfever,  that  there  is  a  somewhat 
smaller  side  force  and  yawing  moment  on  the  propeller  which  are 
perpendicular  to  the  pitching  motion. 

The  above  results  show  that  the  q  damping  of  the  ducted  pro¬ 
peller  is  due  almost  entirely  to  the  duct  normal  force  arising  from 
rotation  about  the  hull  centroid.  Equations  (119),  (120),  and  (122) 
show  that  the  disk  loading  has  a  small  effect  on  pitch  damping  about 
the  hull  centroid  for  the  present  example,  even  for  high  disk  loading. 
Plunging  (d)  Derivatives 

For  the  duct,  the  net  hydrodynamic  reaction  due  to  d  is  a 
normal  force  acting  through  the  midchord  plane,  and  from  Equations 


(122) 

(123) 


50- 


(30)  and  (31)  we  obtain 

(%)d(P) 


and 


0.126 


(124) 


(125) 


The  pitching  moment  about  the  hull  centroid  due  to  the  duct  is 
then 


M' 


D(P) 


a 


-  ^  N 
2  ^D(P), 


(126) 


and  the  corresponding  stability  derivative  is 

(^m  -  n  (0-126)  -  -1.26  (127) 

V  "d  ^(P)  ° 

From  the  results  of  Section  3.4.3,  (Eq.  (72)  and  (74))  the 
estimated  d  reaction  from  the  apparent  mass  of  the  propeller  is 
negligibly  small  compared  with  that  of  the  duct. 

3.5.3  Comparison  of  isolated  hull  derivatives  and  isolated  ducted 
propeller  derivatives  about  hull  centroid 

The  stability  derivatives  which  have  been  obtained  for  the 
hull  and  ducted  propeller  about  the  hull  centroid  are  summarized 
in  Table  V  for  comparison.  It  should  be  recalled  that  these 
results  were  obtained  for  the  configuration  in  Figure  4,  and  for 
small  oscillations  from  axial  flight  at  low  frequency  as  assumed 
in  the  d  analysis .  The  hull  and  ducted  propeller  are  con¬ 
sidered  isolated  from  one  another  but  are  in  thrust  equilibrium. 

An  aircraft-type  propeller  is  located  near  the  duct  exit 


1 
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plane.  The  hull  boundary  layer  is  assumed  to  be  attached 

£ 

and  axially  symmetric.  \ 

The  following  conclusions  are  drawn  from  Table  V  and  the  * 

foregoing  analysis. 

(1)  The  q  damping  moment  about  the  hull  centroid  ') 

from  the  propeller  is  small  compared  with  that  due  to  the  duct, 
even  though  the  propeller  moment  about  its  own  center  was  compar- 
itively  large. 

(2)  The  upwash  over  the  duct  arising  from  the  q  motion 
causes  a  stabilizing  moment  about  the  hull  centroid  which  is  greater 
than  the  destabilizing  moment  from  the  hull  boundary  layer. 

(3)  The  d  damping  )  is  due  mainly  to  the  apparent 

™d 

mass  of  the  hull  boundary  layer. 

(4)  Both  C  '  and  C  '  for  the  hull  are  due  entirely  to 

'"a 

a  q 

its  boundary  layer  for  a  hull  with  a  pointed  base.  Since  these 
damping  components  are  comparatively  large,  refinement  of  the 
present  analysis  is  needed.  Both  of  these  components  are  probably 
over-estimated  by  the  present  analysis  for  small  q  and  d  motions 
of  the  hull,  since  the  analysis  assumes  that  the  hull  boundary 
layer  remains  axially  symmetric. 

(5)  For  a  torpedo-like  configuration,  all  of  the  propeller 

forces  and  moments  (except  thrust)  are  small  compared  with  those  ] 

of  the  duct.  However,  some  of  the  propeller  reactions  are  perpen¬ 
dicular  to  the  motion  which  causes  them.  This  is  not  true  of  the 
duct .  I 

I 

I 

ft 

I 
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3.5.4  Static  stability  of  torpedo-like  configuration 

Before  we  can  make  a  dynamic  stability  analysis  of  a  complete 
torpedo-like  configuration,  we  must  insure  that  the  vehicle  is 
trimmed  and  is  statically  stable.  For  this  purpose  we  shall  now 
add  fins  to  the  configuration  (Fig.  4)  and  determine  the  static 
derivatives  of  the  entire  configuration  consisting  of  a  hull, 
ducted  propeller,  and  fin.  We  shall  not  consider  the  comparatively 
small  propeller  forces  and  moments  (Table  m)  .  The  only  interfer¬ 
ence  effect  which  is  to  be  included  is  propeller  interference  on 
the  duct,  so  that  the  total  normal  force  and  pitching  moment  on 
the  entire  configuration  are  given  by 

and 

"huff"-  V*  «D(P)''-"  V 

The  double  prime  superscript  on  the  moment  indicates  that  it  is 
about  the  center  of  gravity  of  the  hull.*  We  shall  use  the  moment 
transfer  equations  of  page  400,  Reference  10  to  shift  the  centers 
of  moments  and  pitch  from  the  hull  centroid  (previously  considered) 
to  the  hull  center  of  gravity.  The  normal  force  derivatives  due  to 
a  and  d  are  not  affected  by  this  transfer  and  can  be  taken 
directly  from  Table  V.  All  the  other  stability  derivatives,  how¬ 
ever,  are  affected. 

*It  is  assumed  that  the  center  of  gravity  of  the  hull  is  the  same 
as  that  of  the  complete  vehicle. 
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The  center  of  gravity  of  the  entire  configuration  Is  assumed 
to  be  at  0.4L  (Fig.  4).  Using  Table  V,  we  find  that  the  center  of 


pressure  of  the  hull  and  duct  Is  located  a  distance  x 
of  the  hull  centroid  as  given  by 


c.p. 


ahead 


MV 


+  M' 


H. 


a 


D(P) 


a 


CL  a 


(28.0  -  19.4)D/2L 
0.10  +  1.95 


0.21 


(130) 


Thus  the  center  of  pressure  Is  located  0.29L  from  the  nose  In 
Figure  4  and  Is  ahead  of  the  assumed  center  of  gravity.  Since  this 
results  In  static  Instability,  we  must  add  fins  to  make  the  config¬ 
uration  statically  stable.  We  shall  add  fins  near  the  base  of  the 
hull  (Fig.  4)  which  are  sufficiently  large  to  place  the  center  of 
pressure  at  the  hull  centroid  (0.5L)  and  thereby  obtain  a  "10- 
percent  static  margin."  We  shall  use  a  fin  with  chord  c  and  ex¬ 
posed  span  S ,  such  that  the  net  moment  about  the  hull  centroid  Is 
zero.  That  Is, 

M'f  +  -  0  (131) 

Ct  0^  CL 

Thus,  using  Table  V,  we  find  that 

«'f  ■  -  2 

a 
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Now  we  shall  neglect  interference  effects  on  the  fin,  assume  that 
its  lift  curve  slope  is  lir  and  make  no  correction  for  its  aspect 
ratio.  Then  we  can  write  the  normal  force  on  the  fin  as 

2M’f  2 

Np - JpS  -  e.6  £  Aq^a„  .  0.86  -  2ircSq  q 

(133) 

Thus,  if  the  fin  chord  is  taken  equal  to  the  duct  chord  (c  ■  0.2D) , 
we  find  that  the  total  required  exposed  fin  span  is 

S  -  0.537  D  (134) 

With  this  fin  added  for  static  stability,  the  static  a  de¬ 
rivatives  about  the  hull  center  of  gravity  (0.4L)  as  obtained  from 
Table  V  using  the  transfer  Equations  (10-119)  of  Reference  10  are  given 
in  Table  VI . 

3.5.5  Dynamic  stability  of  complete  configuration  with  no  inter¬ 
ference  between  hull  and  ducted  propeller 

The  dynamic  stability  derivatives  for  the  hull  and  duct  about 
the  vehicle  center  of  gravity  can  be  found  by  transferring  their  de¬ 
rivatives  about  the  centroid  (Table  V)  using  Equations  (10-119)  of 
Reference  10.  The  results  are  given  in  Table  VI. 

For  the  fin,  we  shall  neglect  the  moment  about  its  own  axis, 
any  interference  effects,  and  aspect  ratio  corrections.  The  upwash 
over  this  idealized  fin  due  to  pitching  about  the  hull  center  of 
gravity  produces  the  following  derivatives : 
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(\)  r  -  (X) 


(136) 


The  normal  force  and  moment  on  the  fin  due  to  rotation  about  its 
own  axis  are  comparatively  small  and  have  been  neglected  here. 

Due  to  d,  there  is  a  normal  force  on  the  fin  (Section  3.6) 
given  by  Equation  (166) 


N., 


a 


N, 


b(p) . 
a 


2  S 
TT  D 


^  -  0.298 


(137) 


thus  we  find  that 


('c„  \  -  0.298(0.126) 

V 

{'c  "'s  -  0.298(-1.51) 

V  ”'a)F 


0.038 


-0.45 


(138) 

(139) 


The  resulting  stability  derivatives  about  the  hull  center  of  gravi¬ 
ty  with  ducted  propeller  and  fins  are  tabulated  in  Table  VI  for 
the  entire  torpedo- like  configuration. 

The  following  comments  should  be  made  concerning  the  results 
summarized  in  Table  VI : 

(1)  For  the  hull,  both  the  stabilizing  d  moment  and  the 
destabilizing  q  moment  are  due  entirely  to  its  boundary  layer, 
which  has  been  assumed  to  remain  axially  symmetric.  For  the  small 
oscillations  from  steady  axial  flight  considered  here,  these 
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moments  are  probably  overestimated  by  this  assumption,  since  the 
boundary  layer  oscillation  will  actually  differ  from  that  of  the 
hull. 

(2)  For  the  duct,  the  stabilizing  moments  about  the  vehicle 
center  of  gravity  (with  a  and  q)  are  due  primarily  to  the  normal 
force  which  Is  generated  by  the  upwash  over  the  duct.  The  moments 
about  the  midchord  diameter  of  the  duct  are  negligible  In  compari¬ 
son  with  these  moments.  In  this  analysis,  however,  we  have  allowed 
no  Interference  between  the  hull  and  ducted  propeller,  and  we  have 
neglected  the  effect  of  the  hull  boundary  layer  on  the  propeller 
since  It  Is  reasonably  small  compared  with  the  duct  radius  (6*^^  ■ 
1.34  In. ,  R  ■  6  In. ) . 

(3)  We  have  assumed  an  aircraft  type  propeller  (of  low  solid¬ 
ity  with  slender  blades)  which  Is  located  near  the  exit  plane  of 
the  duct.  For  such  a  configuration  the  only  appreciable  effects  of 
the  propeller  are  Its  thrust  force  and  Its  Interference  upon  the 
duct.  However,  small  forces  and  moments  are  generated  on  the  pro¬ 
peller  which  are  perpendicular  to  the  motion  which  causes  them. 

The  following  conclusions  are  drawn  from  Table  VI  for  the 
overall  configuration: 

(1)  The  static  stability  "  Is  Influenced  appreciably  by 

OL 

each  component;  the  hull,  the  ducted  propeller,  and  the  fin. 

(2)  The  q  damping  "  from  the  ducted  propeller  and  fin 
are  both  strongly  stabilizing  and  outweigh  the  destabilizing  effect 
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of  the  hull  boundary  layer,  even  though  this  has  been  overestimated 
here. 

(3)  The  d  damping  "  is  due  almost  entirely  to  the  hull 

a 

boundary  layer  and  has  been  overestimated  here. 

(4)  The  deunping  in  pitch  ")  is  stabilizing  and  is 

q  d 

produced  mostly  by  the  ducted  propeller.  The  fins  provide  about  30 

percent  of  the  damping  and  the  contribution  of  the  hull  is  negligible 

with  the  boundary  layer  producing  opposite  effects  on  "  and  " . 

Ot  q 

Now  we  shall  estimate  the  undamped  natural  frequency  and 

the  damping  ratio  (or  percent  critical  damping)  for  oscilla¬ 

tion  of  the  entire  configuration,  the  undamped  natural  frequency  is 
given  by  (Eqs.  (8-85)  ,  Ref.  10) 


n 


■  -  C. 


m 


a 


q^A(D/2) 

•"cW" 


(140) 


and  the  damping  ratio  for  impulsive  pitch  control  in  steady  flight 
is  (Eq.  (8-86),  Ref.  10) 


:  -  Tc  "  +  c 

Na  "d  ; 


2K 


V 


2m^V  a)_ 
con 


(141) 


where  m  is  the  mass  of  the  entire  configuration  and  K  "  is  the 
c  y 

radius  of  gyration  of  entire  configuration  about  the  hull  center  of 
gravity. 

We  shall  take  m^  and  K„"  to  be  those  for  the  hull  alone 

c  y 

from  Appendix  E  as  follows: 


58 


m. 


Vol 


irr^^tp 


(142) 


Ky"®  -  0.0457 


(143) 


Substitution  of  these  values  and  those  In  Table  VI  and  Figure  4 


Into  Equation  (140)  yields  for  the  natural  frequency 

2 


0)  ^  -  C  "  % «  /  A 

n  m  32(0 

Cw 


iS _ rZaV 

0.456)  \  L  y 


-  -  16.4  C. 


m 


d 


-  95.0 


(0  “9.75  rad/sec 

n 


The  undamped  period  of  oscillation  Is  then 


T  -  —  -  0.66  sec 


and  the  damping  ratio  Is 

C  "  +  C 
in_ 


X  4.57  /\32  Uo^J 

^2.  91  +  — 1^  X  4 — ^ 

8  X  4.57y\l5  X  9.75y 


(144) 


(145) 


(146) 


'•  1.67 


(147) 


so  the  torpedo  Is  heavily  damped. 

This  ratio  Is  Influenced  appreciably  by  each  component  of  the  con¬ 
figuration  (the  hull,  the  ducted  propeller  and  the  fin)  as  seen 


from  Table  VI 
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3.5.6  Interference  between  hull  and  ducted  propeller 

In  the  above  calculative  example,  the  stability  derivatives 
of  an  Isolated  slender  hull  were  used  along  with  those  for  an  iso¬ 
lated  ducted  propeller  with  which  it  was  in  thrust  equilibrium, 

X  ■  '^D(p)  '^p(D)  *  other  words,  the  mutual  interference  between 

the  hull  and  ducted  propeller  was  neglected.  To  obtain  some  in¬ 
sight  into  the  order  of  magnitude  of  this  interference  effect,  we 
shall  estimate  here  the  thrust  interference  for  steady  axial  flight. 
This  will  be  done  by: 

(1)  Assuming  the  hull  drag  to  be  independent  of  hull  shape 
and  given  by  the  flat  plate  boundary  layer  of  Section  3.5.1  (Eq. 
(84))  . 

(2)  Approximating  the  duct  bound  vorticity  distribution 
(Eq.  (A. 50))  by  the  expression 

-  Y  +  2K^  tan  j  (148) 

where  y  is  the  actuator  disk  vorticity  and  is  a  constant  to 

be  determined . 

(3)  Representing  the  hull  shape  by  the  dividing  stream 
surface  about  a  concentrated  point  source  and  sink  (Sketch  C) . 

(4)  Choosing  the  value  of  such  that  there  is  no  lateral 

velocity  at  the  midchord  point  of  the  duct. 

(5)  Solving  for  the  duct  thrust  from  the  leading-edge  singu¬ 
larity  represented  by  K^. 
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In  this  analysis,  the  hull  shape  is  allowed  to  change  accord¬ 
ing  to  the  Inflow  of  the  ducted  propeller  (Sketch  D) ,  and  the  effect 
of  the  hull  boundary  layer  entering  the  ducted  propeller  will  be 
neglected.  Thus,  the  results  will  give  only  an  order-of-magnltude 
estimate  of  the  mutual  Interference. 


Sketch  C.  -  Hull  and  ducted-propeller  nomenclature. 

The  above  analytical  approximations  will  first  be  tested  by 
solving  for  the  value  of  duct  thrust  (for  the  isolated 

d 

ducted  propeller)  using  the  above  approximate  bound  vorticlty 
distribution  and  comparing  this  value  with  Tp(p)  obtained  from 
the  more  exact  theory  (Eq.  (11)).  For  the  approximate  flow  model, 
the  duct  thrust  is  found  to  be  expressible  as 
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and  from  Equation  (11)  we  have 

-  f  p  ^  ^ 

For  cases  in  which  c/D  <  1/8,  we  find  that 

■^DCP) 

0.89  <  = - a  <  1  (151) 

''D(P)e 

Thus,  the  approximate  value  of  duct  thrust  is  reasonably  good  for 
very  short  ducts.  Also,  the  total  circulation  about  the  duct  chord 
is  found  to  be  within  8  percent  of  the  more  exact  value  for  c/D  < 
1/8. 

In  the  absence  of  the  ducted  propeller,  the  hull  shape  (Ref. 
12,  p.  461)  is  axially  symmetric  and  is  symmetric  fore  and  aft. 
Also,  for  a  very  slender  hull,  we  have 

L  =  2S 

and 

Q  =  irrjiv/l)  (152) 

m  o 

In  the  presence  of  the  ducted  propeller,  the  hull  is  no  longer 
symmetric  fore  and  aft,  but  the  forward  half  of  the  hull  is  not 
significantly  distorted  (see  Sketch  D) .  The  duct  thrust  in  the 
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presence  of  the  hull  and  propeller  (as  calculated  by  the  method 
described  above)  Is  given  approximately  by 


'^D(HP).  "  "^dCP) 


ln(8D/c)  -  2 


(153) 


where  It  has  been  assumed  that  s  «  D  «  L  and  r^  «  L.  The 
duct  thrust  In  the  presence  of  the  propeller  Is  therefore  In¬ 
creased,  due  to  hull  Interference,  by  an  amount 


''’h-^(P)  “  '^D(HP)  "  ’’d(P) 

o  a 


(154) 


On  the  other  hand,  the  pressure  drag  of  the  hull  due  to  Inter¬ 
ference  from  the  ducted  propeller  Is  found  from  the  Blaslus  formula 
(p.  168,  Ref.  12)  to  be 


^(P) 


°  2d®v^  \r 


Thus  we  find  that  the  ratio  of  Increased  thrust  to  Increased  drag 
Is 

.  2  ln(8D/e)  -  4  .  r„VD^  _ 

^(P)  -►  H  ln(8D/c)  +  f 

(156) 

and  there  Is  a  net  loss  of  thrust  for  the  entire  configuration  due 
to  duct-hull  Interference.  One  would  expect  this  effect  to  be  even 
greater  without  hull  distortion.  One  can  Imagine  that  the  inter¬ 
ference  forces  are  generated  by  moving  the  ducted  propeller  from  far 
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downstream  to  a  point  near  the  base  of  the  hull  as  shown  in  Sketch 
D.  Now,  since  the  propeller  thrust  (or 

7/V^)  must  be  Increased  to  maintain  thrust  equilibrium.  If  there 
were  no  propeller  in  the  duct,  that  is,  if  y  ■  0,  then  the  hull 
and  duct  interference  forces  would  of  course  be  equal  and  opposite 
but  thrust  equillbritun  would  not  be  possible  without  the  propeller 


y 


*  ’^D(P)  ^P(D)  ^H-^D(P) 


Sketch  D  -  Hull  and  ducted  propeller  interference  forces. 


Using  the  values  estimated  in  Section  3.5.1  for  a  different  (para¬ 
bolic)  isolated  hull  shape  "  0*050,  c  “  0.2D,  D  ■  2r^t  and 

L  "  IOd'^  ,  we  find  from  Equation  (153)  that 


^D(HP)a  “  '^D(P)a 


1  + 


0.050(4)1.69 


106  T, 


D(P)a 


(157) 
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Thus,  the  approximate  duct  thrust  Increases  tremendously  because 
of  interference  with  a  blunt  base  hull.®  However,  this  is  due  to 
the  fact  that  the  ducted  propeller  is  very  lightly  loaded  when  it 
is  isolated  from  the  hull  (far  downstream  in  Sketch  D) .  When  they 
are  brought  together,  the  hull  Induces  a  large  radial  velocity 
component  across  the  duct  as  compared  with  the  lightly  loaded  pro¬ 
peller  (Sketch  D) .  Thus,  the  duct  thrust  due  to  hull  Interference 
is  much  larger  than  that  required  to  overcome  the  hull  drag  for 
the  present  configuration.  If  the  propeller  had  been  more  highly 
loaded  so  that  7  “  exaunple,  then  the  duct  thrust  would 

have  increased  by  a  factor  of  only  2.15  rather  than  106. 

The  above  analysis  indicates  that  large  Interference  effects 
may  exist  between  the  ducted  propeller  and  hull,  and  that  more 
precise  analysis  should  therefore  be  made  to  obtain  valid  quantita¬ 
tive  results.  Since  both  the  hull  and  the  propeller  induce  a  radial 
in-flow  which  results  in  a  leading-edge  suction  on  the  duct,  the 
effect  of  hull  interference  on  the  duct  stability  derivatives  should 
be  similar  to  the  effect  of  an  increase  in  propeller  loading  or 
y/V^.  One  aspect  of  this  observation  is  that  duct  shapes  which 
would  actually  be  used  in  practice  to  propel  hulls  with  blunt  bases 
would  likely  have  the  duct  diameter  decreasing  in  the  streeunwise 
direction  in  order  to  prevent  flow  separation  (in  the  absence  of 
boundary- layer  control).  Thus,  in  view  of  the  large  effect  of  pro¬ 
peller  loading  on  hull-duct  Interference  forces,  a  reasonable  es¬ 
timate  of  these  interference  forces  would  require  that  the  actual 

^ote  that  Equation  (151)  is  considerably  less  accurate  for 
c/D  "0.2  than  for  c/D  ^  1/8. 
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hull  and  duct  shape  and  the  appropriate  propeller  loading  be  approx¬ 
imated  more  closely  than  Is  possible  within  the  scope  of  the  present 
project*  This  Is  particularly  true  for  the  nonaxlal  and  unsteady 
flight  conditions. 

3.6  Comparison  of  Duct  Normal  Force  with  that  of  an  Equivalent  Fin 

It  would  be  of  Interest  to  know  whether  a  ducted  propeller  Is 
more  or  less  effective  In  damping  than  a  fin  having  the  same  pro¬ 
jected  area.  For  this  purpose,  we  shall  compare  the  normal  forces 
on  a  duct  In  the  presence  of  a  propeller  with  those  on  such  a  fin 
for  small  a,  q,  and  d  motions.  The  propeller  forces  themselves 
are  known  to  be  comparatively  small  (Section  3.4)  and  will  there¬ 
fore  not  be  considered. 

The  duct  normal  force  ^^(p)  predicted  In  Section  3.1  (Eq. 

(9))  will  be  compared  with  Np  for  a  fin  with  chord  c  and  total 
exposed  span  D  where  c  "  0.2D  for  both  the  duct  and  the  fin. 

Thus  the  fin  planfozm  Is  the  projected  side  area  of  the  duct.  The 
fin  will  be  considered  as  a  section  of  an  Infinitely  long  plate; 
that  Is,  no  aspect  ratio  corrections  will  be  made.  For  c“  0.2D 
(aspect  ratio  of  5)  these  corrections  are  small. 

First,  let  us  consider  the  effect  of  a  small  angle  of  attack 
Oq.  For  the  duct,  we  have  from  Equation  (9) 

^C„'\  -  1.91  +  0.845  ^  (158) 

\  “a/D(P)  ^o 


and  for  the  fin,  using  a  lift  curve  slope  of  2ir,  the  normal  force  Is 
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-  Np^  -  cDq^ZTO^ 

Hence 


(159) 


D(P} 

s - -  'if 1'91  +  0.845  5?-  1-  1.19  +  0.528  ^  (160) 

8  V^J 

Thus,  even  without  a  propeller  (y  ■  0),  we  see  that  is  19 

percent  greater  than  Np.  This  result  is  nearly  independent  of 

c/D  when  c  «  D  as  can  be  seen  from  Equations  (9)  and  (159) .  in 

Reference  7  an  identical  expression  for  is  derived  for  the 

duct  alone  (a  “  0)  and  it  is  shown  therein  that  is  exactly 

twice  the  lift  of  a  fin  of  elliptical  planform  with  a  total  span 
equal  to  D  and  a  maximum  chord  equal  to  c.  Increasing  the 
propeller  loading  (or  further  increases  the  effectiveness 

of  the  ducted  propeller  over  the  fin. 

Now  consider  due  to  pitching  q  about  the  midchord 

duct  diameter.  From  Equation  (18)  we  find  that 


^  "  n  ^  +  TT  —  -  0.675  ^  +  0.126 

(161) 


for  c  ■  0.2D  and  a  “  0. 

o 

For  the  flat  fin  pitching  about  its  midspan  with  ■  0, 

N_  corresponds  to  that  for  an  equivalent  cambered  fin  (Sketch  A. 7) 
so  that  the  normal  force  is 
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where 


Np  ■  cDq  27ra' 

q 


a' 


2h  .  ccL 

c  4V^ 
o 


(162) 


Thus,  we  find  for  pitching 


N, 


D(P), 

"nI 


-  1.57  +  8.45 


(163) 


The  increased  effectiveness  of  the  ducted  propeller  over  the  fin  In 
pitch  therefore  increases  with  and  decreases  with  c/D. 

Finally,  let  us  consider  the  normal  force  due  to  plunging, 
d  ■  -v/V^  at  ■  0.  For  the  duct,  we  have  from  Equation  (28) 

cx 

For  the  fin,  we  have  from  Equation  (27) 


N. 


(165) 


Thus  for  any  value  of  c/D 


IT 

“f.  2 

a 


1.57 


It  is  therefore  concluded  that  the  duct  normal  force  due  to  a,  q,and 
d  is  larger  than  the  fin  normal  force  for  all  values  of  c/D  and 
y/V^.  Furthermore,  the  effectiveness  with  q  and  a  Increase  as 
the  propeller  loading  is  increased. 
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4.  COMPARISON  OF  EXPERIMENTAL  DATA  AND  THEORY 
4.1  Introduction 

In  this  section,  the  predicted  forces  and  moments  on  the  duct 
in  the  presence  of  the  propeller  will  be  compared  with  every  known 
source  of  applicable  ducted  propeller  data.  The  sources  of  data 
will  be  selected  from  the  references  listed,  which  are  summarized 
in  Reference  2.  We  shall  use  only  the  data  for  which  the  forces 
and  moments  on  the  propeller  and  on  the  duct  were  measured  inde¬ 
pendently  and  where  the  duct  was  reasonably  short,  straight,  and 
thin  as  assumed  in  the  present  theory.  Wherever  a  choice  is  possible, 
only  the  shortest  duct  tested  will  be  considered  here.  Most  of 
these  data  are  for  steady  hovering  flight  or  for  axial  flow,  in  which 
case  only  thrust  forces  exist  on  the  duct  and  propeller  and  for  which 
the  theory  predicts  only  the  duct  thrust  coefficient  in  terms  of 
the  propeller  thrust  coefficient. 

Two  sources  of  data  for  steady  flow  at  angle  of  attack  (Refs.  13 
and  14)  will  be  used  for  comparison  with  the  predicted  normal  force 
coefficient  of  the  duct,  but  only  one  of  these  (Ref.  13)  contains 
measurements  of  the  duct  pitching  moment  coefficient  which  will 
also  be  compared.  The  most  recent  ducted  propeller  data  at  angle 
of  attack  (Ref.  15)  will  not  be  used  here,  because  tne  duct  chord 
to  diameter  ratio  (c/D  ■  0.61)  was  not  small  as  assumed  in  the 
present  theory. 

The  only  known  source  of  dynamic  data  (Ref.  16)  will  not  be 
used  here  for  comparison,  because  the  ducted  propeller  was  mounted 


near  a  large  hull  and  only  the  forces  and  moments  on  the  entire 
configuration  were  measured.  Since  the  hull  interference  and 
boundary-layer  effects  upon  the  ducted  propeller  were  probably  not 
small  enough  to  neglect,  it  does  not  appear  worthwhile  to  attempt 
the  extraction  of  isolated  ducted  propeller  data  from  these  measure 
ments.  Thus  no  dynamic  data  will  be  compared  with  the  theory. 

Even  for  the  data  which  are  used,  the  validity  of  the  compar¬ 
ison  with  the  present  theory  is  questionable  because  of  one  or 
more  of  the  following  factors : 

(a)  Flow  separation  from  the  duct  inlet 

(b)  Duct  chord/diameter  ratio  too  large 

(c)  Nonuniform  propeller  loading 

(d)  Slipstream  swirl  generated  by  a  single  propeller 

(e)  Too  much  duct  taper 

(f)  Large  centerbody 

The  present  theory  predicts  the  forces  and  moments  on  the 

duct  in  the  presence  of  the  propeller  by  representing  the  propeller 

as  a  uniformly  loaded  actuator  disk.  The  duct  force  and  moment 

coefficients,  C-  ,  C„  ,  and  C  ,  are  predicted  as 
’^D(P)  *^D(P)  “d(P) 

functions  of  the  independent  variables  a^,  c/D,  and  y/V^.  There 
fore,  in  order  to  predict  the  duct  data  for  a  given  test,  we  must 
first  express  yA^  in  terms  of  the  measured  data.  For  the  ideal, 
uniformly  loaded  actuator  disk  and  the  corresponding  ideal  slip¬ 
stream  with  constant  velocity  and  no  swirl,  this  is  done  quite 

simply  by  use  of  Bernoulli's  equation  far  ahead  and  behind  the 
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propeller  disk.  However,  this  model  is  evidently  too  idealized 
for  a  satisfactory  comparison  of  theory  and  experiment.  Therefore, 
in  order  to  improve  the  comparison  of  the  present  theory  with  data 
for  real  propellers,  we  shall  now  consider  real  slipstream  effects. 
These  can  be  categorized  as: 

(a)  Slipstream  swirl  generated  by  the  propeller 

(b)  Nonuniformity  of  the  jet  velocity  component  Vj ,  which  is 

aligned  with  the  axis  of  the  slipstream  (Fig.  6) 

(c)  Frictional  losses  which  appear  as  propeller  blade  wakes 

Dual,  counter-rotating  propellers  or  propeller-stator  combina¬ 
tions  are  seldom  used  in  the  data  described  below.  When  such  com¬ 
binations  are  used,  the  slipstream  swirl  should  be  comparatively 
small  conqpared  with  that  for  a  single  propeller. 

For  the  purpose  of  improving  the  comparison  of  the  present 
theory  with  experimental  data,  the  data  will  be  selected  for  com¬ 
binations  of  blade  pitch  P  and  advance  ratio  a  for  which  the 
real  slipstream  effects  are  minimum.  A  method  for  accomplishing 
this  end  will  be  described  below,  and  the  real  slipstream  effects 
upon  the  present  theory  will  be  discussed.  The  purpose  of  the  latter 
discussion  is  to  give  some  insight  into  how  the  predicted  duct 
coefficients  are  influenced,  but  no  attempt  is  made  to  give  quanti¬ 
tative  predictions  of  real  slipstream  effects. 

4.2  Data  Reduction  Method 

The  present  theory  predicts  the  forces  and  moments  on  the  duct 
in  the  presence  of  a  propeller  by  representing  the  propeller  as  a 
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uniformly  loaded  actuator  disk,  and  the  propeller  slipstream  vorticity, 

Y/V^  in  Figure  2,  is  the  only  direct  effect  of  the  propeller  which 

appears  in  the  analysis.  The  purpose  of  this  section  is  to  develop 

expressions  for  Y/V^  in  terms  of  the  propeller  coefficients  which 

are  usually  measured  in  experimental  investigations. 

First,  we  shall  retain  the  assumption  of  an  ideal  actuator 

disk  and  find  y/V«  terms  of  the  propeller  disk  loading  C„ 

°  ^P(D) 

Then  we  shall  include  two  real  propeller  slipstre^un  effects,  the 

swirl  and  the  blade  wake  frictional  losses,  and  derive  a  second 

expression  for  y/V^-  The  purposes  of  the  second  derivation  are 

to  demonstrate:  (1)  a  method  for  selecting  data  which  minimizes  the 

real  slipstream  effects,  and  (2)  the  way  in  which  real  slipstream 

effects  influence  the  predicted  duct  forces  and  moments. 

4.2.1  Evaluation  of  y/V^  for  an  Actuator  Disk 

The  diameter  of  the  actual  stream  tube  flowing  through  the 

ducted  propeller  at  angle  of  attack  changes  continuously  as  shown 

in  Figure  6.  However,  the  duct  forces  and  moments  have  been  pre¬ 
dicted  by  assuming  that  the  slipstream  vorticity  y  could  be  put 
on  the  dashed  cylindrical  surface  of  diameter  O.  In  the  actual 
slipstream,  there  is  a  counter-rotating  swirl  distribution  due  to 
the  cross  flow  of  the  free  stream  over  it  (Fig.  6) .  This  swirl 
distribution  (which  is  generated  by  the  trailing  vortex  filaments 

in  the  duct  wake  dr  in  Fig.  2)  is  already  incorporated  in  the 

% 

present  analysis.  The  kinetic  energy  associated  with  this  swirl 


distribution  causes  an  induced  drag  due  to  lift,  which  appears  in 
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the  calculated  duct  thrust  and  duct  normal  force.  However,  since 
we  have  assumed  that  the  cross  flow  is  small  (V^  sin  «o  « 
we  shall  neglect  the  swirl  velocity  component  in  evaluating  the 
total  pressure  rise  across  the  actuator  disk.  Thus,  we  can  write 
Bernoulli's  equation  ahead  of  and  behind  the  propeller  disk  as 

ip-  f  (167) 

SO  that 


^P(D) 


-  1 


(168) 


Furthermore,  since  in  Figure  6,  we  can  write 

7  ““  Vj  -  Vo  COS 


(169) 


Thus,  by  equating  the  two  expressions  for  Vj  given  by  Equations 
(168)  and  (169),  we  can  evaluate  the  propeller  slipstream  vorticity 
in  terms  of  the  propeller  thrust  coefficient  as 


^  -  '\/cI  -  cos  a  (170) 

^o  V  ^P(D)  ° 

For  the  special  case  of  hovering  flight.  Equation  (15)  replaces 
Equation  (170) . 


4.2.2  Real  slipstream  effects  and  data  selection  procedure 

The  above  evaluation  of  Y/V^  assumes  that  the  slipstream 
is  ideal;  that  is,  the  propeller  is  an  actuator  disk,  which  is 
uniformly  loaded  so  that  it  generates  no  swirl,  the  slipstream 
velocity  is  constant,  and  there  are  ho  frictional  losses 
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(blade  wakes) .  Here  we  shall  evaluate  7/V^  with  these  real  slip¬ 
stream  effects  included  insofar  as  possible  within  the  context  of 
the  present  theory.  To  do  this,  we  shall  first  assume  that  the 
actuator  disk  adds  angular  as  well  as  axial  momentum  to  the  slip¬ 
stream  as  shown  in  Figure  6. 

To  account  for  the  swirl  qualitatively,  we  shall  assume  that 
far  downstream,  at  station  j,  the  slipstream  or  jet  rotates  as  a 
solid  body  with  a  swirl  velocity  component 


The  additional  vorticity  shed  by  the  propeller  to  produce  this 
swirl  is  a  uniform  distribution  of  axial  vortex  filcunents  within 
the  jet  and  a  uniform  distribution  of  filaments  of  opposite  sign 
along  the  jet  boundary.  Within  the  context  of  the  present  theory, 
the  effect  of  this  vorticity  upon  the  duct  is  found  by  placing  the 
slipstream  on  the  dashed  cylindrical  extension  of  the  duct,  as 
indicated  in  Figure  6.  It  can  be  shown  that  when  this  is  done 
the  additional  vorticity  induces  no  velocity  outside  of  the  semi¬ 
infinite  dashed  cylinder.  Thus,  the  equations  for  the  duct  force 
and  moment  coefficients  given  by  the  present  theory  in  terms  of 
Y/Vq  are  not  affected  by  the  additional  vorticity  due  to  solid 
body  rotation  of  the  slipstream. 

The  theory  assumes  that  the  slipstream  vorticity  y ,  associated 
with  the  velocity  component  aligned  with  its  axis,  is  concen¬ 

trated  along  the  jet  surface.  Thus,  within  the  framework  of  this 
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theory,  we  must  retain  the  assumption  that  Vj  is  constant.  How¬ 
ever,  for  single  propellers  operating  at  peak  efficiency,  the  net 
effect  of  propeller  swirl  Vg  is  probably  larger  than  the  non¬ 
uniformity  of  V j .  Hence,  we  can  account  for  the  dominant  real 
effect. 

With  this  representation  of  the  propeller  slipstream,  the 
propeller  disk  loading  Ap  is  no  longer  uniform,  since  the  total 
pressure  in  the  slipstream  is  not  constant.  We  shall  now  assume 
that,  in  addition  to  imparting  axial  and  angular  momentum  and 
energy  to  the  slipstream,  the  propeller  introduces  frictional 
losses  in  the  form  of  blade  wakes.  The  head  loss  H^  associated 
with  these  wakes  caui>es  an  increase  in  propeller  power.  However, 
we  shall  assume  that  the  wedces  mix  rapidly  and  that  the  head  loss 
is  distributed  uniformly  within  the  slipstream.  Thus,  there  is 
no  vorticity  far  downstream  associated  with  H^^  and  no  effect  upon 
the  duct  or  propeller  forces.  The  only  effect  of  H^  then  is  to 
increase  the  propeller  torque  and  power. 

Steady  Flight  at  Angle  of  Attack 

Now  let  us  consider  the  angular  momentum,  the  power,  and  the 
thrust  associated  with  the  slipstreaun  swirl  Vg.  For  a  single 
propeller  and  no  stator,  we  can  equate  the  propeller  torque  to  the 
angular  momentum  of  the  jet,  since  Oj  «  1  (Fig.  6).  Thus,  we 


have 
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2irpV  V  r  ®  ^  d  ^ 

"  i  'in,  ®”> 

o  “'m 
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Aj  -r  -t 

m  m 


m 


where  the  mass  flow  rate  in  the  slipstream  is 


W.  -  pA.V. 
D  D  D 


(172) 


(173) 


For  dual  propellers,  or  propeller-stator  combinations,  t  is  the 
net  torque  on  the  combination. 

The  steady-flow  energy  equation  for  the  stream  tube  in  Figure  6 
(between  stations  far  upstream  and  far  downstream)  can  be  written 
as 


p  +  V  f  p  H,  N 

?  +  -J-J-S-jpV.  d».  -  w .  j  (174) 


where  P  is  the  power  supplied  to  the  propeller  and  H^^/p  is  the 
frictional  head  loss  due  to  blade  wakes.  In  Equation  (174),  we 
have  neglected  the  comparatively  small  swirl  velocity  component  in 
the  jet  which  is  induced  by  the  duct  trailing  vorticity  (see  Fig.  6) 
just  as  we  did  previously  in  calculating  Ap  across  the  actuator 
disk  (see  Eq.  (167)). 


Let  us  now  find  the  variation  of  static  pressure  across  the 
jet  so  that  we  can  evaluate  Equation  (174) .  This  can  conveniently 
be  done  in  a  coordinate  system  which  translates  along  the  jet  with 
velocity  .  In  this  coordinate  system,  we  can  write 


(175) 


Making  use  of  Equation  (171),  we  find  by  integration  that 


Substitution  into  Equation  (174)  then  gives 


and  upon  integration  over  the  jet,  we  find  that 


m 


(177) 


(178) 


where  tcu  is  summed  for  the  propellers  if  there  are  more  than  one. 
Thus,  the  propeller  power  is  increased  by  an  amount 


due  to  swirl  and  blade  wakes 
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Now  let  us  find  the  effect  of  Vg  upon  the  total  thrust  force 
of  the  ducted  propeller,  '^p(p)  ^p(d)’  write  the  momentum 
equation  along  the  duct  axis,  for  a  large  control  surface  surround¬ 
ing  Figure  6,  as 


'^D(P)  ’^P(D) 


-  Wj  cos  Oj  -  cos  "  cos  Oj  I  (pQ  -  p)  dA 


(179) 

Substitution  of  Equation  (176)  and  use  of  cos  »  1  then  gives 


Td(P)  *  ’'p(D)  -  -  ''o  ““  “<.>  -  PX  * 


m 


(180) 


'^D(P)  +  '>’P(D)  "  “o  -  4V 


m 


(181) 


so  that  the  total  thrust  is  reduced  by  an  eunount 


^jV  4V 


due  to  swirl.  Now  we  shall  combine  Equations  (172),  (178),  (180), 
and  (181),  and  put  the  results  in  dimensionless  coefficient  form. 
Thus,  equating  the  two  expressions  for  t  as  given  by  Equations 
(172)  and  (178)  yields  for  a  single  propeller 


1  +  2 


s 


(182) 
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Dlvlslon  of  Equation  (181)  by  (178)  givea  a  propulsive  efficiency 
defined  as 

Vj  *  V 

^  -  cos  o  -  - —  ^ 

,  ^P(f|  ;  ^P(D)  ■  2  ^ (183, 

r  u  2  JJ 

-4-  - 1  +  2  — 4- 

pV„ 


and  the  total  thrust  coefficient  of  the  ducted  propeller  is,  from 
Equation  (181) 

5  4(2)  ^  '^PlDl  .  2  i.  (^4  .  ec  a  -  ^4-  ^  (184, 

^DP  *<><.  ''o  Vo  °  ■‘Vj  J 

where  V  is  the  averaged  axial  velocity  component  through  the 
duct,  given  by 


v-i/v 


dA 


In  addition  to  the  above  set  of  equations,  we  can  write  from 
continuity 


-is  ][o 


(185) 


and,  from  Figure  6,  (where  *  0) 

7  !i 

— •  ai  — ^  —  COS  QL 

V  ^  o 

o  o 


(186) 
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Thus,  we  have  a  set  of  five  equations,  (182)  through  (186),  for 
the  six  dimensionless  unknowns 


V_ 


in  terms  of  three  measured  coefficients 

(a)  Propeller  advance  ratio,  X  -  V^/coR 

(b)  Propulsive  efficiency,  rjp 

(c)  Total  thrust  coefficient,  C„ 

^DP 

We  could  solve  this  set  of  equations  by  using  the  present  theory 
to  evaluate  V/V^  in  terms  of  c/D  and  y/V^.  For  example. 
Figure  3(a)  gives  values  of  V/y  for  hovering  flight  and  for 
c  ••  0.2D.  For  axial  forward  flight,  V  is  simply  increased  by 
The  present  theory  also  gives  the  thrust  ratio 


^D(P)  f  c  JL 
^P(D)  ’  '^o  ’ 


(187) 


so  that  either  of  these  thrust  coefficients  could  be  used  instead 

of  C_  in  Equation  (184).  However,  our  purpose  here  is  to  find 
^DP 

an  independent  expression  for  y/V^  in  terms  of  measured  coefficients 
to  supplement  the  present  theory.  Therefore,  we  shall  drop  H,/pV  ^ 
in  Equation  (183) ,  assiuning  it  to  be  zero,  rather  than  use  the 
present  theory  to  evaluate  V/V^. 
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It  should  be  noted  that  the  coefficients  and  t)- 

^DP 

used  above  are  related  to  the  torque  coefficient  Cq  and  the 
power  coefficient  Cp  as  follows: 


(188) 


C 


P 


T(JD 


(189) 


Axial  Flight 

If  we  now  restrict  the  angle  of  attack  to  zero  and 

eliminate  V_  /V_,  r.  /R,  and  y/V_  with  Equations  (183),  (185), 
m  “'m 

and  (186) ,  respectively,  then  we  can  write  the  remaining  Equa¬ 
tions  (182)  and  (184)  in  the  form 


1 


o: 


and 


(190) 


(191) 


Then,  elimination  of  V/V^  gives  the  expression 
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(192) 


if  H,  ■  a  “0.  From  Equation  (186)  with  a  ■  0,  we  have 
I  o  o’ 


Thus,  Equations  (192)  and  (193)  give  the  propeller  vorticity 

7/V  in  terms  of  three  measured  coefficients  (C„  ,  rjp,  and  ^) 

DP 

for  the  case  in  which  the  angle  of  attack  and  propeller  frictional 
losses  are  zero. 

To  illustrate  the  effect  of  swirl,  we  shall  use  Equation  (192) 

for  y/Vq  with  one  source  of  experimental  data  (Ref.  19)  for 

comparison  with  Equation  (170) .  However,  in  general,  we  shall 

use  Equation  (170)  to  evaluate  Y/V^.  We  shall  do  this  for 

selected  values  of  k  and  ^  which  minimize  the  real  slipstream 

effects.  From  Equations  (170)  and  (183)  it  can  be  seen  that,  for 

fixed  values  of  a  and  C„  (or  V./V) ,  these  effects  (as 

o  ^D(P)  ^  ° 

represented  by  Vg  and  H^^)  are  minimum  when  rjp  is  maximum. 

It  is  noted  in  Equation  (183)  that  both  propeller  swirl  and 
frictional  losses  tend  to  decrease  T]p  as  would  be  e}q>ected  from 
physical  reasoning  for  axial  flow.  It  should  be  recalled  that 
Equation  (183)  is  based  upon  the  following  assumptions: 
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(1)  Vj  is  constant  and  *  0 

(2)  The  jet  slipstrecun  rotates  as  a  solid  body 

(3)  The  frictional  head  loss  is  uniformly  distributed  over 

the  jet 

Using  the  above  results,  we  shall  select  combinations  of  X 

and  p  for  which  t]-  is  a  maximum  for  a  fixed  value  of  C„ 

^  ^P(D) 

Since  ducted  propeller  data  are  usually  presented  in  terms  of  ^ 

and  P,  this  will  be  done  by  use  of  data  plots  as  shown  below  for 

a  fixed  angle  of  attack. 


Sketch  E.-  Example  data  plots. 

It  can  be  seen  that  the  combination  of  and  (3  for  maximum 

Ti„  with  C_  fixed  cannot  be  obtained  by  the  usual  procedure 

^  ^P(D) 

of  maximizing  T]p  with  either  p  or  ^  held  constant. 
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Hovering  Flight 

Carrying  through  the  above  analysis  for  hovering  flight 
(V^  *  0),  with  the  appropriate  definitions  based  on  rather 

then  V^,  we  find  instead  of  Equations  (183)  and  (192) 


and 


=  '^D(P)  1  !p(D). 

%  ■  P 


1  - 


1  + 


% 


^D(P)^  '^P(D)y^ 

^(2 


(194) 


(195) 


It  can  be  seen  that  is  maximum  for  minimum  swirl  and  minimum 

frictional  loss  and  that  could  be  evaluated  from  Equation  (194) 

by  use  of  measured  data  and  Equation  (195)  to  eliminate  V ^ .  How¬ 
ever,  is  the  "compressor  efficiency"  which  is  customarily 

defined  as 


\ 


'‘’d(P)^  * 
p 


/'^P(D) 

V 


h 


for 


0 


(196) 


by  neglecting  swirl  and  evaluating  from  the  Bernoulli  Equa¬ 
tion  (15).  We  shall  retain  this  system  for  evaluating  It 

can  then  be  shown  that 


for 


H, 


% 


0 


(197) 
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where  is  the  usual  figure  of  merit  for  hovering  flight  as 

defined  in  the  List  of  Symbols.  Now,  to  select  hovering  data  for 
which  the  swirl  and  head  loss  are  minimum,  l^he  value  of  p  will 
be  chosen  which  gives  a  maximum  value  of  as  defined  in  Equa¬ 

tion  (197) .  For  no  swirl,  we  have  the  ideal  case  of 


-  1  and 


(198) 


For  a  very  long  duct,  simple  momentum  theory  gives 
and  Equation  (198)  then  gives  for  c  »  D 

tIh  1  and  -  ^|T 


'^P(D)h» 

(199) 


Pitching  or  Plunging  Flight 

For  either  pitching  or  plunging  flight,  the  same  procedure 

of  data  extraction  as  outlined  above  for  steady  flight  at  angle 

of  attack  will  be  used  for  small  pitching  and  plunging  rates.  At 

the  present  time,  however,  no  applicable  dyneunic  data  are  available 

for  comparison  with  the  theory. 

4.3  Thrust  Ratio  in  Hovering  and  Axial  Flight 
Hovering  Flight 

The  division  of  thrust  in  hovering  flight  as  predicted  by 
theory  (Eq.  (17))  is 
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This  expression  is  plotted  in  Figure  7  and  is  compared  with  six 
sources  of  data  chosen  from  Table  I  of  Reference  2.  The  asymptotic 
value  for  large  chord  (Eq.  (199))  is  also  shown.  Equation  (200) 
is  based  upon  the  assumption  that  c  «  D.  However,  this  assumption 
does  not  introduce  much  error  into  the  approximate  expressions  for 
the  velocity  induced  by  the  actuator  disk  if  c  ^  (D/4)  (see  Eq. 
(A.12)).  Furthermore,  the  additional  approximation  that  the  flow 
is  two  dimensional  over  each  chordwise  strip  of  the  duct  is  expected 
to  be  reasonably  accurate  if  c  <  (D/4).  Therefore,  we  expect 
Equation  (200)  to  be  valid  within  this  range  as  indicated  by  the 
solid  line  in  Figure  7. 

Data  bv  Gill 

Shown  in  Figure  7(a) are  data  by  Gill  taken  from  Figures  25, 

26,  and  28  of  Reference  13  for  three  duct  shapes  designed  at  Hiller 
Aircraft  Corp.  These  ducts  were  all  relatively  straight,  thin, 
and  short  and  had  an  inside  diameter  of  24  inches  as  indicated. 

The  ducts  were  used  with  counter-rotating  propellers,  each  having 
three  twisted  and  tapered  blades  with  only  0.04-inch  clearance 
between  the  blade  tips  and  the  inner  surface  of  the  duct.  The 
thrust  ratio  was  found  to  depend  only  very  slightly  on  the  mean 
blade  pitch  setting  p,  as  shown  in  Figure  7(a).  With  independent 
variation  of  p,  neither  the  figure  of  merit 


“h  ^ 


D(P), 


+  T 


P(D), 


/^D(P).  '^P(D) 

V  2pA 


P 


(201) 
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nor  the  compressor  efficiency 


^D(P), 


+  T 


P(D), 


%  = 


P(D), 

2pA 


(202) 


was  a  maximum  when  the  thrust  ratio  '^D(p)/'^p(D)  maximum,  as 

can  be  seen  from  the  data  of  Figure  7(a).  This  behavior  is  believed 
to  be  associated  with  flow  separation  from  the  duct  leading  edge, 
which  was  actually  observed  during  hovering  flight  for  all  the 
duct  shapes  except  for  one  having  a  large  inlet  lip  radius . 

Because  of  this  fact,  it  is  not  surprising  that  the  measured  thrust 
ratio  was  only  about  half  of  that  calculated.  That  is,  the  thrust 
carried  by  the  duct  was  evidently  greatly  diminished  by  flow  separa¬ 
tion  from  the  duct  leading  edge. 

Data  bv  Horn 

Shown  in  Figure  7(b)  is  a  datum  point  by  Horn  taken  from  Fig¬ 
ure  13  of  Reference  17 .  The  configuration  was  of  the  shape  indicated 
and  was  tested  in  water.  A  single  propeller  with  four  blades  was 
used  without  a  stator,  so  that  swirl  was  introduced  into  the  slip- 

t 

stream.  Furthermore,  the  duct  is  too '.long  for  comparison  with  the 
theory,  and  the  blade-pitch  angle  was  r^ot  varied.  Propeller  loca- 
tion  was  varied  in  a  still  longer  duct  (which  is  therefore  not 
shown  in  the  figure)  and  there  was  little  effect  of  propeller 


location  upon  either  ^p(D)  '^D(P)  ‘ 

^  h  ^  h 
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Data  bv  Moaer 

Also  shown  in  Figure  7  (a)  are  data  by  Moser,  taken  from  Fig¬ 
ure  12(b)  of  Reference  18.  The  duct  shape  was  varied  by  adding  a 
sheet-metal  ring  after  the  duct,  as  indicated,  such  that 

0*33  1  i  0.89 

It  was  found  that  the  thrust  ratio  was  virtually  invariant  with 
duct  length  as  indicated  by  the  arrow  in  the  figure.  The  single 
propeller  had  two  blades.  The  blade  pitch  was  not  varied,  and 
stators  were  not  used. 

Data  bv  Kruger 

Finally  shown  in  Figure  7(b) are  data  by  Kruger,  taken  from 
Figure  21  of  Reference  19  for  one  of  fifteen  shroud  shapes  tested 
with  a  relatively  large  centerbody.  The  duct  shape  used  here  had 
the  smallest  camber,  thickness,  and  chord-to-diameter  ratio  of 
all  the  ducts  tested  therein,  A  single  propeller  with  eight 
blades  and  no  stator  was  used  with  this  duct,  and  the  blade-pitch 
setting  p  was  varied  independently.  It  can  be  seen  that  the 
thrust  ratio  in  this  case  was  very  dependent  upon  p  and  again 
was  not  maximum  for  peak  efficiency  (either  or  q^) .  These 

peaks  occurred  at  the  lowest  blade-pitch  angle  which  was  tested, 
namely  15°.  However,  it  can  be  seen  from  the  measured  values  of 
duct  and  propeller  thrust  in  Figure  21  of  Reference  19  that  both 
the  propeller  and  the  duct  were  stalled  at  the  higher  pitch  angles 
so  that  the  measured  thrust  ratio  at  high  ^  is  not  actually 


-88 


comparable  with  the  theory.  For  a  longer,  thicker,  converging 
shroud,  hovering  thrust  ratios  as  great  as  two  were  obtained. 
However,  the  duct  length,  c  ■  0.625D,  and  convergence  were  too 
large  for  confer ison  with  the  present  theory. 


Data  bv  Grose 

Shown  in  Figure  7  (c)Eu:e  data  by  Grose,  taken  from  Figures  18 

to  21  of  Reference  14.  The  shroud  shown  in  Figure  7(c)  (Grose's 

"high-speed"  shroud),  was  relatively  thin  and  straight,  as  indi-, 

cated;  however,  it  was  rather  long  for  comparison  with  the  present 

theory  (c  ■  0.475D).  The  flow  was  observed  to  separate  from  the 

leading  edge  and  to  reattach  near  the  trailing  edge  of  this  duct 

during  hovering.  The  variation  of  thrust  ratio  /^p(D) 

with  blade  pitch  is  shown  in  Figure  7(c).  In  this  case  Mj^, 

and  Tjj/pj  /Tp/jjj  were  all  maximum  for  p  •  15°.  A  second 

^h  ^h 

"static"  shroud,  designed  with  a  large  radius  inlet  lip  to  prevent 
separation,  had  thrust  ratios  as  large  as  unity.  Both  ducts  were 
tested  with  a  single  propeller  having  four  blades  and  no  stators. 
Propeller  tip  speeds  were  between  500  and  1200  feet  per  second 
during  these  tests . 


Data  bv  Platt 

Also  shown  in  Figure  7(c)  are  data  by  Platt,  teOcen  from  Figure  19 
of  Reference  20  for  the  shortest  of  three  ducts  tested.  This  duct 
was  tested  with  counter-rotating  propellers  with  blade  tip  speeds 
between  400  and  600  feet  per  second.  The  front  propeller  had  five 
blades  and  the  rear  seven  blades.  The  variation  of  thrust  ratio 
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with  blade  pitch  is  indicated  in  Figure  7(c).  Here  and 

’^D(P)x  /Tp(D)  were  all  highest  for  the  highest  pitch  angle  tested, 
namely  p  ■  45°.  It  will  be  noted  that  the  values  of  ^^(p)  /''^p(d) 
measured  by  Platt  were  far  in  excess  of  the  values  predicted  either 
by  the  present  theory  or  by  simple  momentum  theory.  However,  neither 
of  these  theories  is  expected  to  be  valid  for  this  intermediate  range 
of  c/D  (c  ■  0.67D).  Thus,  the  need  for  a  theory  in  the  intermediate 
range  of  chord-to-diameter  ratio  is  readily  apparent.® 

Flow  separation  and  high  noise  level  were  observed  by  Platt  only 
at  low  propeller  speed  with  this  duct.  At  the  higher  rotational  speeds 
tested  (corresponding  to  the  data  shown  in  Fig.  7(c)),  the  flow  was 
attached  and  the  duct  thrust  was  about  twice  as  great  as  for  separated 
flow.  This  configuration  is  the  only  one  known  for  which  the  duct 
was  essentially  straight  and  thin  and  for  which  the  flow  was  known  to 
be  unseparated  in  hovering.  The  fact  that  the  Grose  duct  was  always 
stalled  and  the  Platt  duct  generally  unstalled  during  hovering  tests 
apparently  was  not  due  to  Reynolds  number,  because  the  maximum  Reynolds 
numbers  based  on  tip  speed  and  chord  length  were  nearly  equal  for  the 
two  cases . 

Axial  Flight 

The  thrust  ratio  for  axial  flight  can  be  written  (from  Eqs.  (11), 
(168),  (169),  and  (200))  as 

(203) 


P(D) 


1  +  2  ^ 


1+2 


^Subsequent  analysis  has  shown  that  for  hovering  flight  the  theoretical 
thrust  ratio  is  greater  than  unity  for  the  intermediate  range  of  c/D. 
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Thus,  the  thrust  ratio  in  axial  flow  is  evidently  given  by  the 
value  for  hovering  divided  by  1  +  2(V^/y).  Furthermore,  the  ratio 
Vq/Y  is  a  function  only  of  the  actuator  disk  loading  when  there 
are  no  real  slipstream  effects.  That  is,  Equations  (169)  and  (170) 
give  (for  ■  0) 


JL 

^o 


+  1 


P(D) 


+  1 


(204) 


Thus,  Equations  (203)  and  (204)  give  the  predicted  thrust  ratio  in 
terms  of  c/D  and  the  measured  propeller  thrust  coefficient.  This 
predicted  value  will  now  be  compared  with  five  sources  of  measured 
data.  The  data  will  be  restricted,  whenever  possible,  to  those 
combinations  of  propeller  advance  ratio  k  *  V^ouR,  and  blade-pitch 
setting  p  which  give  maximum  propulsive  "efficiency," 

Ip  -  ^  pHP)J  ''o  (205) 


with  fixed  propeller  loading  C„  .  This  combination  should 

^P(D) 

minimize  the  real  slipstream  effects  (see  Section  4*2)  which  have 

been  assumed  to  be  zero  in  the  present  theory.  Thus,  the  predicted 

value  of  propulsive  efficiency  (Eq.  (183))  with  V_  -  H,  "  0)  is 

m 

the  "ideal  Froude  efficiency." 


^P  “ 


(206) 
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Data  bv  Gill 

Of  the  three  Hiller  ducts  shown  in  Figure  7(a),  only  DPS 

4  3 

was  tested  in  axial  flight.  Measured  values  of  T)p  and  values  of 
Vj/V^  calculated  from  Equation  (204)  are  plotted  in  Figure  8, 
using  data  from  Figures  178  and  179  of  Reference  13,  and  the  method 
of  data  selection  is  indicated  in  Figure  8  for  values  of  Vj/V^ 
of  1.75,  2.0,  and  3.5. 

For  the  values  of  V^/V^  indicated  in  Table  VII,  this  pro¬ 
cedure  (A  and  ^  in  each  case  selected  for  the  highest  value  of 
Tjp)  gives  the  results  shown  therein.  The  values  of  V^/V^  in 
Figure  8  differ  from  those  in  Table  VII  because  the  former  v«re 
chosen  to  illustrate  that  p  varies  with  V^/V^  for  maximum  rjp. 
The  measured  values  of  T)p  in  Table  VII  are  considerably  smaller 
than  the  predicted  ideal  Froude  efficiency  (Eq.  (206)).  Therefore, 
real  slipstream  effects  are  probably  not  negligible  as  assumed  in 
the  above  comparison  of  data  and  theory. 

It  was  not  reported  in  Reference  13  whether  this  duct  was 
stalled  at  the  above  test  conditions,  but  it  was  observed  to  be 
stalled  during  hovering.  Therefore,  the  low  values  of  measured 
thrust  ratio  as  compared  with  theory  could  be  due  to  flow  separation 
and/or  real  slipstream  effects. 

Data  bv  Horn 

The  duct  shape  of  Horn,  shown  in  Figure  7(b),  was  also  tested 
in  axial  flow,  and  the  values  in  Table  VIII  were  obtained  from 
Figure  13,  Reference  17,  by  the  method  described  above. 
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As  indicated  in  Figure  7(b),  this  duct  is  actually  too  long  for 
accurate  prediction  of  the  thrust  ratio  by  the  present  theory.  A 
single  propeller  was.  tested  with  3  fixed.  The  axial  flight  data 
indicate  a  wide  variation  of  thrust  ratio  with  advance  ratio.  At 
the  advance  ratio  for  which  T]p  was  highest,  there  is  little 
difference  between  the  small  values  of  measured  and  predicted 
thrust  ratio.  Real  slipstream  effects  were  probably  large,  since 
T]p  as  measured  was  much  smaller  than  the  predicted  ideal  Froude 
efficiency. 

Data  bv  Allen 

In  Reference  21,  data  are  presented  for  axial-flow  tests  of 
two-bladed  propellers  with  and  without  a  duct.  Three  different 
propellers  were  used  which  had  pitch  settings,  p  >  20°,  25°,  and 
30°.  The  duct  shape,  propeller  location,  and  tip  clearance  are 
indicated  in  Sketch  F,  below. 


Sketch  F.-  Ducted-propeller  configuration 
tested  by  Allen  (Ref.  21). 
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It  can  be  seen  that  more  duct  divergence  and  "internal  diffusion" 
were  incorporated  into  this  configuration  than  in  those  of  Figure  7. 

From  a  ploc  of  measured  data  (taken  from  Fig.  2,  Ref.  21) 
similar  to  Figure  8,  the  values  in  Table  IX  were  obtained  for  maximum 

at  the  selected  values  of  C„ 

^  ^P(D) 

Both  the  measured  pressure  distribution  over  this  duct  surface 

and  the  good  agreement  between  the  measured  values  of  duct  thrust 
and  those  predicted  by  the  approximate  potential  flow  analysis 
(of  Ref.  21) ,  indicate  that  this  duct  was  not  appreciably  stalled 
during  high-speed  forward  flight  with  any  of  the  three  propellers 
tested.  However,  the  measured  value  of  gradually  increased 

and  then  suddenly  dropped  appreciably  when  V^/nD  was  decreased 
to  values  slightly  less  than  those  corresponding  to  the  data  used 
above.  It  is  therefore  concluded  that  the  duct  was  not  stalled 
for  the  data  used  above.  At  high  advance  ratios,  beceune 

negative  due  to  the  large  external  duct  divergence  and  friction 
drag.  Negative  values  of  ^^(p)  also  were  measured  without  the 
propeller.  It  is  observed  that  high  values  of  T)p  were  attained 
and  that  the  measured  thrust  ratio  for  p  -  25°  is  greater  than 
the  predicted  ratio  (for  a  straight  duct) .  It  seems  reason2d)le 
to  attribute  the  latter  fact  mainly  to  internal  diffusion. 

Data  bv  Kruaer 

The  ducted  propeller  of  Kruger  shown  in  Figure  7b  was  also 
tested  in  axial  flow,  and  the  blade  pitch  setting  was  varied  inde¬ 
pendently  over  the  range  15°  <  ^  <  55°.  The  procedure  described 
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above  was  employed  for  selecting  "K  and  p  for  maximum  T)p  at 

fixed  C_  ,  and  the  values  in  Table  X  were  obtained  from  a 
^P(D) 

data  plot  similar  to  Figure  8.  (Data  from  Fig.  21,  Ref.  19,  were 
used  for  this  purpose) . 

It  can  be  seen  that  at  the  highest  value  of  C-  the 

P(D) 

measured  value  of  T)p  was  not  much  smaller  than  the  ideal  Froude 

efficiency  and  the  agreement  between  the  measured  and  predicted 

value  of  thrust  ratio  is  fair. 

However,  comparatively  high  values  of  3  and  were  used 

during  these  tests.  From  the  data  (Fig.  21,  Ref.  19),  it  appears 

that  stalling  of  both  the  duct  and  the  propeller  occurred  at  low 

advance  ratios,  but  that  neither  was  stalled  at  the  advance  ratios 

in  Table  X.  Therefore,  at  the  lower  values  of  C-  the  poor 

P(D) 

agreement  between  the  theory  and  data  for  rjp  and  ''^i5(p)/^p(D) 
is  attributed  mainly  to  real  slipstream  effects.  However,  a  small  part 
of  the  duct  thrust  was  generated  by  interference  with  the  compara¬ 
tively  large  centerbody.  This  effect  was  measured  (Fig.  7,  Ref.  19) 
and  found  to  give  ^d(h)  "  0.063  Aq^. 

To  illustrate  that  the  above  data  reduction  procedure  gives 
better  overall  agreement  between  theory  and  experiment,  we  shall 
now  use  the  more  conventional  method  where  for  each  pitch  setting 
the  advance  ratio  ^  which  gives  maximum  propulsive  efficiency 
is  chosen.  The  results  are  given  belowt 
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1.0 

.039 

.52 
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Upon  comparison  with  Table  X  it  i,s  evident  that  the  overall  agree¬ 
ment  for  thrust  ratio  is  not  as  good  by  this  method. 


Ve  shall  now  attempt  to  account  for  the  real  slipstreeun  effects 

in  these  data  by  using  the  procedure  described  in  Section  4.2.2. 

With  this  method,  we  choose  X  and  p  for  highest  T]p  for  each 

fixed  value  of  C-  as  was  done  for  Table  X.  However,  we 

^P(D) 

shall  use  Equations  (192)  and  (193)  to  estimate  with  pro¬ 

peller  swirl  present  in  the  slipstreeun,  rather  than  using  Equa¬ 
tion  (170)  which  is  valid  only  for  an  ideal  slipstream  with  no 
swirl.  Using  this  methbd,  we  find  the  values  given  below  for  two 

values  of  C,-  which  correspond  to  the  first  and  fourth  values 

’^P(D) 

in  Ted>le  X. 


'■  Vain'S-"  ■ 
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(juR 

^T 

^P(D) 

X 

^o 

^D(P) 

’’P(D) 

^D(P)  . 

^P(D) 

(pred. 

0-75 

1  +  2(V^/7) 

)(Eq.  (203)) 

15 

0.14 

3 

2.05 

0.3 

0.38 

25 

.32 

0.5 

0.75 

.6 

.21 

It  can  be  seen  from  these  results  that  the  comparison  for  thrust 

ratio  is  improved  by  including  the  effect  of  propeller  swirl  in 

evaluating  t/V  .  The  experimental  values  of  T^p,  C„  ,  and  ^ 

DP 

were  used  in  evaluating  from  Equations  (192)  and  (193). 

The  relatively  poor  agreement  for  the  lower  value  of  C_  may 

^P(D) 

be  due  to  frictional  head  loss  which  is  not  included  in  Equation  (192) , 
since  we  had  set  ■■  0. 

Data  bv  Grose 

The  ducted  propeller  of  Grose  shown  in  Figure  7(c)  was  also 

tested  in  axial  flow  with  independent  variation  of  blade  pitch. 

Using  the  data  reduction  method  illustrated  in  Figure  8  (for 

maximum  Tip  at  fixed  C„  ) ,  we  find  the  values  in  Table  XI 
^  ^P(D) 

by  using  data  from  Figures  33  through  37  of  Reference  14.  It 
can  be  seen  that  relatively  high  values  of  T)p  were  attained  in 
these  tests,  so  that  real  slipstream  effects  were  probably  com¬ 
paratively  small. 

Although  this  duct  was  apparently  unstalled  in  axial  flow 
(from  the  discussion  in  Ref.  14),  particularly  at  the  higher 
advance  ratios,  agreement  between  theory  and  experiment  for  the 
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thrust  ratio  cannot  be  expected  because  the  duct  was  too  long 
(c/D  -  0.475)  for  accurate  prediction  with  the  present  theory. 
Furthermore,  the  frictional  drag  of  the  duct  was  relatively  large 
for  these  tests  as  compared  with  the  duct  thrust.  In  fact,  it 
was  found  that  for  each  blade  pitch  setting  tested,  except  ^ 
there  was  a  net  drag  force  on  the  duct  at  the  advance  ratio  for 
which  T)p  was  maximum.  Using  the  data  at  peak  efficiency  for 
3  «  17°,  we  find  the  following  values. 


OlQ 
>  1  C 

^T 

’^P(D) 

y/^o 

‘^D(P)/^P(D) 

'’^D(P)  ..  0^070 

^P 

Ip  *  v/ 

(pred.) 

*P(D) 

(pred. ) 

17 

0.66 

0.41 

0.19 

0.0049 

0.061 

0.72 

0.91 

Thus,  much  smaller  than  predicted,  in  fact  generally 

negative,  by  the  data  reduction  procedure  using  T]p  for 

max 

fixed  p. 

4.4  Forces  and  Moments  at  Angle  of  Attack 

Here  we  shall  compare  the  present  theory  with  all  the  available 

data  for  the  static  stability  derivatives  of  ducted  propellers 

having  relatively  short,  straight,  thin  ducts.  There  are  evidently 

two  sources  of  such  data:  Reference  14  for  only  (C„  )  at 

\  “a/D(P) 

relatively  low  angles  of  attack  (0  <  <  6.3°),  and  Reference  13 

for  both  {c„  ^  and  f \  at  relatively  large  angles 

V  Na^(P)  \  ”*a/D(P) 

of  attack  (50°  <  <  80®) . 
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The  static  derivatives  of  the  duct  forces  and  moments  as 

predicted  by  the  present  theory  will  be  compared  with  these  two 

sources  of  data  as  selected  by  the  technique  used  in  Section  4.3 

above.  That  is,  combinations  of  advance  ratio  and  blade  pitch 

will  be  chosen  which  give  highest  tIq  at  fixed  values  of  C_ 

^  ^P(D) 

Data  bv  Grose 

The  comparison  of  the  Grose  data  will  be  made  for  the  "  high 
speed"  duct  (shown  in  Fig.  7(c))  at  the  highest  angle  of  attack 
tested,  namely  6.3°.  The  data  are  from  Figures  69  through  80  of 
Reference  14. 

The  predicted  value  of  the  duct  lift  force,  for  «  1, 
can  be  expressed  as 

S(P)  “  *^D(P)  *  ’’d(P)  “o 


which  gives  upon  substitution  of  Equations  (9)  and  (11)  with 
c/D  -  0.475 


S(P) 


a 


D(P)  ‘^o^^o 


3.42  +  1.47 


+  0.009 


(208) 


The  value  of  will  be  obtained  from  Equation  (170)  using 


the  measured  value  of  C_ 

P(D) 

maximum  ti-,  at  fixed  C„ 

^  ^P(D) 

gives  the  values  in  Table  XII. 


.  Choosing  values  of  ^  and  k  for 
from  a  data  plot  similar  to  Figure  8 


The  duct  was  apparently  unstalled  during  the  test  conditions 


used  above.  All  of  the  data  used  here  were  with  a  free-streeun 
Mach  number  of  0.2,  the  lowest  value  which  was  tested.  The 
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agreement  between  theory  and  experiment  for  the  duct  lift  curve 
slope  is  excellent,  despite  the  rather  large  value  of  c/D. 

Data  bv  Gill 

The  three  ducts  of  Gill  shown  in  Figure  7  (a)  were  tested  at 
large  angles  of  attack.  However,  only  the  data  for  the  shortest 
duct  (D  P^S)  will  be  compared  with  theory  for  nonaxial  flow. 

All  of  the  ducts  were  probably  stalled  at  the  large  angles  of 
attack  tested. 

The  predicted  coefficients  of  this  duct  (from  Eqs.  ,v(9) ,  i.dl) , 
and  (13),  with  c  ■  0.15D)  are 

(S)d(P)  “  “o  “o  <209) 


'^P(D) 


0.680  +  1.25  sin® 


(210) 


('m)D(P)  -  1-30  ‘in  “o  <211) 

and  from  Equations  (169)  and  (170),  we  have 

^  ■  t  *  co=  %  ■  V%,D)  "  ^ 

Using  the  data  reduction  method  of  Section  4.2.2  and  the 

measured  data  from  Figures  58,  60,  61,  196,  198,  and  199  of 

Reference  13,  we  obtain  the  values  in  Table  XIII  for  two  angles 

of  attack,  50°  and  ■ 

»  o  o 


80°. 
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It  can  be  seen  that  the  measured  values  of  thrust  ratio  and 
duct  pitching  moment  are  much  smaller  than  predicted.  This  is 
attributed  to  flow  separation  at  the  duct  leading  edge,  since  both 
Td(p)  and  ^(p)  depend  primarily  on  leading-edge  suction  (see 
Eg.  (13)).  The  agreement  is  considerably  better  for  (^n)d(P) 
because  it  is  less  sensitive  to  leading-edge  separation  than  is 
duct  thrust.  The  experimental  coefficients  (^n)d(P)  from 

53  to  86  percent  of  the  predicted  values . 

4.5  Sumitiary  and  Conclusions  for  Comparison  of  Data  and  Theory 

For  hovering  flight,  the  predicted  values  of  the  duct-to- 
propeller  thrust  ratio  are,  in  general,  eibout  twice  the  measured 
values.  This  is  probably  due  mainly  to  flow  separated  from  the 
duct  leading  edges.  However,  in  the  case  of  the  duct  tested  by 
Platt  (Ref.  20),  which  was  comparatively  long,  the  measured  values 
of  thrust  ratio  were  much  higher  than  predicted  either  by  the 
present  theory  (for  small  c/D)  or  by  simple  momentum  theory  (for 
large  c/D) .  Additional  uncertainty  is  introduced  into  the  results 
by  real  slipstream  effects  as  described  in  Section  4.2.2. 

For  axial  flight,  the  data  of  Gill  (Ref.  13)  and  Grose  (Ref.  14) 
are  probably  most  comparable  with  the  present  theory,  even  through 
flow  separation  probeds ly  occurred  for  the  duct  used  by  Gill;  and 
the  duct  used  by  Grose  was  actually  too  long  for  accurate  comparison 
with  the  present  theory.  The  measured  thrust  ratios  for  these  ducts 
are  about  half  the  predicted  values. 
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The  data  reduction  procedure  described  in  Section  4.2.2  to 
account  for  real  slipstream  effects  was  compared  with  the  usual 
procedure,  and  it  was  shown  that  the  comparison  of  theory  and  the 
data  by  Kruger  (Ref.  17)  was  considerably  improved  when  the  real 
effects  were  included  in  the  theory. 

Excellent  agreement  is  obtained  between  measured  and  predicted 
values  of  duct  lift  force  at  small  angles  of  attack  (a^  <  6°)  using 
the  data  by  Grose  (Ref.  14),  despite  the  large  chord/diameter  ratio 
(c/D  -  0.475).  Only  large  angles  of  attack  were  tested  for  the 
duct  used  by  Gill  (Ref.  13),  and  this  duct  was  probably  stalled. 
Hence,  the  duct  thrust  and  pitching  moment  were  much  smaller  than 
the  calculated  values.  Nevertheless,  the  normal  force  was  more 
than  half  the  computed  value. 

The  following  conclusions  are  drawn  from  the  comparison  of 
data  with  the  present  theory: 

(1)  For  the  only  source  of  data  for  (C„  \  at  small 

V  ”ajD(P) 

angles  of  attack  with  a  thin  duct  to  which  the  flow  was  attached 
(Ref.  14) ,  this  derivative  was  in  excellent  agreement  with  the 
present  theory.  Hence,  the  dominant  contribution  of  the  ducted 
propeller  to  the  dyneunic  stability  of  a  torpedo-like  configuration 
is  well  predicted  by  the  present  theory  (see  Section  3.5). 

(2)  At  high  angles  of  attack  with  flow  separation  from  the 

duct  (Ref.  13),  experimental  duct  coefficients  (^ci)d(p)  from 

53  to  86  percent  of  the  predicted  values.  The  measured  values  of 
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duct  thrust  and  pitching  moment  were  much  smaller  then  predicted. 
This  is  due  in  part  to  real  slipstreaum  effects  but  probably  more 
to  effects  of  flow  separation  upon  duct  leading-edge  suction. 

(3)  There  is  an  urgent  need  for  experimental  data  on  ducted 
propellers  which  do  not  exhibit  leading-edge  separation. 


5.  OUTLINE  OF  TEST  PROGRAM 


It  is  evident  from  the  previous  comparison  of  ducted  propeller 
data  with  the  present  theory  that  additional  data  are  needed  for 
all  flight  conditions  with  short  ducts  and  without  flow  separation 
from  the  duct  leading  edge.  There  are  no  data  whatever  for  the 
dynamic  stability  derivatives  of  Isolated  ducted  propellers.  Fur¬ 
thermore,  there  is  an  apparent  need  for  parametric  experimental 
studies  of  such  design  variables  as  propeller  position  in  the  duct, 
duct  chord-to-diameter  ratio,  duct  camber,  etc.  Basic  flow  studies 
should  be  directed  toward  prevention  of  flow  separation  from  the 
duct  leading  edge  and  direct  measurement  of  real  slipstream  effects. 

The  following  experimental  program  is  specifically  designed 
to  provide  ducted  propeller  stability  derivatives  for  comparison 
with  the  present  theory.  A  model  configuration  is  suggested, 
specific  tests  are  outlined  for  each  flight  condition,  and  the 
stability  derivatives  are  predicted  for  the  suggested  configuration. 

5.1  Model  Design 

The  ducted  propeller  shown  in  Figure  9  is  suggested  for  ob¬ 
taining  data  comparable  with  the  present  theory.  The  inner  dia¬ 
meter  is  constant  so  that  tip  clearance  does  not  change  with  pro¬ 
peller  location.  The  radius  of  the  inlet  lip  is  relatively  large 
and  is  provided  with  boundary- layer  control  to  prevent  leading- 
edge  flow  separation.  Boundary- layer  control  is  considered 
essential  to  prevent  flow  separation  from  the  duct  for  all  the 
test  conditions.  Suction  through  either  a  slotted  or  porous  tube 
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seems  suitable  for  boundary- layer  control,  and  the  precise  shape  *  ^ 

of  the  leading  edge  may  have  to  be  investigated.  Counter-rotating, 
variable-pitch  propellers  should  probably  be  used  to  minimize  swirl. 

The  Hiller  rig  (at  DTMB)  having  three-bladed  propellers  with  a 
diameter  of  24  inches  and  a  spacing  of  1.5  inches  (Ref.  13)  might 
be  used  with  a  modified  duct.  The  propellers  should  be  movable  in 
the  duct  to  at  least  two  well-separated  positions  (fore  and  aft)  . 

If  the  propeller  axial  location  is  varied  in  the  duct,  the 

inflow  velocity  distribution  to  the  propeller  disk  will  vary  with  j 

axial  position  of  the  propeller.  The  present  analysis  assumes  a 

uniform  disk  loading,  and  it  may  be  necessary  to  vary  both  blade  .1 

setting  and  blade  twist  to  obtain  a  uniform  disk  loading  experi-  -j 

mentally.  On  the  other  hand,  changes  in  blade  setting  alone  may 

give  a  good  enough  approximation  to  uniform  loading.  In  order  to  1 

investigate  this  point,  we  have  calculated  the  radial  profile  of 

axial  velocity  at  the  entry,  exit,  and  central  duct  planes  for 

c  ■  0.2D  and  a  “  0.  Tl  results  (Fig.  3(a))  show  the  axial 

o  1 

i 

velocity  profiles  induced  by  the  bound  vorticity  on  the  duct  and 
the  trailing  vorticity  shed  from  the  disk.  It  can  be  seen  that  the 
profiles  are  relatively  flat,  except  in  the  range  0.8  <  r/R  <  1. 

This  is  particularly  true  when  the  axial  component  of  the  free 
stream  velocity,  cos  a^,  is  added  to  these  profiles  to  obtain 
the  total  axial  velocity  profile.  Thus,  the  existing  dual  pro¬ 
pellers  used  with  the  Hiller  rig  seem  adequate  for  the  study  out-  j 

lined  here.  (These  propellers  were  designed  to  operate  in  a 

I 
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unlfomi  flow.)  Furthermore,  actual  duct  thickness  and  boundary- 
layer  effects  tend  to  decrease  the  velocity  near  the  duct  leading 
edge  from  the  computed  value,  which  is  singular  as  shown  in  Figure 
3(a)  . 

In  the  Hiller  rig,  duct  forces  and  total  configuration  forces 
were  measured  (Ref.  13) .  By  subtraction,  the  forces  on  the  propeller 
plus  motor  housing  were  obtained.  By  use  of  a  dummy  motor,  the 
aerodynamic  force  on  the  driving  motor  housing  was  shown  indirectly 
to  be  small.  However,  it  would  be  desirable  to  measure  this  force 
directly,  particularly  at  angle  of  attack,  by  use  of  a  separate 
fairing  or  wind  shield  over  the  motor.  In  this  way,  any  appreciable 
disturbance  to  the  slipstream  and/or  error  in  the  measured  aero¬ 
dynamic  reaction  on  the  propeller  could  be  measured  directly. 


5.2  Hovering  Tests 

The  following  specific  procedures  are  suggested  for  the  hover¬ 
ing  case:  Use  only  the  highest  practical  propeller  speed  (and 
Reynolds  number) .  Observe  flow  separation  and  measure  minimum 
boundary- layer-control  flow  required  to  prevent  it.  Using  this 

amount  of  boundary- layer-control,  measure  '^£|(p)  »  "^pdj)  • 

h  h 


Compare  thrust  ratio  with  theory  by  selecting  data  as  described  in 
Section  4.2.  For  each  value  of  traverse  the  slipstreeun  and 
measure  the  magnitude  and  direction  of  the  time  average-velocity 
vector  to  determine  whether  the  swirl  is  actually  a  minimum  When 
is  maximum.  The  traverse  probably  should  be  about  one 
diameter  downstream  of  the  duct  exit  plane  in  order  to  minimize 
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both  blade  Interference  with  the  probe  and  excessive  mixing  of  the 
slipstream  with  entrained  fluid.  Total  pressure  should  also  be 
measured  to  determine  whether  it  is  nearly  constant  across  the 
wake,  as  assumed.  Use  at  least  two  propeller  positions  in  the 
duct ,  varying  (3  at  each . 


5.3  Axial  Flow  Tests 

If,  as  predicted  by  theory,  there  is  found  to  be  little  effect 

of  propeller  axial  position  on  /Tp^jjj  for  hovering 

h  ^ 

(measured  at  optimum  for  each  position)  ,  then  only  one  central 

position  might  be  used  in  the  axial  flow  tests.  However,  at  least 

one  check  case  should  be  made  to  determine  whether  propeller 

position  also  has  little  effect  in  forward  flight,  as  predicted. 

Repeat  the  procedure  described  for  a  sufficient  number  of  advance 

ratios  to  establish  the  values  of  C_  0»^)  >  ti„  (PA)  »  and 

^P(D)  P 

Td(p)/Tp(d)  at  maximum  T^p.  Compare  Tj,  (p)/Tp  with  theory. 

Use  the  duct  alone  to  see  if  Tj^  =  0  as  predicted  for  zero  duct 
thickness  and  to  determine  whether  the  use  of  boundary-layer  control 
actually  has  little  effect  upon  T^^  (aside  from  preventing  flow 
separation  in  the  presence  of  the  propeller) . 


5.4  Tests  at  Angle  of  Attack 

Repeat  the  procedure  in  Section  5.3  for  a  sufficient  number 
of  values  of  in  the  range  0  <  <  90°  to  determine  the 

effect  of  on  the  aerodynamic  reactions  (including  the  moment 

and  normal  force  on  the  duct  and  on  the  propellers) .  Use  the  duct 
alone  for  comparison  with  theory  (at  C„  ■  0)  and  also  for 

S(I» 
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comparison  with  results  obtained  for  the  value  of  ^  for  which  the 
propeller  is  free-wheeling.  The  latter  comparison  will  indicate 
the  effect  of  propeller  losses.  Traverse  the  slipstream  to  deter¬ 
mine  swirl,  variation  in  total  pressure,  and  deflection  of  the  slip¬ 
stream  boundary  by  the  free  stream  (if  feasible  in  the  presence  of 
mixing) .  For  each  run,  use  sufficient  boundary-layer  control  to 
prevent  flow  separation  everywhere  on  the  duct.  Compare  the  static 
stability  derivatives  and  thrust  ratios  with  the  theory  for  which 

each  of  these  is  given  as  a  function  of  c/D,  C-  ,  and  a 

^P(D)  ° 

at  optimum  r]p. 

Many  of  the  tests  suggested  herein  are  actually  being  planned 
at  DTMB  and  should  be  coordinated  with  that  facility. 

5.5  Pitching  (q)  Tests 

In  evaluating  stability  derivatives  for  pure  pitch  (as  on  a 
rotating  arm)  use  the  same  blade  pitch  settings  which  gave  optimum 
T)p  in  Section  5.4  for  each  value  of  Determine  the  boundary- 

layer  control  required  to  prevent  flow  separation  for  each  run  (or 
determine  this  in  Section  5.4  versus  Reynolds  number  and  use  the 
results  here).  Use  maximum  arm  radius  and  speed  to  maximize  Reynolds 
number  and  minimize  boundary-layer  control  requirement.  Use  the 
duct  alone  to  obtain  additional  Information  (as  in  Section  5.4)  at 
angles  of  attack  for  which  separation  can  be  prevented.  Measure 
the  aerodynamic  reactions  on  the  duct  and  propellers,  propeller 
power  P  and  propeller  rotational  speed  a>.  Compare  with 
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theory  which  gives  these  (as  function  of  c/t), 
at  optimum  'n„)  •  If  a  rotating  arm  of  radius 

Sr 


%  (D)’ 

R  is  used,  then 

cL 


q  -  ^  (213) 

^'a 

and  the  dimensionless  pitch  parameter  is 

Ea  .  ^  (214) 

o  a 

Several  values  of  this  ratio  should  be  tested  in  order  to  deter¬ 
mine  the  maximum  value  to  which  each  aerodynamic  reaction  remains 
proportional  to  Dq/V^  as  predicted  by  theory. 

5.6  Plunging  (d)  Tests 

Test  with  pure  sinusoidal  plunging  with  boundary -layer  con¬ 
trol  as  required  to  prevent  separation.  Cover  a  sufficient  range 
of  and  >  to  establish  the  functional  relationship  of  each 

variable.  In  addition  to  the  complete  configuration,  test  the 
duct  alone  as  in  Section  5.4. 

The  present  theory  (for  low  frequency,  small  angle  of  attack, 
and  low  amplitude  oscillation)  predicts  that  the  duct  normal  force 
due  to  d  is  independent  of  k,  a^,  and  frequency;  that  is. 
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This  theory  should  be  compared  with  data  of  low  frequency,  but 
higher  frequencies  should  also  be  tested.  All  the  aerodynamic 
reactions  on  the  duct  and  propeller  should  be  measured.  Both  the 
steady  values  and  the  time  variations  should  be  measured  for  com¬ 
parison  with  previous  tests  (Sections  5.3  and  5.4)  and  with  theory. 

5.7  Predicted  Coefficients  and  Derivatives 

The  static  thrust  coefficients  and  dynamic  stability  deriva¬ 
tives  for  the  proposed  ducted  propeller  configuration  shown  In 
Figure  9  are  given  below.  They  are  obtained  from  Sections  3.1  and 

3.4  (with  c  ■■  0.2D)  as  functions  of  the  disk  loading  C_  and 

^P(D) 

the  flight  conditions  (a^,  q,  &) .  In  the  equations  given  below, 

^  ^  ®o’  **D(P)  Is  a  positive  pitching 

moment  about  the  duct  midchord  diameter,  and  Is  about  a 

propeller  diameter. 

Hovering  Plight  (V^  ■  0) 

For  hovering  flight,  the  predicted  thrust  ratio  Is,  from 
Equation  (17) , 

'^D(P).  2  a 

■= - S  ^  (In  80  -  2)®  -  0.724  (216) 

^P(D)^ 

Axial  Flight  (a^  ■  0) 

For  axial  flight,  the  predicted  thrust  ratio  Is,  from  Equations 


(11)  ,  (168) ,  and  (169)  , 
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^PlP)  .  0.724  (217) 

’'P(D)  1  +  2(Vo/Y> 


Steady  Flight  at  Angle  of  Attack 

For  this  flight  condition,  the  static  coefficients  of  the  duct 
in  the  presence  of  an  actuator  disk  are  from  Equations  (9) ,  (11) , 
and  (13)  . 


sin  a 


91  cos  a  +  0.845 


t) 


(218) 


z 

-  0.724  ^ —  +  1.46  sin^  a  (219) 

V  V  D(P)  V  ^  ° 

o 


(c\  “  1.45  ^  sin  a  (220) 

V  D(P)  ^o  ° 


The  only  static  derivative  of  an  isolated  dual  propeller  at 
small  angle  of  attack,  with  constant  blade  pitch  and  chord,  and 
with  'K  ■  V^/oD  ■  0.1,  and  o  ■  4Bb/3TrD  “0.1  is  from  Section  3.2.1, 
a  normal  force  derivative  given  by 


“  0.05,  0.10,  0.16  for  p 

V  “/  P 


since 


10,  20,  30°,  respectively, 

(221) 

(222) 


For  the  Hiller  propellers,  we  have  a 0.126  (from  Reference  13). 

Due  to  duct  interference  there  is  a  pitchihg  moment  on  the 
propeller  which  is  estimated  (from  Eg.  (65))  to  be  given  by 
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(^”*a)D(P) 


«#,  0.451,  0.015 


(223) 


with  the  propeller  located  at  the  entry,  central,  and  exit  planes 
of  the  duct,  respectively. 

Pitching  Flight 

For  small  rates  of  pitch,  the  duct  q  derivatives  in  the 
presence  of  an  actuator  disk,  are  from  Equations  (18),  (20),  and 
(22), 

/'C„  A  -  0.675  +  0.126  cos  a„  (224) 

®qjD(P)  ’'o  ° 


C'*^q)D(P) 

-0.0460  sin  a 

o 

(225) 

-  0.0376  ^ 

(226) 

V  "*qyD(P) 

The  q  derivatives  for  an  isolated  propeller  at  »  0 
are,  for  the  previously  assumed  propeller  parameters:  a  side 
force  derivative  given  by 


-  (%)p  ■  ° 


and  for  a  single  propeller  a  pitching  moment  derivative  given, 
respectively,  by 


^^m  )  “  ^ 


Note  that  ^  for  a  dual  propeller  is  not  given  in  Reference  5 

nor  is  it  calculated  herein. 
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Due  to  q  Interference  from  the  duct,  there  is  a  pitching 
moment  on  the  propeller  which  is  estimated  (from  Eq.  71)  to  be 
given  by 


(S) 


+  00,  -12.0,  -10.2 


(229) 


with  the  propeller  located  at  the  entry,  central,  and  exit  planes 
of  the  duct  respectively. 

Plunging  Flight 

For  small  angles  of  attack,  the  estimated  normal  force  deriva¬ 
tive  of  the  duct  due  to  plunging  (Eq.  (30))  is 


0.126 


(230) 


There  is  no  thrust  force  or  moment  due  to  d;  i.e.. 


The  estimated  reaction  on  the  propeller  due  to  plunging,  for 
constant  blade  pitch  and  chord  and  Oq  «  is  only  a  normal 
force  with  (see  Eq.  41) 


sin  ^ 


2 


(231) 
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6 .  CONCLUSIONS 

A  theoretical  analysis  has  been  presented  for  predicting  the 
forces,  moments,  and  stability  derivatives  of  an  isolated  ducted 
propeller.  The  analysis  is  an  extension  of  the  steady-flow  anal¬ 
ysis  of  Reference  1  in  which  the  duct  is  represented  as  a  straight, 
short,  thin  ring  surrounding  a  uniformly  loaded  actuator  disk 
which  produces  a  slipstream  velocity  much  greater  than  the  lateral 
component  of  the  flight  velocity  (V^  »  sin  a^) .  This  repre¬ 
sentation,  of  course,  requires  that  the  minimum  disk  loading  (for 
which  the  theory  is  valid)  increase  with  angle  of  attack.  In  the 
present  work,  the  approximations  of  Reference  1  have  been  investi¬ 
gated  in  detail  and  the  effect  of  moving  the  propeller  in  the  duct 
has  also  been  taken  into  account.  In  addition,  the  interference 
effect  of  the  duct  on  a  real  propeller  with  finite  blades  has  been 
incorporated,  and  the  cases  of  pitching  and  plunging  motions  have 
been  treated.  An  experimental  program  specifically  designed  to 
check  the  theory  has  been  outlined. 

6.1  Theoretical  Results  for  Isolated  Ducted  Propeller 

The  results  quoted  here  are  applicable  to  all  angles  of  attack 
to  the  extent  that  the  propeller  can  be  represented  by  an  actuator 
disk.  Any  results  pertaining  to  the  forces  on  a  propeller  with 
finite  number  of  blades,  however,  are  restricted  to  small  angles 
of  attack.  With  this  restriction,  the  following  conclusions  can  be 
drawn  from  the  theoretical  analysis  of  an  isolated,  moderately 
loaded,  ducted  propeller  of  solidity  0.1: 
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(1)  Except  for  thrust  force  and  pitching  moment  due  to  pitch 
(about  the  duct  center) ,  the  forces  and  moments  on  the  duct  are 
considerably  larger  than  those  on  the  propeller. 

(2)  When  it  is  located  within  the  duct  the  propeller  produces  a 

positive  damping^  f-C^  )  which  is  greater  than  the  negative  damping  due 

q 

to  the  duct.  This  fact  is  significant  for  STOL  airplanes  with  ducted 
propellers  mounted  on  a  lateral  axis  through  the  center  of  gravity. 

(3)  The  forces  and  moments  acting  on  the  duct  are  unaffected 
by  the  location  of  the  propeller  inside  the  duct. 

(4)  The  duct  interference  on  the  propeller  produces  a  large 
positive  d2unping,  unless  the  propeller  is  placed  near  the  duct  lead¬ 
ing  edge  where  it  produces  large  negative  damping. 

(5)  The  normal  force,  provided  by  a  ducted  propeller,  is  from 
20  to  60  percent  greater  than  for  a  flat  plate  of  the  same  projected 
area,  even  for  zero  disk  loading. 

(6)  All  of  the  static  and  pitching  derivatives  for  the  duct  in 
the  presence  of  the  propeller  are  found  to  increase  with  propeller 
disk  loading.  Some  derivatives  increase  nearly  proportionately  to 
disk  loading,  and  others  are  nearly  independent  of  it. 

6.2  Theoretical  Results  for  Torpedo-Like  Configurations 

From  the  analysis  of  a  torpedo-like  configuration  (in  axial 
flow)  consisting  of  a  body  of  revolution  with  rear-mounted  ducted 
propeller,  the  following  conclusions  are  drawn: 

(1)  The  dominant  contribution  of  the  ducted  propeller  to  the 
dancing  is  produced  by  the  normal  force  on  the  duct  due  to  its 
induced  angle  of  attack  arising  from  the  q  motion.  This  damping 
is  increased  at  higher  disk  loadings . 
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(2)  Both  C  and  C  due  to  the  hull  are  large  and  are 

m_  m. 

g  (S 

associated  with  the  hull  boundary  layer,  which  has  been  assumed  to 

oscillate  as  a  solid  body  with  the  hull.  Therefore,  a  more  refined 

analysis  seems  warranted,  even  though  and  due  to  the 

g  ct 

hull  tend  to  cancel  one  another. 

(3)  The  thrust  interference  between  the  ducted  propeller  and 
the  hull  appears  to  be  large,  depending  upon  the  geometry  and  the 
propeller  disk  loading.  Therefore,  a  more  refined  analysis  of  this 
effect  appears  to  be  warranted. 

(4)  Because  of  the  radial  inflow  produced  at  the  base  of  a 
symmetrical  hull  with  attached  flow,  such  a  hull  is  expected  to 
affect  the  stability  derivatives  of  a  rear-mounted  ducted  propeller 
in  the  same  direction  as  would  an  increase  in  propeller  disk  load¬ 
ing. 


6.3  Comparisons  with  Experiment 

The  comparisons  attempted  herein  between  the  present  theory 
and  available  experimental  data  have  indicated  the  following  con¬ 
clusions: 

(1)  Only  one  reference  (Ref.  14)  provided  data  free  of  flow 
separation  for  a  check  on  any  of  the  calculated  stability  deriva¬ 
tives.  This  check  showed  that,  for  small  angles  of  attack,  the 
present  theory  gave  excellent  predictions  of  the  static  stability 
derivative  .  The  dominant  factor  in  the  dynamic  stability  of 

Ct 

the  torpedo-like  configuration  at  zero  angle  of  attack  was  found 
to  be  the  lift  curve  slope  of  the  Isolated  ducted  propeller. 
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(2)  Only  one  reference  (Ref.  21)  provided  data  for  a  short 
duct  free  of  flow  separation  as  a  check  on  the  thrust  ratio  in 
axial  flow.  For  this  case,  the  theory  over-predicted  ’^j3(p)/^p (D) 
by  about  20  percent  at  the  lowest  advance  ratio.  The  error  was 
larger  at  higher  advance  ratios  because  of  duct  drag  not  included 
in  the  theory. 

(3)  In  one  case  of  a  rather  long  unseparated  duct  in  hovering 

flight  (Ref.  20),  the  measured  ratio  Tj^  was  far  above 

that  predicted  by  either  the  present  (small  c/D)  theory  or  the 
asymptotic  simple  momentum  theory  for  infinitely  long  ducts.  There 
thus  appears  to  be  a  real  need  for  a  theory  valid  for  ducted  pro¬ 
pellers  of  intermediate  chord/diameter  ratios. 

(4)  There  appear  to  be  no  e^qperimental  data  for  pitching 

moments  on  unseparated  ducts  during  steady  flow  at  angle  of  attack 

or  for  either  normal  force  or  pitching  moments  on  unseparated  ducts 

at  high  angles  of  attack.  However,  the  static  stability  derivative 

Cjj  predicted  by  the  present  theory  at  angles  of  attack  of  50°  and 
g 

80”  agreed  within  about  15  to  45  percent  with  the  data  of  Reference  13 
for  which  the  flow  was  probably  separated  at  the  duct  leading  edge. 

(5)  There  are  at  present  no  experimental  data  for  the  dynamic 
stability  derivatives  of  isolated  ducted  propellers. 

(6)  The  data  reduction  technique  presented  herein  succeeded 
in  bringing  theory  and  experiment  into  closer  agreement  than  would 
be  indicated  by  comparisons  made  at  the  same  advance  ratio  and/or 
blade  setting. 
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7.  RECOMMENDATIONS  FOR  FUTURE  WORK 

On  the  basis  of  the  foregoing  conclusions  and  the  investiga¬ 
tions  undertaken  In  the  present  report >  the  following  recommenda¬ 
tions  are  made  regarding  future  theoretical  and  experimental  work 
on  ducted  propellers : 

(1)  There  Is  an  urgent  need  for  systematic  experimental  data 
on  ducted  propellers  having  relatively  straight,  thin  ducts  of 
small  chord/diameter  ratio  without  leading-edge  separation  at  all 
flight  conditions  to  check  the  theory.  A  program  designed  for 
this  purpose  has  been  suggested  herein. 

(2)  A  theory  should  be  developed  for  treating  ducted  propel¬ 
lers  of  larger  chord/  diameter  ratios  than  are  treated  in  the 
present  theory.  In  addition,  the  effects  of  duct  camber  and  thick¬ 
ness  should  be  investigated. 

(3)  The  present  theory  should  be  extended  to  cover  nonuniform 
disk  loadings  and  wake  swirl. 

(4)  The  effects  of  interference  between  a  ducted  propeller 
and  a  wing,  ground  plane,  or  adjacent  ducted  propellers  should  be 
studied  theoretically,  especially  for  VTOL  configurations. 

(5)  The  theory  should  be  extended  to  study  the  speed  deriva¬ 
tives  (V)  of  a  ducted  propeller  at  angle  of  attack,  which  may  be 
significant  in  the  transitional  flight  regime  for  a  VTOL  aircraft. 

(6)  A  more  refined  analysis  should  be  carried  out  to  study 

the  mutual  interference  between  the  hull  and  rear-mounted  ducted  pro¬ 
peller  for  a  specific  modern  torpedo  and  for  a  high-speed  submarine. 
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Thls  analysis  should  Include  a  further  investigation  of  the  aver¬ 
aging  process  used  herein  because  of  its  possible  significant  effect 

on  the  duct  derivative  Cjj  Which  evidently  dominates  the  damping 

OL 

in  pitch  of  the  configuration. 

(7)  Engineering  calculations  should  be  made  to  investigate 
the  relative  Importance  of  the  ducted  propeller  stability  deriva¬ 
tives  on  the  overall  stability  of  vehicles  in  the  following  cate¬ 
gories  : 

(a)  VTOL  aircraft  employing  tilting  ducts. 

(b)  Vehicles  employing  multiple  ducted  propellers  in  which 

mutual  interference  may  be  important  but  ground 
effect  is  not. 


(c)  Ground-effect  machines. 
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APPENDIX  A 


DUCT  SURFACE  VELOCITY  DISTRIBUTION  W  FOR 
PURE-PITCH  PLIGHT  CONDITION 


As  explained  in  Section  2,  the  complete  velocity  distribution 
W  for  pure  pitch  is  composed  of  three  parts  (see  Fig.  2)  :  w.  +  w. 


due  to  hovering,  due  to  free  stream  at  angle  of 

attack  a  ,  and  w  +  w  due  to  pitching  rate  q  about  the 
duct  midchord  diameter.  Here  we  shall  determine  successively  these 
three  parts  of  W.  The  derivation  of  the  first  two  parts  will  be 
compared  with  Reference  1,  from  which  it  is  taken,  in  order  to 
point  out  certain  detailed  errors.  However,  the  final  result  of 
Appendix  B  (after  integration  of  and  cancellation  of  higher 

order  terms  in  c/D)  gives  net  duct  forces  and  moments  which  are 
identical  with  the  results  given  in  Reference  1  for  both  hovering 
and  steady  flight  at  angle  of  attack. 


A.l  SURFACE  VELOCITY  DUE  TO  HOVERING,  Wj^ 

For  no  motion  of  the  duct  (hovering  flight),  the  duct  surface 
velocity  distribution  w^^  is  induced  by  the  vorticity  shed  by 
the  actuator  disk  (y  in  Fig.  2)  and  by  the  bound  vorticity  7^^ 
which  is  generated  on  the  duct  to  satisfy  the  boundary  and  Kutta 
conditions  (considering  7^  to  lie  in  a  flat  surface) .  Finally, 
as  described  in  Section  2,  w^^  is  self-induced  by  7^  due  to 
the  fact  that  this  vorticity  actually  lies  on  a  ring  rather  than 
on  a  flat  surface. 

The  complete  surface  velocity  distribution  on  the  duct  for 
hovering  flight,  t  be  obtained  from  the  results  of 

this  section  by  adding  Equations  (A.  15)  and  (A. 47)  for  w^^  and 
Equation  (A. 53)  for  .  This  gives 
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On  the  inside  and  outside  of  the  duct  trailing  edge, 
is  3/4  and  -(1/4),  respectively.  For  obtain  the 

following  values 


c/D 

Wj^,/y 

0 

0 

0.1 

0.136 

.2 

.189 

.4 

.240 

Thus  is  somewhat  smaller  than  w^^  at  the  trailing  edge. 

It  is  assumed  that  the  disk  loading  C_  or  pressure 

^P(D) 

rise  across  the  disk  Ap,  is  constant  over  the  disk,  which  is 
located  at  the  exit  plane  of  a  thin  duct  of  constant  diameter 
D.  The  flow  field  generated  by  the  actuator  disk  is  induced  by 
the  shed  vorticity  y  which  surrounds  its  slipstream.  The  y 
vorticity  distribution  is  approximated  by  a  sheath  of  vorticity 
uniformly  distributed  over  a  semi-infinite  cylinder  as  shown  in 
Figure  2  and  Sketch  A.l. 
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Sketch  A.I.-  Ducted  actuator-disk  nomenclature. 


A.  1.1  Velocity  Induced  by  y  Shed  From  Actuator  Disk 

The  velocity  components  induced  directly  by  y  at  the  duct 
surface,  are  (p.  319,  Ref.  6) : 

[f  Vl  -  k^  K(k)j  (A.2) 

where  K(k)  and  E(k)  are  the  complete  elliptic  integrals  of  first 
and  second  kind,  respectively,  as  defined  in  the  List  of  Symbols. 

The  series  expansions  of  these  functions  are  (p.  307,  Ref.  6) : 
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1 


-  - 


IfZDi 


(A. 5) 


(A. 6) 


(A. 7) 


and  for  Equations  (A. 2)  and  (A. 3)  we  have  ^  as  shown  in 

Sketch  A.l.  Now,  if  j  «  |d  |,  then  k^  *  1,  and  the  leading 
terms  in  the  series  expansion  of  E  and  K  are  simply 


K(k) 


in 


4D 

e 


(A. 8) 


E(k)  “  1 


(A. 9) 


The  duct  lies  on  the  surface  r  ■  D/2,  -c  <  <  0.  Since 

we  shall  assume  here  that  c  «  D,  the  velocity  components  induced 
by  the  actuator  disk  on  the  duct  surface  are  given  approximately 
by  the  leading  terms  in  the  series  expansion  of  Equations  (A. 2) 
and  (A. 3)  which  are  (for  ^ 


u 


y 


a 


(A. 10) 


I 
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y 


a 


2^t  ' 

(2  -  In  8)  +  In 

D 

(A. 11) 


These  equations  agree  with  Equations  (A. 2)  and  (A. 3)  within 
1  percent  and  10  percent,  respectively,  when 


(A.  12) 


A.1'2  Complex  Potential  F  Required  for  Boundary  Condition 

In  order  to  make  the  duct  a  streamline,  a  distribution  of 

vorticity  must  be  added  along  the  duct  chord  which  induces  a 

radially  inward  velocity  component  v^  equal  to  the  radial 

component  v  (Sketch  A.l).  To  do  this,  we  shall  find  a  complex 
"^s 

potential  function  F(z)  which  gives  v^  ■  v  by  treating  a 
chordwise  strip  of  the  duct  as  a  two-dimensional  airfoil.^  It  is 
convenient  for  this  purpose  to  introduce  the  Joukowski  transfor¬ 
mation  (p.  159,  Ref.  12,) thus 

z  *  4  (c  +  ■  Xg  +  iy  (A. 13) 

and  with  the  actuator  disk  at  the  trailing  edge  of  the  duct,  we 
have  from  Sketch  A.l 


Xg  -  ^t  +  f  (A. 14) 

The  chordwise  strip  (-  c/2  <  x  <  +  c/2)  maps  on  to  the  unit 

s 

16 

circle  C  ■  e  such  that  on  the  duct  surface  z  ■  x  ■  c/2  cos  9. 
a  a 

Thus,  we  can  write  Equations  (A. 10)  and  (A. 11)  in  terms  of  6  as 
follows : 


^This  approach  is  justified  with  the  restriction  that  c  «  D. 
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y 

U  “  T" 
Ts  27r 


r  a:  _  ^ 

1  2  2D 


(1  -  cos  0) 


[in  ^  -  in 


(1  -  cos  6) 


} 


Tg  27r 


^2  -  in 


8D 


+  ln(l  -  cos  0) 


(A. 15) 

(A. 16) 


Let  the  required  complex  potential  function  in  the  ^  plane 


F  •  *  <l>  *  +  * 


27r 


I 


00 

n=c»  ^ 


(A. 17) 


Then  the  complex  velocity  in  the  ^  plane  is  given  by  (p.  150, 
Ref.  12) 

dP ' 

“p'  -  ”f'  ■  *  aT 


Hi 

27rC 


la  n 
n 

•  (n+i) 


iOy 


On  the  circle  (C  *  e  ),  we  have 
s 


ir 

"f  '  ■  '  “  iTT 

s  s 


—  e”^®  -  ^®n”® 

n*6 


- (n+i) i0 


(A. 18) 


(A. 19) 


The  complex  velocity  at  a  corresponding  point  in  the  physical 
plane  is 


“p  -  ^''p  ■  ^  al*  <"p‘  -  i''p'> 


-  t) 


(A. 20) 


iO, 


so  that  on  the  flat  plate  (C  *  e  )  we  have,  in  terms  of  0, 

s 


s  s 


(u,;  - 


2e 


i0 


ci  sin  0 


(A. 21) 
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Substituting  Equation  (A. 19)  into  (A. 21)  gives 


u  —  iv  »  . .  ■"  ^ 

F  TTC  sin  e 
s  s 


-i„e 
e 


c  sin  6 


(A. 22) 


Therefore,  Vp  is  finally  found  to  be 


n»o 


+v-  sin  0  •  Im 
F 


(cos  n0  -  i  sin  n6) 


(A. 23) 


Now  let  us  express  v  in  the  same  form  as  Vp  .  The 
following  identity  can  be  found  on  page  202  of  Reference  22) 


ln(l  -  cos  0)«-ln2-27  ~  cos  nd 


(A. 24) 


By  use  of  this  and  the  identity 

2  sin  9  cos  n0  =  sin(l  +  n)0  +  sin(l  -  n)0 
Equation  (A. 16)  can  be  written  in  the  form 


(A. 25) 


sin  0  • 

7s  2ir 


{0  - 

00 

TU- 


sin  6 


sin(n  +  1)0  -  sin(n  -  1)0 


1} 


(A. 26) 


Now,  setting  v  (Eq.  (A. 26))  equal  to  Vp  (Eq.  (A. 23))  in  order 
's  's 

to  satisfy  the  boundary  condition  of  zero  normal  velocity  through 
the  duct  surface,  we  find  that  (a^)  is  real,  so  that  Equation  (A. 23) 


can  be  written  as 


Vp  sin  0 
‘a 


^  2a  n 

-  y  -  .i 


sin  n0 


(A. 27) 
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Substitution  of  v.  for  v_  (Eq.  (A. 26)  into  Eq.  (A. 27))  then 
gives 


00 

^  -  m  -^^^sin  e  + 


4Trna, 


yc 


sin  n@ 


00 

Esin(n  +  1)8  ^  sin(n  -  1)0  _  q 
n  n  ” 

n“i 


(A. 28) 

By  reindexing,  this  equation  can  be  written  in  the  following  form; 


sin  nd 


sin  nd  f. 
’  2^  n  -  1 
n“a 


(A. 29) 


Consequently,  the  coefficients  a^^  are  found  to  be 


■»  0 
o 


(A. 30) 


a  - 

n 


Yc 


for  n  >  1 


27rn(n  -  1) 

Thus,  Equation  (A. 17)  for  the  complex  potential  can  be  written  as 


ir 


F'  «  (f  +  if  •>  I 


n  ^  +  1  )  - ; -  +  ~T~ 

^  n(n^  -  1)  ^ 


(A. 31) 
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Now,  on  the  circle  (C  *  e  ) ,  we  have 

8 


p  *  »  $  *  +  * 

S  S  ^8 


£hi 

27r 


*  T  ‘*n 

n»6 


(cos  n©  -  i  sin  n©)  (A. 32) 


Thus,  substitution  of  Equation  (A. 30)  gives 


*  '  ■  -  — fl— 
*s  2Tr 


oa 

+  Is.  V'  sin  n©  . 
2^  ^  n(n^  -  1) 


sin  © 


(A. 33) 


In  order  to  evaluate  the  series  the  following  technique  is  used 
(from  Ref.  1):  We  use  the  identity  from  Reference  22 ,  page  188 


£ 

n“i 


sin  n©  TT  -  © 


0  <  ©  <  2ir 


(A. 34) 


We  then  form  the  product 
sin  n©  _ 


(cos  ©  -  1) 


n 


1  sin(n-H)0  +  sin(n 


-1)0  -  2  sin  n© 

(A. 35) 


Now,  by  reindexing,  we  can  write  Equation  (A. 35)  as 
a» 

-SAa,  -  «  I— 7— Q  (pQg  e  _  2)  +  ^  sin  ©  ,  0  <  0  ^  27r 

^  n(n*  -  1)  2  4 


(A. 36) 


Thus,  substitution  into  Equation  (A. 33)  yields 

r^.© 


1  V  r*  * 

^  (cos  0-l)+§-sin6  +  a  sin© 

s  2ir  2Tr  L  2  ’  J  1 


(A. 37) 
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Finally,  substitution  of  Equation  (A. 30)  and  differentiation  with 
respect  to  6  gives  the  surface  velocity  distribution  on  the  unit 
circle  due  to  F.  That  is, 


^h  7c  r  1  -  cos  9  e  -  IT  .  ^  .  3 

-  27  ^  27  L - 2 -  ^  ~2  0  +  4  cos  e 

+  In  -  J)  cos  (A. 38) 


A. 1.3  Fulfillment  of  the  Kutta  Condition 

The  value  of  the  total  circulation  around  the  duct  in  hovering 
Tjj  will  now  be  determined  such  that  the  Kutta  condition  of  finite 
velocity  is  satisfied  at  the  trailing  edge  of  the  duct.  Mathemati¬ 
cally,  this  requires  that  the  trailing  edge  be  a  stagnation  point 
in  the  transformed  plane.  That  is,  w^'  ■  0  at  0*0.  From 
Equation  (A. 38),  with  0  ■  0,  we  have  at  the  duct  trailing  edge 


("f' 


I 


= 


f) 


Ih 

27r 


-  if 


This  is  only  the  part  of  w^'  which  is  due  to  the  potential 
function  F.  However,  w  ',  the  part  due  to  y  is  zero,  as  shown 
below.  From  Equations  (A. 15)  and  (A. 16)  with  0  ■  0,  we  find 
that  the  velocity  components  induced  by  7  at  the  trailing  edge 
of  the  plate  are 


u 


1 

4 


JL 

2ir 


[2  -  in 


8D 


+  ln(l  -  cos  0) 


(A. 40) 
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The  complex  velocity  in  the  circle  (C)  plane  is  related  to 
that  in  the  physical  plane  by 


u’-iv‘-fu  -iv'i^ 

y  y  \  y  y/^fi 

(A. 41) 

Evaluation  at  0  ■ 

0  gives 

u  '  -  iv  '  •  ( 
7t  Yt  V 

"^u  -  iv  '^^fsin®  S  +  i  sin  6  cos  6) 

^  Yt  Yt^  2  V  J 

L-o 

(A. 42) 

Hence,  substitution  into  Equations  (A. 40)  gives 

ft 

1 

^  sin  d  cos  0  ln(l  -  cos  0) 

®  '0-0 

(A. 43) 

V  ' 

'  -  ^  sin^  0  ln(l  -  cos  0)1 

10-0 

(A. 44) 

Evaluation  of  Equations  (A. 43)  and  (A. 44)  shows  that  they  are 

each  equal  to  zero;  therefore,  w  *  is  zero. 

't 

Thus,  if  we  set  w„  '  equal  to  zero  in  Equations  (A. 39) 

t 

to  satisfy  the  Kutta  condition,  we  find  that 

^  0"  ^  -  0 

Finally,  substitution  of  Equation  (A. 45)  into  Equation  (A. 38) 
gives  the  surface  velocity  on  the  circle 

Wp  '  ■  ^  8)  +  (6  -  ir)sin  ©j  (A. 46) 

8 

Hence,  we  find  from  Equation  (A. 41)  that  the  surface  velocity 
distribution  on  the  plate  due  to  y 
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W  '  _ _ 2 - 

F  dx„  c  sin  6 
8  8 

■  i  [G  -  2 


(A. 47) 


We  have  now  found  the  hovering  velocity  distribution  (wj^  ■  u  +  Wp  ) 

^s  s 

induced  by  y  and  This  result  (Eq.  (A. 15)  plus  (A. 47))  is 

identical  with  Equation  (12)  of  Reference  1.  Now  we  seek  the 


averaged  self-induced  component  >  as  discussed  in  Section  2. 

A.  1.4  Self-Induced  Velocity  Wj^ 

As  described  in  Section  2,  there  is  a  surface  velocity  compo¬ 


nent  Wj^  which  is  self-induced  by  the  bound  vorticity  due 

to  the  fact  that  it  lies  on  a  ring  rather  than  on  a  flat  surface. 


We  shall  find  >  i^s  average  value,  by  summing  the  elementary 
velocity  components  induced  by  the  vorticity  (y^^  each 

elementary  ring  d^^  of  the  duct  (Sketch  A. 2). 


Sketch  A. 2.-  Velocity  components  induced  by  a 
vortex  ring  of  constant  strength. 
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Each  elementary  ring  consists  of  a  vortex  ring  of  constant  strength 
(with  respect  to  0)  equal  to  7^^  where  is  the  axial 

distance  from  the  ring  to  a  fixed  point  on  the  duct  surface.  Thus 
-c  <  <  c*  For  c/D  «  1,  we  can  write  (see  pp.  306  and  307  of 

Ref.  6) 


Thus,  dVj^  ,  the  radial  velocity  component,  is  equal  to  that 
“s 

for  a  straight  line  filament  of  infinite  length  (D  =  <»)  if 

«  D.  Consequently,  v,  and  the  boundary  condition  have  been 

”s 

accounted  for  correctly  in  the  foregoing  (flat  plate)  analysis. 

However,  we  must  account  for  Uj^  ,  the  axial  component  of  velocity 

s 

which  is  self-induced  by  the  bound  vorticity  on  the  duct,  since 
Uy^  is  zero  for  a  straight  line  filament. 

Following  the  method  used  in  Reference  1,  both  the  bound 
vorticity  distribution  7^^  and  the  self-induced  velocity  w^ 
will  be  averaged  over  the  chord  as  follows-  Corresponding  to  the 
surface  velocity  w^^,  we  obtain  the  bound  vorticity  distribution 

7j^  -  ~  (A. 49) 

Thus,  using  Equations  (A.  15)  and  (A. 47)  for  Wj^,  we  have 


'Yh 


tan  J  +  T 


-  e 


(A. 50) 
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The  average  value  of  over  the  chord  is  therefore 


h  c 


+  *  TT 

J  '^h  "  2  J  '^h  ®  "  2  C"  ^ 

-  f  ,  o  > 


For  this  constant  vorticity.  Equation  (A. 48)  gives 


Ijl. 

%  “  2TrD 


(in  4D  -  1  -  in  |) 


■5? 
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-  ^  +X 
2  *S 


Ih_ 

2TrD 


{<ln  -  1)=  -  (f  -  \)  [l"  (f  -  “s)  ■ 
-  (f  *  *.)  [^"(f  *  -n)  -  ']} 


(A. 52) 


Averaging  this  expression  over  the  chord  gives  the  averaged  self- 
induced  velocity  Wj^  as 


w, 


hi  c 


"h  ‘^^s  -  27rD 


^(in^.l) 


c 

2 


Yc  A  16D  4D  ^  1  \ 

■  4T5  -  vO"  T  ^  i; 


(A. 53) 


This  equation  is  identical  with  Equation  (13)  of  Reference  1, 
except  for  the  last  term  (t  1/2)  which  is  incorrectly  given  as 
(-  1/4)  in  Reference  1.  The  accuracy  of  the  averaging  process  is 
examined  in  Appendix  D. 
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A. 2  SXJRPACE  VELOCIIY  DISTRIBUTION  FROM  AT 

o  o 

Here  we  shall  find  the  second  part  of  W  shown  in  Figure  2, 

neunely  free  stream  at  angle 

of  attack  a^.  The  first  component  w^  is  induced  by  the  bound 

vorticity  y  which  is  generated  to  satisfy  the  boundeury  condi- 

tion  and  the  Kutta  condition,  and  by  the  wake  vortex  fileunents 

dr  which  are  shed  from  the  duct  at  angle  of  attack.  Finally, 

“w 

as  described  in  Section  2,  w_  is  self-induced  by  which  lies 

’  '  a 

on  a  circle. 

A. 2.1  Determination  of  w^ 

The  elementary  velocity  induced  by  a  semi-infinite  trailing 

vortex  filament  in  the  duct  wake  (dr  in  Fig.  2  and  Sketch  A. 3) 

“w 

is  given  by  the  Biot-Savart  law  (p.  304,  Ref.  6) .  Consider  the 
radial,  axial,  and  tangential  components  of  the  induced  velocity. 

The  variation  of  the  radial  component  along  the  chord  is  negligible 
because  the  trailing  vortex  fileunents  extend  to  +  « .  Consequently, 
we  can  set  c  >  0  as  in  Sketch  A. 3.  Also,  there  is  no  axial  com¬ 
ponent  of  induced  velocity,  and  the  effect  of  the  tangential  com¬ 
ponent  is  considered  to  be  negligible  here  (as  in  Ref.  1) . 

The  radial  conponent  of  velocity  induced  at  (p  by  the  vortex 

filament  dr  at  0'  is  given  by  (p.  304,  Ref.  6) 

“w 

OL 

- - /  I"  T' ^-7-^ 

^s  4TrD  tan  ^  j 


(A. 54) 


Sketch  A. 3.-  Trailing  vortex  filament. 


The  strength  of  the  trailing  vortex  filament  dr  is  equal  to 

“w 

the  decrease  in  r^,  the  flight  induced  bound  vorticity  on  the 
duct,  as  one  moves  around  the  ring.  Thus,  we  can  write 


(A. 55) 


Now,  for  a  flat  plate  in  uniform  flow  at  angle  of  attack,  the  Kutta 
condition  gives  the  circulation  as 


r-  “  TTCV^  sin 
£  o  o 


(A. 56) 
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where  sin  is  the  normal  velocity  component.  The  normal 
velocity  at  each  chordwise  strip  of  the  duct  due  to  both  the  free 
stream  and  the  waOce- induced  velocity  (Sketch  A. 3)  is  given  by 

sin  cos  A  -  v. .  (A. 57) 

o  o  w 

8 

Thus,  the  circulation  about  the  duct  chord  is,  from  Equations  (A. 55) 
and  (A. 56), 

“o  ■  fe  /  (di^)  ^-=-^  d0  (A. 58) 

o 

The  method  used  in  Reference  1  to  solve  this  integral  equation  for 

r~(0)  is  to  assume  initially  that 
u 

cos  <f>  (A. 59) 

where  is  a  constant,  such  that 

dr„ 

^  sin  0*  (A. 60) 


Substitution  of  this  value  and  the  trigonometric  identity 


_4.  <P  -  4>'  1  t  cos  (0  -  0' ) 

2  sin(0  -  0') 


(A. 61) 


into  the  integral  in  Equation  (A. 58)  yields 
?Tr 


j  sin  0'  cot  ^  (j0'  -  ”27r  cos  0 


(A. 62) 


Substitution  of  this  value  into  Equation  (A. 58)  then  gives 


^k  “  Cyo  “o  “  2d) 


(A. 63) 
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which  satisfies  the  initial  assumption  that  does  not  depend 

on  0.  Solving  for  Fj^  gives  the  0  variation  of  duct- circula¬ 
tion  from  the  free  stream  and  wake  when  the  Kutta  condition  is 
satisfied  everywhere  at  the  trailing  edge.  Thus,  we  have  from 
Equations  (A. 59)  and  (A. 63), 


sin  a 


IS. 

2D 


cos  0 


(A.64) 


The  term  7rc/2D  in  the  denominator  is  due  to  the  wake.  The 
Surface  velocity  distribution  w^  induced  on  each  chordwise  strip 
of  the  duct  by  the  free  stream  and  wake  is  assumed  to  correspond 
with  that  for  a  flat  plate  at  incidence  so  that  (from  Eqs.  (A. 78) 
and  (A. 80)) 


& 

V  •  cos  a  +  — “  tan  -r 
a  o  o  TTC  2 


(A. 65) 


and,  with  Equation  (A.64) ,  we  find  that 

V  sin  a  g 

w^  -  cos  cos  0  tan  j  (A. 66) 

This  equation  is  identical  with  Equation  (17)  of  Reference  1. 

A. 2. 2  Determination  of  w^ 

Finally,  the  self-induced  velocity  W  due  to  the  circula- 

“i 

tion  Fq  about  the  duct  must  be  included.  To  do  this,  we  shall 
first  find  an  ei^ression  for  the  elementary  axial  velocity  com¬ 
ponent  induced  at  radius  r  -  D/2  by  an  elementary  vortex  ring 
(Sketch  A. 4)  having  diameter  D  and  circulation  dF^  ■  dFj^  cos  0', 
in  accordance  with  Equation  (A. 59). 
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Sketch  A. 4.-  Velocity  components  induced  by  a  vortex  ring 
with  variable  strength. 


Then,  by  using  the  results  on  pages  305  through  307  of  Reference  6, 
we  can  write 


-  -  2iD 


ZTT 

in 


[cos(0  -  0')  -  iJ  cos  <t>' 


4^. 


+  2—2  cos(0  "•  0*) 


3/2 


d0' 


(A. 67) 


Further,  by  letting 


= 


(A. 68) 


and 


_  <J)  -  (t)'  +  TT 

~  2 


(A. 69) 
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we  can  rewrite  Equation  (A. 67)  in  the  form 


du 


2 


a 


27rD 


s  _  k' 


ar^  2 


2 


2  cos  0  co8*y  -co8  <t>  cos^X  +  2  sin  0  sin  X  cos^y 
/ 

(l  -  k®  sin®  X  ) 

(A. 70) 

The  limits  of  integration  can  now  be  shifted  to  0  tt,  and,  since 
the  last  term  integrates  to  zero,  Equation  (A. 70)  becomes 

cos®X  -  sin®x[(l  -  k®sin®X)  +  (k®  -  1)] 


dX 


du 


a 

-2irD  ^  ■  k®cos  0 


(l  -  k®  sin®  X  ) 


3/2 


dX 


(A. 71) 


Now,  using  the  relations  given  on  page  306  of  Reference  6  for  the 
complete  elliptic  integrals,  we  find  that  the  above  equation  can 
be  written  as 


du 

a 


(A. 72) 


where  K(k)  and  E(k)  are  defined  in  the  list  of  symbols.  Finally, 
using  the  leading  terms  in  the  series  expansions  of  K  and  E 
(Eqs .  (A.  8)  and  (A. 9)  for  «  D,  k®  “  1,  we  obtain  the  desired 
expression 


du. 


-dr 


a 


a. 


2wD 


A  , 

4D 

(3  -  in 

(A. 73) 


This  result  is  four  times  smaller  than  the  corresponding  expression 
given  by  Equation  (19)  of  Reference  1.  In  a  similar  manner,  it  can 
be  shown  that  the  radial  velocity  component  dv 


a 


for  the  ring 
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filament  is  the  same  as  for  a  line  filament  for  «  D.  Thus, 
it  is  not  considered  here,  since  it  has  already  been  accounted 
for  in  Equation  (A.66) . 

Let  us  now  find  w  by  the  averaging  process  described  in 

“i 

Section  2.  The  average  vorticity  (or  circulation  per  unit  length) 
along  the  chord  length  is,  from  Equation  (A. 64), 


—  ^ 

/w  s  — m 


'a 


1  + 


IS. 

2D 


cos  </) 


SO  that  each  elementary  vortex  ring  (Sketch  A. 4)  has  strength 


(A. 75) 


At  X  ,  a  fixed  point  on  the  duct  surface,  the  velocity  induced 


by  the  bound  vorticity  y  distributed  over  the  chord  is  (from 
Eqs.  (A.  73)  and  (A. 75)): 
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2TrD 


/ 

-  S-  +  X 
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^  (f  ^ «.)  [‘“  (f  ='.)  - 1]} 


(A. 76) 
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Averaging  this  expression  over  the  chord  gives  the  average  induced 
velocity 


1 

c 


2 


dx_  ■ 
a  s 
8 


sin  a 

O _ O 

1  +  2D 

2D 


£sL  fl  4D\ 

”  2TrD  V2  "  c  J 

Cf  " 


(A. 77) 


This  result  is  four  times  smaller  than  the  corresponding 
Equation  (20)  of  Reference  1,  and  the  term  3/2  appears  here  within 
the  parentheses  rather  than  the  9/4  of  Reference  1.  The  accuracy 
of  the  averaging  process  used  both  here  and  in  Reference  1  is 
examined  in  Appendix  D. 


A. 3  SURFACE  VELOCITY  DUE  TO  PITCHING  q 

Here  we  shall  find  the  third  and  last  part  of  W  indicated 
in  Figure  2,  that  is,  w  +  w  which  is  due  to  the  pitching  motion 

q  q^ 

of  the  duct  q  about  its  midchord  diameter.  Due  to  this  motion, 
each  chordwise  strip  of  the  duct  translates  horizontally  and  rotates 
about  a  midchord  axis  (Sketch  A.5) . 
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Thus,  will  be  considered  in  three  parts: 

w  due  to  the  rotation  of  each  chordwise  strip  about  its 
midchord  axis 

w  due  to  the  horizontal  translation  of  each  chordwise  strip 

w  induced  by  the  wake  shed  from  the  duct  as  described  in 

"a  Section  2 . 


Finally,  w  is  self-induced  by  the  bound  vorticity  y  which 
‘’i  " 

corresponds  to  w^  and  which  lies  on  a  circle.  We  will  neglect 

the  effect  of  the  yawing  motion  of  the  chordwise  strips,  so  that 

for  each  strip  w^  is  due  only  to  the  pitching  component  of  its 

rotation.  This  is  consistent  with  Section  A. 2  where  only  the 

component  of  at  which  was  normal  to  each  strip  was 

considered  to  be  effective. 

The  final  results  show  that  the  predominant  effect  of  q  is 

caused  by  w  (the  horizontal  translation)  when  c  «  D.  Thus, 
‘*2 

the  net  force  and  moment  on  the  duct  due  to  q  are  nearly  the 
same  as  for  a  duct  which  distorts  with  each  chordwise  strip 
moving  parallel  to  itself  with  a  velocity  (qD/2)  cos  <p  (see 
Sketch  A. 5).  With  the  initial  assumption  that  c  «  D,  the  flow 
over  each  chordwise  strip  is  nearly  two  dimensional  so  that  we 
can  obtain  w_  and  w  from  flat-plate  formulas.  We  shall, 
therefore,  start  the  analysis  by  obtaining  these  formulas  which 
will  be  useful  later  when  we  consider  the  flow  over  the  chordwise 
strips  of  the  duct. 
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A. 3.1  Flat-Plate  Formulas 

Consider  the  flat  plate  in  Sketch  A. 6  which  is  in  a  pure-pitch 
flight  condition  with  not  necessarily  small, 


\ 

\ 


Sketch  A. 6.-  Flat  plate  in  pure  pitch. 


We  wish  to  find  the  surface  velocity  distribution  ''f  *  ''f  +  w- 

a  q 

and  the  bound  vorticity  distribution  flat 

plate  which  are  generated  by  the  two  parts  of  its  motion,  transla¬ 
tion  V  at  angle  a  ,  and  rotation  q,  about  its  midchord  axis, 
o  o  r 

For  the  first  part,  at  a^,  with  the  Kutta  condition  fulfilled, 

we  have  (Ref.  8,  p.  38) 


f 

w,  -  cos  a  +  sin  a  tan 

o  o  o  0  2 

a 


(A. 78) 
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f 

-  2Vq  sin  tan  (A. 79) 

Tf  -  TTcV^  sin  (A. 80) 

and  the  thrust  force  due  to  leading  edge  suction  is  given  by 

■  Trpc(V^  sin  (A. 81) 


To  find  the  second  part,  we  shall  use  the  well  known  "equi¬ 
valent  camber"  rule  as  follows;  The  bound  vorticity  y,  must 

q 

induce  a  velocity  normal  to  plate  equal  to  the  normal  velocity 
component  of  the  plate  (v^  •  “qx^  in  Sketch  A. 6)  in  order  to 

q 

satisfy  the  boundary  condition  of  no  flow  through  the  plate.  The 
Scune  distribution  of  normal  velocity  is  caused  by  a  circular-arc 
camber  (Sketch  A. 7)  for  which 


y  -  h  ^1  -  4 

^  c 


) 


(A. 82) 


Thus,  the  slope  of  the  cambered  airfoil  is 


dx. 


2 

C 


and  the  normal  velocity  is 


-  X  V  cos  a 

2  C  O  O 

c 


(A. 83) 


(A. 84) 


Thus,  by  equating  v  and  v-  ,  we  obtain  the  equivalent  camber 

q 

h  corresponding  to  the  pitch  rate  q^  of  the  plate.  That  is 


£l  .  _ _ 

c  8V^  cos 

o  o 


(A. 85) 
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Sketch  A. 7.-  Cambered  airfoil  nomenclature. 


Now,  since  circular-arc  camber  produces  a  surface  velocity  distri¬ 
bution  given  by 


w 


cos  a  sin  9 
o  c 


(A. 86) 


when  the  boundary  and  Kutta  conditions  are  satisfied  (see  Equa¬ 
tion  (10.7)  of  Reference  8),  we  find  for  the  equivalent  canOser  of 
Equation  (A. 85) 

w^  -  — j-  «  sin  (A. 87) 

Now,  with  Equations  (A. 78)  through  (A.  87)  for  a  flat  plate,  we  are 
ready  to  find  the  bound  vorticity  distribution  and  surface  velocity 
distribution  induced  on  the  duct  by  its  pitching  motion  about  its 


midchord  diameter. 
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A. 3. 2  Determination  of 

The  pitching  rate  of  each  chordwise  strip  caused  by  pitching 
of  the  duct  about  its  midchord  diameter  (Sketch  A. 5)  is  q  cos  0. 
Therefore,  using  this  value  in  Equation  (A. 87)  for  q^  gives  the 
first  component  of  w^;  namely, 


w. 


qc 

—  cos  0  sin  6 


(A.  88) 


so  that 


2 


A. 3. 3  Determination  of  w 

*^2 

The  translational  velocity  of  each  chordwise  strip  of  the 
duct  caused  by  its  pitching  about  the  midchord  diameter  (Sketch  A. 5) 
is  -q (D/2) cos  0.  This  chordwise  translation  produces  the  surface 
velocity  distribution 

w  ■  +  ^  cos  0  (A. 90) 

%  ^ 

which  is  clearly  equal  on  inner  and  outer  surface.  Hence, 

7q  -  Tq  -  0  (A. 91) 


A. 3. 4  Determination  of  w^ 

Since  varies  with  0,  there  must  be  vortex  filaments 

dr  vdiich  are  shed  parallel  to  the  duct  axis,  as  shown  in  Figure  2. 
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The  effect  of  the  duct  wake  which  is  generated  by  the  q  motion 
can  therefore  be  found  by  analogy  with  the  results  obtained  in 
Section  A. 2.1,  since  the  wake  due  to  angle  of  attack  shows  a 
similar  character.  The  net  duct  circulation  generated  by  the  q 
motion  is  r  «  r  since  r  is  zero;  and  r  (Eq.  (A. 89))  is 
analogous  to  the  first  term  on  the  right  side  of  Equation  (A. 63) 
(where  qc/4  now  replaces  sin  a^) • 

Thus,  we  obtain  the  total  circulation  including  the  wake 
effect  by  analogy  with  Equation  (A.  64) : 


r  +  r  .  £os^  / K  2) 

q  q_  1  +  K£.  ^  / 

^  2D 


(A. 92) 


The  circulation  due  to  the  wake  is  (by  subtraction  of  Eq.  (A. 89)) 

r„ 


«  IS. 
2D 


1  +  IS 
^  2D 


(A. 93) 


The  q  wake  induces  a  nearly  constant  normal  velocity  over 
each  chordwise  strip  of  the  duct,  since  each  vortex  filament  in 
the  wake  extends  to  +  <»  (Fig.  2) .  Thus,  the  chordwise  distribu¬ 
tions  of  surface  velocity  w  and  bound  vorticity  7  (i.e., 

"a  "3 

2w  )  which  compose  F  are  (from  Eqs.  (A. 78)  and  (A. 80)  with 

%  % 

“o  “  ) . 


ly 

2 


ire 


tan  I 


■n-qc^  cos  0 

eo  (1  +  f§) 


(A. 94) 
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A. 3. 5  Determination  of  w_ 

The  self  induced  velocity  due  to  the  bound  vorticity  in 

pitch  (y  “  Yq  Yq  )  can  be  found  by  analogy  with  the  results 
in  Section  A. 2. 2  for  y^.  If  we  replace  sin  by  qc/4,  as 

in  Section  A. 3. 4,  we  obtain  from  Equation  (A. 77) 


_  3si 

8D 


1  _ 


1  + 


In 


IS. 

2D 


cos  0 


(A. 95) 


APPENDIX  B 


INTEGRATION  OP  TO  OBTAIN  DUCT 

COEFFICIENTS  AND  DERIVATIVES 


Expressions  for  the  net  force  and  moment  on  the  duct  will  be 
developed  here  for  the  pure-pitch  flight  condition  (Fig.  2)  by 
integration  of  the  pressure  distribution  p_  over  its  surface. 

These  expressions  are  then  used  to  obtain  the  static  coefficients 
and  q  stability  derivatives  of  the  duct  in  the  presence  of  the 
propeller. 

The  pressure  distribution  P-(d|0)  is  related  to  the  surface 

s 

velocity  distribution  W  by  the  pressure  equation  in  moving  coor¬ 
dinates  (p.  87,  Ref. 12)  as  applied  to  the  pure-pitch  motion  (Fig.  2) 


P.  *  fo  ♦  f  (V  -  -  PsO 


where  V^  is  the  velocity  (in  a  fixed  frame  of  reference)  of  the 
point  in  question  on  the  duct  surface.  Since  the  first  two  terms 
on  the  right  side  are  constant  over  the  duct  surface,  they  do  not 
contribute  to  the  net  force  ot  moment  Further¬ 

more,  the  velocity  V  is  the  seune  on  the  inner  and  outer  surface 
and  hence  does  not  contribute  to  Np^pj  or  Thus,  the 

duct  force  and  moment  can  be  found  by  substituting  -(p/2)W^  for 

p  in  the  following  expressions: 

8 


**D(P) 


dx^  j  cos  0  d0 


aTT  aTT 


sin  9  cos  0  dO  d0 


(B.2) 


(«n). 


PgXg  dXg  -J  cos  <t>  <30 


sir  sir 


I 


p  cos  9  sin  6  cos  0  <36  <30  (B.3) 

s 


The  duct  force  and  moment  due  to  leading-edge  suction  T  are 
given  by 


a  n 

T  -  r  D 

J  2 


(T)  d0 


(B.4) 


(”t)d(p)  “  /  I 


(T)  COS  0  d0 


(B.5) 


where  T  will  be  evaluated  subse<2uently.  The  net  moment  on  the 


duct  in  the  presence  of  the  propeller  (Fig.  la)  is  then 


"d(p)  ■  («n)d  (P)  *  D(p) 


(B.6) 


Let  us  proceed  by  finding  and  integrating  it  over  the 

duct,  surface  to  obtain  Equations  (B.33)  and  (B.35)  for  and 

^Mn^D(p)  functions  of  and  qD/V^.  Finally,  Tj^  and 

^Hj,^D(p)  evaluated  in  Section  B.4. 

B.l  INTEGRATION  OF  TO  FIND  DUCT  FORCE  AND  MOMENT 

The  surface  velocity  distribution  on  the  duct  during  pure- 
pitch  flight  from  Equation  (7)  is 


W  "  w,  +  w.  +  w  +  w  + 

h  h^  a  q  q^ 


(B.7) 


where 


+  V 

^  2  3 


(B.8) 


and  the  components  are  given  by  the  following  equations: 


”h  * 

\ 

w 

a 

w 

“i 

w 

w 

w 

^3 

\ 

(A.l) 

(A.53) 

(A. 66) 

(A. 77) 

(A. 88) 

(A. 90) 

(A.  94) 

(A. 95) 

and  the  c's  do  not  depend  on  0  or  6. 


( 
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The  duct  normal  force  and  pitching  moment,  from  Equations  (B.2) 
and  (B.3)  are  given  by 

£T  2ir 


^D(P)  *  ^0^  j  6  dd  J  cos  0  <30 


(B.15) 


] 

] 


pcA 

D 


27r 


d0  (B.16) 


and 


^  STT  27r 

(**n)d(P)  "  ~  J  cos  0  sin  6  d6  J  W®  cos  0  d0 

o  o 


aTT 


(B.17) 


“  ■  /  ^2^®^  [^F  ^1^®^  ■*■  °i]  ®  ®  (B.18) 

o 

Note  that  terms  which  integrate  to  zero  have  been  omitted  in  the 
above  expressions.  We  shall  now  evaluate  the  remaining  integrals 
by  letting  2^-9  and  making  use  of  the  expressions  which  follow 
them. 

aTT 


[  “  /  c 

"  J 


f  sin  6  dd 
2 


■  c^  y  j^c^  (1  -  cos  d)  +  sin®  d  +  c^  sin  dj  dd 
o 

-  27rc.  (c^  + 


sir 


I  ■  /  c  f  cos  d  sin  d  dd  ■  -  TTC  c 
2  j  1  2  12 


(B.19) 


(B.20) 


i 


B-5 


ZTT 


•/ 


2^  ^1^2  ® 


(B.21) 


^3  "  ^  j^"  ®  '*’  2D  “  cos  9) 

'  0  *  tId)  **"  f]  -  cos  e)+  ^  sin^e  +  c^sin  ©j 


dO 
(B.22) 


"  /  I  ■  Td  sin^p  +  ^  sin®3  in  sin  p 

"  (^  ■*■  TiD)c'os  pj  +  iic)sin®p  -  qc  sin^p 

+  CgSin  p  cos  pj  dp 


(B.23) 


lr 


.  2:1 


dp 


+  §  in  ° 


IT 

^  I  [^°2  ^  qc)sin*P  -  qc  sin*pl 

O  J 


TT 

*  ^  J  [^°2  ■*■  ‘lo>sin‘‘p  -  qc  sin®p1  in  sin  p  dp 
ft  J 


l  “  ife)  /  OgSin^P 


dp 


(B.24) 
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27r 

^4  ~  ^  J  qc)  (sin^p  -  2  sin*p)  -  qc(sin*p  -  2  sin®p) 

o 

+  c^(sin  P  -  2  sin®p)cos  pj  dp  (B.27) 

TT  p 

^4  ‘  ^  J  I  ■  +  (§  ^  ®in  P 

o  ^ 

■  llD)l3i^]  [^°2  +  qc)sin2p  -  (2c^  +  3qc)sin‘‘p 

+  2qc  sin®p  +  c^(sin  P  -  2  sin®p)cos  pi  dp 


(B.28) 


B-7 


-  2 


/  • 


qc)p  sin®p  -  (2c  +  3qc)p  ain^p  +  2qc  3  sin®p 


+  c  O  sin  p  -  2p  sin®p)cos  p  dp 


] 

1  ^ 

■  2^^  i  +  5  ^  j  j^CCg+^IcJsin^P  -  (2c^+3qc)sin®p  +  2qc  sin®pj  dp 


2c 


IT 

J  |^(c^+qc)sin'*p  -  (2c^+3qc) sin®p  +  2qc  sin®pj  In  sin  P 


w 


-  (2  +  In  f  *  ^3  -  2  sin'^p^  dp 


-  *  «=)  [-7*  iId  («  ’)] 


+  7  )  +  (2c  +  3qc) 


*  ' — ^  5  In  j^gjj  -  g  )|  +  qc 


16D 
c 


(- 


“^(2  16D 


)  “  I  ^  ife  (• 


■*■  16D  16D  ■*■ 


(B.29) 
f  3t 

h 

^)i 

(B.30) 


I,  -  Y 


[2  "  2°D  16D  ■*■  3}]  “  3§D  6"  16D  il)"^  2^  (2 


+  In 


16D 


(B.31) 
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The  following  expressions  have  been  used  in  evaluating  the  above 


integrals: 


IT  TT 

f  p  Sin™  p  dp  -  ™  ^  J  p  sin™"®  p  dp 

>  o 


J  p  sin®  p  dp  -  ^ 


n 

J  P  sin"*  p  dp  -  ■— 


j  p  sin®  p  dp  - 


/  P  P  ■W  - 


p  cos  p  sin  p  dp 


P  cos  p  sin®  p  dp  -  - 


p  cos  p  sin®  p  dp  -  - 


‘ft  -'-.  .'<i..  -i ‘'f.-s'-  .'f 'ir'lRfe*''*  s!“  ^  SjWjWBI'^fiPKHipIP 


1 


I 


B-10 


TT 


I 


sin  3  In  sin  3  dp  » 


VTr(i4-i) 


TT 

f 


J 

o 


sin  p  In  sin  p  dp 


[.(!)-  ^<3, ] 


TT 

16 


(-3  in  4  *  i) 


TT 


/ 


sin®  p  In  sin  p  dp  ■ 


TT 


/  “ 


sin®  p  In  sin  p  dp  ■ 


■fe  (■  3  in  4  + 

jfg  (-  35  in  4  4  Iff) 


where 


(n  4  i). 


1  .  (2n  -  1) 


2n 


r  (n)  «  (n  -  1)  I  n 
f(n  +1)  ■  -  €  +  1  • 

■  -  €  -  in  4 
^(f)“  -  e  -  in  4  ■ 
V'(|-)--€-ln4- 


is  an  integer  >  0 


2 

3 


] 

] 

1 


I 

1 
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^(2)"  -  6  -  in  4  +  If 

V'(f)-  -  €  -  in  4  +f§| 

6  ■  Euler's  constant 

TT 

f  sin*"  P(ln  sin  p)cos  p  dp  ■  0  from  symmetry 


cos  0(1  -  cos  6)  ■  2  sin®  P  -  4  sin^  P 

sin®  6  cos  0  ■  4 (sin®  P  -  3  sin*  P  +  2  sin®  P) 
sin  0  cos  0  *  2 (sin  P  -  2  sin®  p)cos  p 


By  substitution  of  the  above  expressions  into  Equations  (B.16)  and 
(B.18)  we  obtain  finally  the  duct  normal  force 


+ 


_ Si 

COS  0 


(B.32) 


*»D(P) 


2pcAV^ 

D 


sin  -  (Trqc  /8VqD) 
1  +  (Trc/2D) 


+  cos 
o  o 


21SL 

8D 


+ 


4V_ 


+  w_ 


cos  <t> 
o 


16d) 


+  cos 
o  o 


ISL 

4D 


(B.33) 
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and  the  duct  moment  due  to  N_ 


(«n)i 


(C  -  TTC 

2  L  1 


yv_  _yc(, 

+  2  2D  v" 


-  L 

32D  V  12y  ®3  2  \2 


_S_+  1 

16D  3 


1  .  _c_ 

2  16D 


(B.34) 


(“n)i 


;  sin  -  (irqc  /8V^D) 

4D  ^  1  +  (Trc/2D)  [_  Y  V'°®  “o  4  y 

§  d"  tId  ■*■  j}]  16V^D  (’■”  il) 


w  +  w 
a.  q. 


-  Vj"coi"0"  )  V"  ife  I) 


(B.35) 


B.2  STATIC  (a)  COEFFICIENTS  OF  DUCT 

As  described  in  Section  2,  we  find  the  static  (a)  coef¬ 
ficients  of  the  duct  by  setting  q  ■  0  in  Equations  (B.33)  and 
(B.35).  This  gives 


(S)d(p)<“o>'’’°>  ■  “  § 


sin  a 


TT  -^  +  TT  cos  +  TT  ^ 
o  o 


^  6"  16D  *  0  -  16d)| 


(B.36) 
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C**n)d 


(P) 


(a„,0,0)  -  ^ 


Sin  a. 


1  +  (7rc/2Dy 


27r 


+  cos  a\  + 

O  Jo 


w 


Equation  (B.36)  is  identical  with  Equation  28  of  Reference  1  except 

for  the  addition  of  the  amount^ 7  -  2  I'O  2)  in  the  square  bracket. 

\16  8  /  ^''o 

Equation  (B.37)  is  identical  with  Equation  30  of  Reference  1  except 
for  the  sign  of  the  w^^  term  and  the  addition  of  the  amount 
(t  -  0  in  the  square  bracket.  All  of  these  differences 

appear  only  in  higher  order  terms  in  c/D  which  will  be  dropped 
later. 


B.3  DUCT  q  DERIVATIVES 

As  described  in  Section  2 ,  we  find  the  q  derivatives  of  the 
duct  by  differentiating  Equations  (B.33)  and  (B.35)  with  respect 
to  Dq/V^  holding  V^  and  constant.  After  substitution  of 

Equations  (A. 53)  and  (A. 77)  for  Wj^  and  w^  this  gives 


B~14 


Only  the  leading  terms  in  Equations  (B.38)  and  (B.39)  will  be  used 
in  the  subsequent  analysis.  These  are  of  lowest  order  in  c/D  and 
are  underlined. 


B.4  DUCT  FORCE  AND  MOMENT  DUE  TO  LEADING-EDGE  SUCTION 

The  thrust  force  Tjj(p)  moment  ^Mt}d(P)  duct 

arising  from  leading  suction  will  now  be  determined.  The  only 
parts  of  W  in  Equation  (B.7)  which  are  singular  at  the  leading 
edge  are  w^,  w^,  and  w  ,  which  are  proportional  to  tan  (0/2). 
Hence  by  analogy  with  the  flat-plate  formulas  (Eqs.  (A. 78)  and  (A. 81)) 
we  find  that  the  thrust  force  on  each  chordwise  strip  of  the  duct 
is  proportional  to  the  square  of  the  coefficient  of  the  singularity 
and  is  given  by 


JL 

2ir 
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(v, 


^sin 
o  o 


_  I3S. 

8D 
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cos  0 


1  + 


TTC 

2D 


T  ■  Tfpc 


(B.40) 
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APPENDIX  C 


VARIATION  OP  ACTUATOR-DISK  LOCATION 


r 

1. 

[ 

I 


To  determine  the  effect  of  actuator-disk  location  upon  the 
net  aerodynamic  reaction  on  the  duct,  we  shall  calculate  the  sur¬ 
face  vorticity  distribution  during  hovering  flight  for  two 

alternate  disk  locations,  the  duct  leading  edge  and  roidchord. 

It  is  recalled  that  the  actuator  disk  was  assvimed  to  be  at  the 
duct  trailing  edge  in  Section  3. 

The  leading-edge  location  is  considered  first.  In  this 
case  it  can  be  seen  (Sketch  A.l)  that  Equation  (A. 14)  is  replaced 
by 


X 


s 


(C.l) 


By  following  through  the  analysis  in  Section  A.l  with  this  modifi¬ 
cation  we  find  that  the  ec[uations  therein  are  changed  slightly  as 
indicated  below. 

Equations  (A. 15)  and  (A. 16)  become 


(cos  0+1) 


) 


(C.2) 


(C.3) 


C-2 


Substitution  of  p  ■  0  +  tt  into  the  identity  used  to  obtain 
Equation  (A. 26)  gives 

00 

ln(l  +  cos  p)--  In  2-2^  ^  (-1)"  cos  np  (C.4) 

n«i 

By  use  of  this,  Equation  (A. 16)  can  be  written 

00 

“  it  (2  -  if^)sin  9  -  ^  Y.  n  ® 

®  '  n-i 

and  Equation  (A. 26)  becomes 


(C.6) 


Thus  the  factor  (-1)*^^  is  inserted  into  the  summation 
symbol  in  Equations  (A. 32)  and  (A.33) .  Using  0  ■  p  +  tt  in  the 
identity  from  Reference  22, we  obtain 

f  "P  .  .  I  (C. 

n"i 


C-3 


when 


-TT  ^  <  TT 


Multiplying  by  (cos  3  +  1),  reducing  and  re-indexing  gives  the 
identity 


E 


n*2 


n+i 


n(n^  -  1) 


sin  p  -  (cos  p  +  1)  I  (C.8) 


when 


-TT  <  P  <  TT 


By  use  of  this,  the  main  equations  in  the  series,  Equations  (A. 33) 
through  (A. 47)  become 


‘*’s'  "  “  2^  ^  L  y  6  +  1)  +  •!■  sin  ej  +  a^  sin  9  (C.9) 

"f  '  ■  -  ^  S  [-  ^  2°°  ^  +  I  sin  0  t  I  COB  e 

s  ^ 


+  in  -i—  -  I")  cos  ej  (C.IO) 

“Yt  ■  *  S 

(2  -  l"  f )  (=-13) 
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^  ^  (in  i|2  -  3)  (C.14) 

^  f  *  ®) 

Finally  Equation  (A. 50)  for  the  bound  vorticity  distribution  along 
the  duct  chord  becomes  (with  the  actuator  disk  at  the  duct  leading 
edge) 

■  T  [C'"  -  2)  in"  f  -  «] 

where 

0  <  0  <  TT 

Now  following  through  the  analysis  in  Section  A.l  with  the 
disk  at  the  midchord  of  the  duct,  we  find  that  the  equations  in 

Section  A.l  are  changed  as  follows.  Equation  (A. 14)  becomes 

Xg  -  (C.17) 

With  this  modification  Equations  (A. 15)  and  (A. 16)  become 

“I's  "  ^  [5  S  '‘I®  ®  (i  e  -  in  (C.18) 

''y,  ■  (2  1  7  in  ■=“*  e  -  in  “)  (0.19) 

Using  the  identity  from  page  188,  Reference  22, 

00  ^ 

-  ^  cos  2n  ■  j  In (4  cos^d)  ■  In  2  +  j  In  coa^d  (C,20) 

n-i 
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we  can  write 


V  sm  0  •  y 

27r 


COS  2n0  sin  9  (C.21) 


n«i 


j  (c.: 


for  even  n  ■  0 

a  m  .yS.(l  _  16^ 

1  47r  \  2  ~  J 


YC ( “1 ) ^ 

®n  “  -  472n(n  -  l)  "  >  1 


)  (C.22) 
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n 

27r  "  47r 


oo 

^  (-J-)"sin(2n  -1)0  .  /3  ,  16D\  .  ^ 

^  2n(n  -1)  ^  (^2  “  ® 

n“2  ' 


(C.23) 


Using  the  series  from  page  188,  Reference  22, 


Y'  —  sin  np  ■  P 
^  n  ^  2 

n“i 

-for 


0  <  8  <  27r 

and  letting  p  -  e  +  7^,  we  obtain 


(-1)”  sin  2n0  .  _  g 

n-i 


(C.24) 


(C.25) 


for 


C-6 


Forming  the  product  -9  cos  6  and  simplifying  gives  for 
-  Tr/2  <  0  <  +  7r/2 


20  cos  0  -  sin  0  -  ~ 


(C.26) 


n-a 


By  use  of  this,  the  main  equations  in  the  series,  (A. 33)  through 
(A. 47)  become 


<D  .  -  _  JL.  _  1£ 
s  iTT  4Tr 


0  cos  0  +  ^  -  In  sin  9 


and 


,  .  -  _  Jl  +  OS 

F  2ir  4'jr 
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i^ln  -  2^  cos  0+0  sin  0 
when  -  Tr/2  <  0  <  +  ir/2,  and  when  ir/2  <  0  <  37r/2 


r.  0  • 

<1,  .  -  _  Jl.  _ 

s  2ir  4ir 


(0  -  7r)cos  0  +  3in  0 


and 


-  27  47  ^  -  2  +  <0  -  ^)sin  ©j 


We  also  find 


(C.27) 


(C.28) 


(C.29) 


(C.30) 
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27r  ^  47r  v"  c  y 

(C.31) 
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(C.32) 
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Finally  Equation  (A. 50)  for  the  bound  vorticity  distribution  along 
the  duct  chord  becomes  (with  the  actuator  disk  at  the  duct  mid¬ 
chord) 

■''h  ■  ?  [0-  ^  -  2)  tan  I  -  e] 

for 

and  (C.35) 

for  f  <  0  <  ^ 

It  can  be  seen  from  Equations  (A.50) ,  (C.16),  and  (C.35)  for  the 
three  alternate  locations  of  the  actuator  disk  that  the  net  vor¬ 
ticity  distribution  7  Yjj  invariant  with  disk  location. 

That  is,  as  the  disk  moves  forward  in  the  duct,  its  trailing  vor¬ 
ticity  y  is  subtracted  from  so  that  the  sum  7  + 

invariant.  Thus,  it  becomes  evident  that  with  the  present  flow 
model  the  flow  pattern  (and  net  duct  reaction)  is  independent  of 
the  actuator -disk  location  for  all  flight  conditions  (not  only 
hovering)  so  long  as  the  disk  is  located  within  the  duct.  The 
static  pressure  rise  through  the  actuator  disk  merely  causes  hoop 
tension  in  the  duct  downstream  of  the  disk  and  does  not  directly 
affect  the  resultant  force  on  the  duct  in  the  present  analysis 
since  its  diauneter  is  constant.  For  diverging  ducts  this  pressure 
rise  does  contribute  directly  to  duct  thrust,  however. 


APPENDIX  D 


INVESTIGATION  OF  AVERAGING  PROCESSES 

Two  averaging  processes  have  been  used  in  Reference  1  and  In 

Appendix  A  herein  to  obtain  the  axial  velocity  distribution  on  the 

duct  surface  (w.  ■  w,  +  w  +  )  which  is  self-induced  by  the 

i  h^ 

bound  vorticity  distribution  on  the  duct  (7,  +  y  +  •  Specif i- 

*1  ct  ^ 

cally:  (1)  the  bound  vorticity  (7)  was  first  averaged  over  the 

chord  length  to  give  7  and  (2)  the  surface  velocity  distribution 
induced  by  7  was  averaged  over  the  chord  length  to  give  w^.  ^ 

To  assess  the  inaccuracy  which  is  introduced  into  the  final 
results  of  the  analysis  by  these  approximations,  we  shall  use  the 
actual  chordwise  distribution  of  bound  vorticity  7  to  calculate 
the  actual  value  of  the  self-induced  surface  velocity  w^  at  the 
duct  • 


D.l.  HOVERING  COMPONENT,  W^ 

^i 

The  bound  vorticity  on  the  duct  during  hovering  (or  axial) 
flight  is  (Eq.  (A. 50)) 

We  shall  now  find  w^  the  velocity  induced  by  7^  at  three 
points  along  the  chord  and  compare  these  three  values  with  the 
average  value  w^^  which  was  used  in  Section  3 . 


^Only  the  axial  velocity  component  was  considered  because  the  rad¬ 
ial  component  induced  by  a  vortex  ring  for  small  c/D  is  nearly 
the  same  as  that  induced  by  a  straight  line  fileunent  (D  ■  ») 
Hence  the  radial  component  of  the  self-induced  velocity  is  zero. 


D-2 


(cos  6  -  cos  0') 


Sketch  D.I.-  Duct  nomenclature. 

From  Equation  A. 48  (with  ■  c/2  (cos  6  -  cos  0')  from 
Sketch  D.l)  we  find  that  each  elementary  vortex  ring 
induces  at  the  axial  velocity  component 


au. 

2TrD 


1  n 

8D  1 

1 

in 

c(cos  6  -  cos 

0') 

-  ± 

(D.2) 


The  total  induced  axial  velocity  at  x^  is  therefore  obtained  by 
integration  over  the  duct.  That  is, 
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2  ^s 


w, 


^i  27r^D 


-  2^  tan  1^  +  TT  -  0'  ^In 


8D 


-  1 


w, 


-  (f  -  -s) 

-  In  j  cos  6  -  cos  »i)«R 

j  jjO”  ~  0'^  +  TT  Sin  0' 


(D.3) 


yc 


-  0'  sin  0'j  ^In  -  1  -  In  |  cos  0  -  cos  0'|  )  d0'  (D.4) 


8D 


i) 
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Evaluating  this  expression  at  three  locations  we  obtain: 
at  midchord  (2x_/c  ■  cos  S  ■  0) 


w  -  ^ 
h^  4VD 


(in  ^  -  l)  (in  +  In  2  (in  ^  -  2)  +  1 


(D.5) 


at  the  trailing  edoe  (2x  /c  ■  cos  6  ■  1) 
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"h^  -  ^  -  0  ('"  “  -  0  *  0"  ^  2  - 1)  *  v"  2] 


(D.6) 


and  at  the  leading  edge  (2x  /c  ■  cos  0  “  -  1) 
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V  •  ™  [('"  0  *  (‘"  -  j)  0"  2  •  1) 


+  2  -  3  In  2 


Now  from  Equation  (A. 53),  we  have  for  the  averaged  value 


(D.7) 


w  -  ^ 
hj^  4liD 


(D.8) 


Note  that  the  first  term  in  the  bracket  of  Equations  D.5  through 
D.8  is  identical,  and  for  very  small  values  of  c/D,  the  remainder 
of  the  bracket  is  negligible.  However,  the  actual  value  of  the 
remainder  can  be  expressed  in  the  form 

A  (in  ^  -  1)  +  B  (D.9) 

in  each  case  with  the  value  of  A  and  B  tabulated  as  follows. 
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Chordwise  location 

A 

B 

Leading  edge 

In  2  +  1  -  1.69 

1  -  4  In  2  -  -1.77 

Midchord 

In  2  =  0.69 

1  -  In  2  -  0.31 

Trailing  edge 

In  2  -  1  -  -0.31 

1 

Averaged  value 

used 

3/2  -  In  2  -  0.81 

Evaluation  of  Equations (D. 5)  through(D.8)  for  two  values  of  c/D  gives 

the  values  for  4t7D  w,  /cv  tabulated  below 

^i 


Chordwise  location 

nmoQ^nii 

Leading  edge 

13.09 

18.92 

Midchord 

11.87 

16.92 

Trailing  edge 

9.08 

13.54 

Averaged  value 
used 

11.86 

17.09 

Thus ,  the  averaging  process  introduces  an  error  in  the  local  value 
of  Wj^  of  about  30  percent  for  c  “  0.2D. 


D.2  ANGLE -OF- ATTACK  COMPONENT,  W 


a. 


From  Equation  (A. 63,  we  have  for -the  bound  vo^ticity  on  the 
duct  induced  by  the  free  stream  at  angle  of  attack 


'a 


Wa  (O')  -  (2Tr  -  e-)  -  tan  ^ 


a 


(D.IO) 


D-5 


Now  sxibstltuting  for  dP  and  using  Co  ■  c/2  (cos  0  -  cos  &•)  in 

Ob  X\ 

Equation  (A. 73)  we  obtain  the  total  induced  axial  velocity  at 
IT 


(1  -  cos  6' ) 


0 


In 


cos 


$  -  cos  6 


(D.ll) 


Evaluation  of  this  integral  gives: 


2TrD 

r 


a 


(D.12) 


where  *  2,3,4  at  the  duct  leading  edge,  mldchord,  and  trailing 
edge,  respectively. 

Also  from  Equation  (A.77) ,  we  can  write  the  averaged  value  of 
as 


w  -  In  —  -  I  -  In  4  -  In  -  2.89  (D.13) 

Oi  c  2  c 

This  value  is  nearly  equal  to  the  value  at  mldchord.  The  following 
values  of  (2TrD/r  )w  are  obtained  from  Equations  (D .  12)  and  (D .  13) 

0^  01  jL 

for  two  values  of  c/D. 


Chordwise  location 

c  -  0.2D 

c  -  O.ID 

Leading  edge 

2.39 

3.07 

Mldchord 

1.39 

2.07 

Trailing  edge 

0.39 

1.07 

Averaged  value 
used 

1.50 

2.18 
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It  can  be  seen  that  the  actual  values  depart  considerably  from  the 
averaged  value  which  was  used. 

D.3.  PITCHING  COMPONENT,  W 


The  bound  vorticity  on  the  duct  7^  induced  by  the  pitching 
motion  q  (from  Eqs.  (A. 88)  and  (A. 94))  is 


Yq  “  7q 


+  Y  ■  qc  cos  0  (  sin  0*  -  tan  )  (D-14) 


Using  Equation  (A. 73)  (with  dP  “  (cos  9  -  cos  6’)) 

we  obtain  the  total  induced  axial  velocity  at  x 


TT 
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w  “  cos 
q.  4TfD 


sin^e*  -  (1  -  cos  0') 


1  + 


2D 


(‘■'f  - 


-  In  I  cos  9  -  cos  ®'|'^  99' 
Evaluation  of  this  integral  gives: 
at  midchord  (2x  /c  “  cos  0*0) 


(D.15) 
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acf 
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(in  f  -  3)  +  in  2  +  i  -  (in  if  -  3)]  cos  ^ 

2D 


(D.16) 


at  the  trailing  edge  (2x  /c  “  cos  0*1) 


w  -  ^ 
8D 
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(D.17) 
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and  at  the  igadinq  {2x  /c  -  cos  ©  -  -  1) 


1  ^0/20 
"  2  " 


6"  ^  -  0 

^  2D 


cos  0 


(D.18) 

Now  from  Equation  (A. 95) ,  we  can  write  for  the  averaged  value 


“q,  ■  ^  77^  [O'  f  -  3)  f  -  in  2I  cos  « 

^  2D 


(D.19) 


By  use  of 

I  — ^ -  jw  a  (D.20) 

\qc®  cos  0/  ‘^i  ® 

evaluation  of  Equations  (D.16)  through  (D.20}  for  two  values  of  c/D 
gives  the  following  values  of 


Chordwise  location 

mam 

mopoiiij 

Leading  edge 

0.31 

1.16 

Midchord 

1.55 

2.30 

Trailing  edge 

0.79 

1.44 

Averaged  value 
used 

1.14 

1.90 

It  is  found  that  the  components  w^  vary  as  much  as  100  per¬ 
cent  (hut  not  by  an  order  of  magnitude)  from  the  average  values 
w^  previously  obtained.  The  self-induced  velocity  components  gen¬ 
erally  have  only  small  effects  upon  most  of  the  final  stability 


D-8 


derivatives.  All  of  them  affect  only  higher-order  terms  in  c/D 

for  the  pitching  derivatives.  The  a  coefficients 

are  unaffected  because  neither  w.  nor  vr  affect  the 
D(P)  h^ 

leading-edge  singularity.  The  error,  from  averaging  w.  and 

”i 

w  ,  in  Equation  (B.14)  for  N-  may  be  appreciable,  however, 
ct  ^  ^ 

Inasmuch  as  the  derivative  (  C„  \  is  found  to  dominate  the 

V  ®*a/D(P) 

dynamic  stability  of  a  torpedo-like  configuration,  this  averaging 
error  would  seem  to  warrant  further  attention. 


APPENDIX  B 


HULL-BOUNDARY-LAYER  CALCULATIONS 


In  order  to  determine  the  center  of  pressure  of  the  hull  at 
angle  of  attack  (a/L  in  Section  3.5.1  and  Fig.  5),  the  volume  of 
the  hull  combined  with  the  displacement  thickness  of  its  boundary 
layer  is  required,  as  seen  from  Equation  (76) .  This  volume  will 
now  be  determined  by  using  the  "effective"  hull  radius.  Thus,  we 
let 


where  r^  is  the  actual  hull  radius.  Now,  for  the  assumed  flat- 
plate  turbulent  boundary  layer,  we  have  for  the  displacement  thick¬ 
ness  5* 

6*  -  0.37  (Ref^/®  (2  "  *h) 


and 


Thus,  we  find  from  Equations  (E.l),  (E.2),  and  (78)  that 


(E.3) 


(E.4) 
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The  volume  of  the  hull  and  boundary  layer  is  then  given  by 


Vol*  -  Trr„® 


-1 


(E.5) 


Thus 


,  letting  i  =  2  obtain 
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"l  15  Vl4  19  j  r^  13  J 
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for  Re.  ■  4X10^  we  have  from  Equation  (86) 


-  ^  m  fl  *  ^  m  0  223 
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(E.7) 
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Vol* 

Trr 
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0.523  +  0.167  +  0.019  -  0.709 


(E.9) 


In  addition  to  the  volume  the  following  integral  is  needed 
in  Section  3.5.1  for  use  in  Equation  (95). 
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-  0.223 
^m 
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The  radius  of  gyration  of  the  hull  will  now  be  determined 

for  use  in  Section  3.5.1.  The  mass  of  the  hull  is  given  by 


mjj  -  p  Vol  -  ^  Trr/Lp 

and  about  the  hull  centroid,  at  Xjj  ■  0,  we  have 
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(E.14) 
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(E.20) 

Now,  about  the  hull  center  of  gravity,  *jj  “  W^O, 

we  find  for 

the  radius  of  gyration 

V"  ‘  ^  ^  So  -  0-0457  L- 

(E.21) 

Ky"  0.214  L  -  2.14  ft 

(E.22) 

TABLE  IV.-  COMPARISON  OF  DUCT  AND  PROPELLER  STABILITY 
DERIVATIVES  DUE  TO  PITCH. 


Stab. 

Deriv. 

D{p) 

P(D) 

c  * 

9.91 

< 

0.8 

(due 

to 

a^') 

q 

c  ' 

0 

< 

.15 

(due 

to 

q 

c 

m 

.364 

.065 

(due 

to 

“o’) 

q 

< 

.5 

(due 

to 

q) 

c  * 

0 

< 

2 

(due 

to 

“o’) 

q 

c  * 

""q 

-98.7 

1 

< 

8 

TABLE  V.-  STABILITY  DERIVATIVES  ABOUT  HULL  CENTROID  (0.5L) 
FOR  THE  CONFIGURATION  OF  FIGURE  4. 


Flight 

Stability 

Hull 

Duct 

Propeller 

condition 

derivative 

H 

D(P) 

P(D) 

"^N 

0.099 

1.95 

<  0.15 

a 

a 

C  ' 

"a 

28.0 

-19.4 

1  1  <  1.5 

C  *  '  C  ' 

n*  '  m 

ct  cc 

c  ' 

.495 

9.91 

<  .8 

q 

-Cy  •  <  .15 

q 

q 

52.6 

-98.7 

<  8 

^n*  ’  «  ’ 

q  a 

14.2 

.126 

Small  compared 

with  duct 

• 

a 

C  ' 

-57.6 

-1.26 

Small  compared 

•"d 

with  duct 

q  +  d 

C  '  +  C  ' 

-5.0 

-100.0 

Small  compared 

_ 3 _ _ 

with  duct 

TABLE  VI . -  STABILITY  DERIVATIVES  ABOUT  THE  VEHICLE  CENTER  OF 
GRAVITY  (0.4L)  FOR  THE  COMPLETE  CONFIGURATION  OF  FIGURE  4. 


Flight 

condition 

stability 

derivative 

Hull  . 

H 

Ducted 

propeller 

D(P) 

Fin 

F 

Total 

HDPF 

'“a 

0.099 

1.95 

0.86 

2.91 

a 

C  " 

“a 

27.8 

-23.3 

-10.3 

-5.8 

C  " 

.594 

11.9 

5.15 

17.6 

q 

"q 

79.4 

-141.8 

-61.9 

-124.3 

^N. 

a 

14.2 

.126 

.038 

14.4 

a 

"‘d 

-86.0 

-1.51 

-.45 

-88.0 

q  +  a 

c  +  c 

m. 

q  a 

*  -6.6 

-143.3 

-61.5 

-212.3 

TABLE  VII.-  COMPARISON  OF  PRESENT  THEORY  WITH  EXPERIMENTAL 
DATA  BY  GILL,  REFERENCE  13 ,  FOR  AXIAL  FLOW. 


V. 

-J. 

'^o 

^T 

^P(D) 

exper . 

JL 

V 

o 

^D(P) 

^P(D) 

exper. 

^D(P)  _  0.78 

^P 

exper. 
Eq. (205) 

2 

’'P(D) 

(pred. ) 

Eq.  (203) 

’ip  (V./V^)+l 

(pred. ) 

Eq. (206) 

4.1 

IG 

3.1 

0.25 

0.47 

0.25  ’ 

0.39 

2.4 

o 

• 

in 

1.4 

.12 

.33 

.44 

.59 

1.8 

2.2 

0.8 

.16 

.22 

.51 

.71 

TABLE  VIII.-  COMPARISON  OF  PRESENT  THEORY  WITH  EXPERIMENTAL 
DATA  BY  HORN,  REFERENCE  17,  FOR  AXIAL  FLOW. 


V 

o 

nD 

^T 

P(D) 

exper . 

V 

o 

’’P(D) 
exper . 

V(P)  -  0.7 

^P 

exper • 

n  “  - _ ^  ^ 

^P(D) 

(pred . ) 

Eq. (203) 

(Vj/V^j+I 

(pred.) 

Eq.  (206) 

0 

0.46 

0.25 

1.14 

0.46 

.32 

0.125 

0.31 

0.81 

.50 

.22 

.11 

.18 

.036 

.51 

.95 

.75 

.068 

.034 

0 

.012 

.59 

,97 

1.0 

.018 

.009 
_ 1 

-.57 

.006 

.27 

1.0 

TABLE  IX.-  COMPARISON  OF  PRESENT  THEORY  WITH  EXPERIMENTAL 
DATA  BY  ALLEN,  REFERENCE  21,  FOR  AXIAL  FLOW. 


^T 

'^^P  (D) 
exper. 

'^D(P)  0.7 

^  -  2 

nD 

^P(D) 

exper. 

*P(D) 

(pred.) 

Eq.  (203) 

^P 

exper . 

^P  (Vj/V^)+1 

(pred.) 

Eq.  (206) 

20 

0.40 

0.92 

0.38 

0.087 

0.11 

0.70 

0.84 

25 

.60 

.40 

.19 

.089 

.061 

.77 

.91 

TABLE  X.-  COMPARISON  OF  PRESENT  THEORY  WITH  EXPERIMENTAL 
DATA  BY  KRUGER,  REFERENCE  19,  FOR  AXIAL  FLOW. 


V 

o 

OM 

P(D) 

exper . 

JL 

V 

o 

"^PCD) 

exper. 

"^DCP)  ^  0.75 

^P 

exper . 

^  -  2 

*P(D) 

(pred.) 

Eq. (203) 

(Vj/V^)  +  1 

(pred. ) 

Eq.  (206) 

15 

0.14 

3 

3 

0.3 

0.45 

0.42 

0.50 

25 

.19 

2 

0.73 

.5 

.25 

.43 

.73 

25 

.26 

1 

.41 

.6 

.13 

o 

in 

• 

.83 

25 

.32 

0.5 

.22 

.6 

.08 

.56 

.90 

35 

00 

• 

.22 

.10 

.8 

.04 

.58 

.95 

TABLE  XI.-  COMPARISON  OF  PRESENT  THEORY  WITH  EXPERIMENTAL 
DATA  BY  GROSE,  REFERENCE  14,  FOR  AXIAL  FLOW. 


V 

o 

nD 

^T 

^P(D) 

exper. 

JL 

^D(P) 
'’^P(D) 
exper • 

^D(P)  ...  0.70 

^P 

exper . 

2 

■fpCD)  ^  1+2(V</Y) 

(pred.) 

Eq. (203) 

"Ip  (Vj/V^)+1 

(pred. ) 

Eq. (206) 

22 

0.72 

0.64 

0.28 

0.044 

0.084 

0.77 

0.88 

22 

.63 

.97 

.39 

.086 

.115 

.70 

.84 

27 

.72 

1.09 

.45 

.095 

.127 

.70 

.82 

TABLE  XIZ.-  COMPARISON  OF  PRESENT  THEORY  WITH  EXPERIMENTAL 
DATA  BY  GROSE,  REFERENCE  14,  AT  ANGLE  OF  ATTACK. 


V 

o 

nD 

^T 

^P(D) 
exper . 

JL 

V 

o 

(^Od(p) 

(per  rad) 

exper.  - 

C^Od(p) 

(per  rad) 

(pred.) 

Eq.  (TOar 

22 

0.25 

0.12 

3.58 

3.60 

32 

H 

.16 

.08 

3.56 

3.54 

TABLE  XIII.'  COMPARISON  OF  PRESENT  THEORY  WITH  EXPERIMENTAL 
DATA  BY  GILL,  REFERENCE  13,  AT  ANGLE  OF  ATTACK. 


“o 

V. 

V 

o 

(=» 

)d(p) 

■Ml 

MM 

)d(p) 

exper. 

pred. 

Dq. (209) 

exper. 

pred. 

Eq.  (210) 

exper. 

pred. 

Eq. (211) 

50 

3.40 

18 

1.76 

2.02 

0.11 

0.56 

0.42 

2.52 

50 

1.80 

9 

.88 

1.32 

.20 

.75 

.23 

1.13 

80 

3.25 

15 

1.12 

2.12 

.15 

.82 

.44 

3.68 

80  j 

2.05 

— j 

9 

1.07 

1.47 

.10 

1.12 

.30 

2.41 

(a)  Ducted  actuator  disk  nomenclature. 


On  inner  surface  0  <  6  <  ^ 

On  outer  surface  tt  <  6'  <  2^ 

X  ■  c/2  cos  6 

■  S 

(b)  ^  plane  transformation. 

Figure  1.-  Mathematical  representation  of  ducted  propeller. 


Figure  2.-  Nomenclature  for  ducted  actuator  disk 
in  pure  pitch  flight. 


ir  of  gravity  at  0.4L 
it  0.5L 

^Ducted  propeller 
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■  D 
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Duct 
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swirl 


thrust  ratio 
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►-1.5" 
(propeller 
spacing) 

'■  '  s 

h 

f 
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Angles  indicated  are  blade  pitch  setting  p. 


(a)  Data  by  Gill,  Reference  13. 

Figure  7.-  Comparison  of  theoretical  and  experimental 
thrust  ratio  in  hovering  flight. 


thrust  ratio 


c/D  chord  to  diameter  ratio 


(b)  Data  by  Horn,  Moser,  and  Kruger. 
Figure  7.-  Continued. 
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(c)  Data  by  Grose  and  Platt, 
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